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FOREWORD 


This  is  the  final  technical  report  on  Contract 
F19628-69-C-0277,  Project  8693.  The  work  reported 
herein  was  accomplished  by  Hughes  Research  Labora¬ 
tories,  Malibu,  California.  This  report  covers  the 
study  period  16  June  1969  through  15  June  1972. 

Principal  investigator  for  the  program  was 
C.  R.  Giuliano.  L.  D.  Hess  and  G.  R.  Rickel  contributed 
to  the  experimental  phase.  Principal  contributor  to  the 
theoretical  aspect  of  the  program  was  R.  W.  Hellwarth, 
with  D.  F.  DuBois  providing  additional  input. 

This  document  was  submitted  by  the  authors 
September  1972.  The  Air  Force  contract  monitor  was 
Erlan  S.  Bliss,  Capt.  USAF. 
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ABSTRACT 


This  report  contains  a  cumulative  summation  of  a  three-year  experi¬ 
mental  and  theoretical  program  on  laser  induced  damage  in  ruby  and 
sapphire  at  6943  A.  Both  bulk  and  surface  damage  were  studied,  using 
a  well-characterized  single  transverse  and  longitudinal  mode  Q-switched 
ruby  laser  and  amplifier.  Temporal  studies  of  the  development  of 
internal  damage  tracks  using  a  high  speed  streak  camera  have  been 
correlated  with  self-focusing  theory  with  good  agreement  and  constitute 
the  first  direct  evidence  for  moving  self-foci  in  solids.  It  has  been 
found  that  beams  with  elliptical  spatial  cross-sections  have  higher 
self-focusing  thresholds  than  circular  beams.  This  result  indicates 
that  problems  that  arise  from  self-focusing  of  high  power  optical  beams 
can  be  avoided  by  proper  choice  of  beam  shape.  An  investigation  into 
the  relation  between  surface  damage  and  surface  plasmas  is  described. 
The  study  of  temporal  and  spatial  surface  plasma  differentiation  shows 
that  the  entrance  plasma  consists  of  two  components:  an  air  plasma 
and  an  explosion  plasma,  while  the  exit  surface  plasma  is  of  the  explo¬ 
sion  type.  The  surface  plasmas  are  shown  to  be  a  result  rather  than  a 
cause  of  the  surface  damage.  Surface  damage  experiments  on  samples 
polished  by  energetic  Ar+  ion  beams  have  shown  a  marked  increase 
(up  to  lOx)  in  damage  threshold  indicating  that  surface  finish  and/or 
crystallinity  are  very  important  in  determining  damage  resistance. 
Theoretical  efforts  have  developed  a  new  model  for  optical  breakdown 
in  materials  based  on  a  transfer  of  optical  energy  to  photoelectrons 
which  undergo  rapid  energy  transfer  to  acoustic  phonons  via  optical 
phonons.  Also,  a  theoretical  framework  has  been  developed  describing 
the  momentum  distribution  of  conduction  electrons  in  a  polar  lattice 
driven  by  strong  optical  fields.  Models  of  self-focusing  mechanisms 
are  discussed,  and  nondamaging  methods  are  suggested  for  measuring 
surface  absorption  which  can  result  in  distinguishing  between  surface 
damage  mechanisms. 
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SECTION  I 


EXPERIMENTAL  STUDIES: 

OPTICAL  DAMAGE 

A.  INTRODUCTION  AND  SUMMARY 

This  report  represents  the  culmination  of  a  three  year  program 
on  laser  induced  damage  in  ruby  and  sapphire.  The  broad  goals  of  the 
program  were  to  perform  experiments  designed  to  elucidate  the  nature 
of  the  damage  mechanism  or  mechanisms  and  to  provide  means  of 
avoiding  or  eliminating  damage.  A  parallel  theoretical  effort  was  to  be 
carried  out,  the  function  of  which  was  to  examine  present  theories  of 
damage  and  to  explore  new  theories  for  damage  mechanisms  in  the  light 
of  the  experimental  observations. 

The  main  experimental  emphasis  in  this  program  has  been 
placed  on  bulk  damage.  In  the  early  stages,  bulk  damage  thresholds 
were  measured  on  a  number  of  samples  from  different  sources.  There 
has  been  much  indirect  evidence  that  bulk  damage  in  inclusion-free 
solids  is  initiated  by  self-focusing.  We  have  studied  the  temporal  evo¬ 
lution  of  damage  tracks  in  sapphire  and  ruby  using  a  high  speed  streak 
camera  and  correlated  the  damage  evolution  with  the  theory  for  a 
moving  catastrophic  self-focus.  These  results,  which  are  presented  in 
Section  I-B-l  through  4,  constitute  the  first  direct  evidence  for  a 
moving  self-focus  in  solids  and  give  direct  confirmation  that  the  ini¬ 
tiating  mechanism  for  bulk  damage  is  self-focusing.  Further  results 
on  damage  threshold  for  different  focusing  conditions  presented  in 
I-B-5  indicate  as  do  the  above  mentioned  results  that  both  electro- 
strictive  and  electronic  self-focusing  mechanisms  are  operative  for  the 
pulse  durations  and  spot  sizes  employed  in  these  studies.  Damage 
thresholds  as  a  function  of  optical  pumping  in  ruby  samples  are  dis¬ 
cussed  in  Section  l-B-9. 
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Surface  damage  threshold  studies  undertaken  in  the  latter  stages 
of  this  program  have  attempted  to  provide  an  understanding  of  the  rela¬ 
tion  between  surface  damage  and  surface  plasmas  (Section  I-C-l)  and 
have  shown  that  the  plasmas  are  a  result  of  rather  than  the  cause  of 
surface  damage;  although  some  features  of  surface  damage  (especially 
on  the  entrance  surface)  are  shown  to  result  from  the  interaction  of  the 
plasma  with  the  surface  (Section  I-C-3).  Also  contained  in  Sec¬ 
tions  I-C-l  and  2  are  the  result  of  surface  damage  threshold  measure¬ 
ments  as  a  function  of  beam  size  and  divergence.  These  results  show 
that  both  entrance  and  exit  damage  mechanisms  are  power  density  (or 
energy  density)  dependent  and  are  not  complicated  by  nonlinear  propa¬ 
gation  phenomena  (except  for  ruby  exit  damage  which  is  discussed  in 
Section  I-C-2).  Surface  damage  morphology  is  also  discussed  in 
Section  I-C,  and  differences  between  entrance  and  exit  damage  are 
illustrated  by  the  use  of  scanning  electron  micrographs. 

In  the  interest  of  avoiding  or  eliminating  damage,  we  have  dis¬ 
covered  two  very  promising  areas  for  pursuit.  First,  for  bulk  damage, 
it  has  been  found  during  this  work  that  the  threshold  for  self-focusing 
is  greater  for  beams  having  an  elliptical  spatial  cross-section  than  for 
circular  beams.  This  effect  is  discussed  in  Section  I-B-6.  For  unfo¬ 
cused  elliptical  beams,  the  self-focusing  threshold  increases  linearly 
as  the  ratio  of  major  to  minor  elliptic  axes.  Hence,  one  can  avoid 
problems  that  arise  from  self-focusing  by  appropriately  choosing  the 
beam  shape.  The  second  promising  area  involves  the  increase  of  sur¬ 
face  damage  thresholds  with  ion  beam  polishing.  At  this  time,  these 
threshold  increases  of  up  to  lOx  and  greater  have  been  realized,  but  the 
nature  of  the  improvement  is  not  understood  in  detail.  The  phenomenon 
can  be  explained  in  terms  of  a  distinct  improvement  in  surface  finish 
or  crystallinity  or  a  combined  influence  of  both.  The  results  of  these 
observations  are  discussed  in  Section  I-C-4. 

The  theoretical  effort  in  this  program  has  given  interpretation 
of  experimental  results,  developed  new  theories  for  damage  mecha¬ 
nisms,  and  suggested  new  experiments  to  further  elucidate  possible 
mechanisms  for  surface  damage.  Early  in  the  program  a  new  model 
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was,  developed  for  the  optical  breakdown  of  a  material  based  on  a,trans- 

1  i  i  ‘  _ 

fer  of  the  optical  energy  to  .photoelectrons  which  undergo  rapid  transfer 
to  acoustic  phonons  via, optical  phonons.  This  treatment  is  presented  in 
Appendix  A. ,  In  Appendix  B  we  present  the  development  of  a  theoretical 

I 

framework  for  the  study  of  the  rpomentum  distribution  of  conduction 

■,  |  I 

electrons  in  a  polar  lattice  being  driven  by  a  strong  optical  field. 

Although  we  have  not  expended  the  fconsiderable  computational  effort 

i 

necessary  to  obtain  numerica4  estimates  of  numbers  of  ihigl^  energy 
electrons,  we  hav,e  been  able  to  show  that  many  ad  hoc  hypotheses  that 
have  been  put  forward  are  without  theoretical  basis. 

In  Section  1I-A  we  discuss  the  development  of  models  of  the  i 
mechanisms  responsible  for  the  self-focusing  that  triggers  bulk  damage. 

We  have  found  evidence, that  electronic  hyperpolarizability  is  not  the 

i  .  ' 

dominant  mechanism  in  laser  glass,  sapphire,  and  YAG  but  probably  is 
in  the  alkali  halides.  We  have  pointed  out  ways  of  establishing  the 

i 

mechanisms  by  standard,  nondsunaging  optical  measurements. 

Surface  damage  is  discussed,  (n  Section  Il-B.  Here  we  have 
studied  the  three  types  of  mechapisms  most  likely  to  be  central  in  spr-  , 
face  damage:  hdating  py  linear  surface  absorption,  heating  by  photo¬ 
electron  absorption,  and  avblanche  ionization.  We  have  isolated 'eight 
methods  which  we  show  can  achieve  nondamaging  measurements  of  sur¬ 
face  absorptidn  and  distinguish  among  the  mechanisms.  i  \ 

•  ‘  . 

i  '  • 

B.  SE-LF- FOCUSING  BULK  DAMAGE  IN  RUBY  AND  SAPPHIRE 

1  '  ;  , 
•  i  ’  t  ! 

The  early  stages  of  this  program  were  devoted  to,  measuring 

i  '  * 

thresholds  for  bulk  damage  in  ruby  and  sap'phire.  We  obtained  samples 

*  j 

from  different  suppliers  and  carried  out  damage  threshold  measure¬ 
ments.  In  the  early  measurements  a  long  focal  length  lens  (48  cm)  was 
used  almost  exclusively  for  creating  damage.  The  results  of  the  early 
threshold  measurements  are  presented  in  Semiannual  Reports  1  and  2. 

Most  of  the  work  done  at  that  time  was  exploratory;  many  of  the  results 
had  to  be  subsequently  modified  in  the  light  of  information  obtained 
later  in  the  program.  The  early  period  was  of  great  value,  mainly  as 

'  1  !  I 
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a  time  for  gathering  clues  to  be  pursued  later  and  for  refining  the 
apparatus  and  experimental  techniques  as  well  as  gaining  experience 
and  familiarity  with  the  various  phenomena  encountered. 

In  the  following  sections  we  present  the  more  significant  results 
obtained  in  our  studies  of  bulk  damage.  Most  of  the  emphasis  is  cen¬ 
tered  on  the  phenomenon  of  self-focusing,  which  is  the  initiating 
mechanism  for  bulk  damage  in  inclusion-free  sapphire  and  ruby  and 
most  other  transparent  materials.  Many  of  the  results  presented  here 
have  contributed  to  a  deeper  understanding  of  the  self-focusing 
phenomenon. 

1.  Streak  Camera  Studies 


Some  of  the  more  significant  results  of  this  study  have  arisen 
from  a  series  of  experiments  in  which  we  observed  the  temporal  devel¬ 
opment  of  damage  tracks  using  a  fast  image  converter  camera  operat¬ 
ing  in  the  streaking  mode.  The  results  of  these  experiments  are  under¬ 
stood  in  terms  of  self-focusing  as  the  initiating  mechanism  for  bulk 
filamentary  damage;  our  experiments  constitute  the  first  direct  obser¬ 
vation  of  a  moving  self-focused  spot  in  solids.  Shortly  after  we  per¬ 
formed  our  first  streak  camera  experiments  we  came  in  contact  with 
J.  H.  Marburger  of  the  University  of  Southern  California  who  has  been 
working  with  the  theory  of  self-focusing  for  some  time.  A  treatment  of 
our  experimental  measurements  in  the  framework  of  his  theory  has 
yielded  a  satisfying  picture  of  the  dynamics  of  self-focusing  and  its 
connection  to  the  damage  tracks.  This  collaboration  has  led  to  publi¬ 
cation  inPhysical  Review  Letters  and  the  proceedings  of  the  1971  Third 
ASTM  Symposium  on  Damage  in  Laser  Materials.  The  latter  paper, 
which  is  somewhat  more  detailed  than  the  former,  is  reproduced  in  the 
following  pages,  as  published. 
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Tim  Evolution  of  Damagt  Track*  In  Sapphlr*  and  Ruby* 


Concatto  R.  Glullano 

Hughea  Raaearch  Laboratorlea 
Malibu,  California  90265 


A  faat  atraaklng  camera  haa  bean  employed  to  obaarve  the  time-evolution  of 
bulk  damage  filament*  in  aapphlre  and  ruby  on  a  nanoaacond  time  acal*.  Exter¬ 
nally  focuaed  light  from  a  mode  controlled  ruby  laaer  and  amplifier  waa  uaed  to 
Induce  damage  ir  the  aampla*  atudlad.  Either  aida-acattarad  light  at  6943  X  or 
blua-green  light  from  the  aelf-lumlnoua  damage  altea  waa  detected  by  the  atreaklng 
camera.  It  waa  found  that  damage  flrat  occur*  at  or  near  the  natural  focu*  of  the 
lena,  and  that  the  damage  track  grow*  in  an  upatream  direction  toward  the  laaer. 

The  damage  filament  la  formed  during  the  rialng  portion  of  th*  incident  laaer  pulae, 
reaching  lta  full  length  at  the  peak  of  the  pulae.  The  relatlonahlp  between  the 
locetlon  of  the  damage  at  a  given  inatant  and  th*  correapondlng  lnatantaneoua  laaer 
power  will  be  dlacuaaed  In  the  context  of  a  poaalble  aalf-focualng  mechanism. 


Key  Word*;  Bulk  damage,  filamentary  damage,  dynamic*,  moving  focua, 
sapphire,  self-focusing,  atreak  photography. 


1.  Introduction 

The  causes  of  laaer  Induced  damage  can  be  divided  Into  three  aeparate  steps.  First,  the  initiating 
step,  In  which  the  Intensity  1*  built  up  In  the  medium  by  a  mechanism  such  aa  self-focusing.  Second, 
the  stadium  responds  to  these  high  Intensities  by  undergoing  any  of  a  number  of  nonlinear  processes  such 
as  multiphoton  absorption  and/or  phonon  generation.  Third,  the  local  power  dissipation  la  high  enough 
to  exceed  the  elastic  limit*  of  the  material,  causing  catastrophic  breakdown. 

Much  Indirect  evidence  haa  existed  for  supporting  self-focusing  as  the  initiating  mechanism  for 
laser-induced  bulk  damage.  [l,2j.**  However,  different  Interpretations  for  the  phenomenon  have  been 
postulated  such  as  self-trapped  filaments  of  light,  repetitive  focusing,  and  moving  foci.  We  have 
observed  the  time-development  of  damage  tracks  In  sapphire  and  ruby  and  found  that  It  Is  most  logically 
explained  In  terms  of  a  moving  focus.  This  la  the  first  direct  observation  In  support  of  a  moving  focus 
In  solids. 

We  will  concern  ourselves  here  with  the  phenomenon  of  bulk  damage  in  sapphire,  although  much  of  the 
discussion  applies  In  general  for  most  transparent  dielectrics  in  which  filamentary  damage  Is  observed. 
The  damage  la  characterized  by  a  long  track,  aa  shown  in  Fig.  1.  Here  we  see  a  typical  damage  track 
formed  by  a  single  pulse  from  a  laser  operating  in  a  single  longitudinal  and  transverse  mode.  The 
damage  filament  la  characterized  by  a  "head"  at  the  part  of  the  track  nearest  the  sample  entrance  and  a 
tapering  "tall,"  the  end  of  which  occurs  at  or  very  near  the  focal  plane  of  the  lens  uaed  to  focus  the 
laaer  light  In  the  material.  The  track  often,  but  not  always,  displays  a  region  of  relatively  dense 
damage  at  or  near  the  tip  of  the  tall,  usually  in  the  form  of  a  crack  or  "damage  star."  as  we  examine 
the  regions  of  the  damage  track  upstream  from  the  tall  (Fig.  1  (b)),  we  begin  to  see  some  subsidiary 
fracturing  alongside  the  damage  filament.  This  fracturing  Is  randomly  located  and  becomes  Increasingly 
dense  toward  the  head  of  the  track  where  the  extent  of  subsidiary  fracturing  1*  maximum  (Fig.  1  (c>). 

If  the  Incident  light  pulse  Is  not  temporally  smooth  the  damage  track  shows  the  same  filamentary  be¬ 
havior,  but  In  addition,  there  are  occasional  regions  of  more  dense  damage  along  the  track.  An  example 
of  this  Is  shown  In  Fig.  2,  where  we  see  a  track  formed  by  a  pulse  in  which  the  laser  was  oscillating 
In  two  longitudinal  modes.  If  the  depth  of  modulation  of  the  Incident  pulse  is  large,  the  filamentary 
character  of  the  damage  track  may  not  be  present.  Instead,  one  may  see  only  a  series  of  fractured 
regions  spaced  along  the  light  path  with  gaps  of  undamaged  material  between  them.  We  also  notice  the 


*Thls  work  was  supported  by  the  Advanced  Research  Projects  Agency  under  ARPA  Order  .'34  with  Air  Force 
Cambridge  Research  Laboratories. 

‘Numerals  In  square  brackets  indicate  literature  references  appearing  at  the  end  of  this  paper. 
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general  feature  that  the  damage  filament  Iteelf  la  narrower  at  the  tall  of  the  track  (~  2  to  5  md)  and 
and  broader  at  the  head  (aeveral  tene  of  mlcrona).  We  will  later  explain  theae  qualitative  featurea  of 
bulk  damage  In  the  light  of  our  recent  obaervatlon  of  the  time-evolution  of  theae  damage  tracke. 


2.  Experimental 

The  experimental  aetup  la  ahown  In  Fig.  3.  The  oaclllator  employa  a  4  inch  long  x  1/4  Inch  diam¬ 
eter  ruby  pumped  by  two  linear  lampa  In  a  double  elliptical  pjimp  cavity.  The  ruby  cryatal  la  water- 
cooled  by  a  cloaed-cycle  refrigeration  ayatem  maintained  at  0°C.  The  high  reflectivity  mirror  la 
coated  with  a  99+X  reflectivity  high  field  damage  coating  from  Perkin  Elmer  Corporation.  Q-awltchlng 
la  accompllahed  with  a  solution  of  cryptocyanine  In  methanol  In  a  1  nxn  path  length  cell  whoae  trana- 
mlealon  la  30X  at  6943  X. 

The  temperature  controlled  (34°C)  reaonant  reflector  that  waa  dealgnad  to  optimize  longitudinal 
mode  control  conalata  of  two  quartz  etalona  and  a  quartz  apacer,  whoae  combined  effect  la  to  enhance 
cavity  modaa  aeparated  by  2  cm*1  and  to  dlacrlmlnate  agalnat  Intermediate  modea. 

Portlona  of  the  laaer  beam  are  apllt  off  In  varloua  waya  (aee  Fig.  3),  ao  that  the  power  output, 
near-and  far- field  patterna,  and  Fabry-Perot  pattern*  can  be  monitored  for  each  ahot.  Thla  la  accom¬ 
pllahed  In  the  following  way.  Light  reflecting  from  wedged  beamaplltter  W,  glvea  two  diverging  beama; 
one  of  them  hlta  the  magnealum  oxide  dlffuaer,  where  the  acattered  light  11  monitored  by  a  blplanar 
photodiode  uaed  aa  our  power  monitor.  The  aecond  beam  from  hlta  ground-glaaa  acreen  G,  where  It 
la  photographed  through  lena  L  and  the  1  m  focal  length  camera ‘focused  at  Infinity.  Thla  glvea  a  mag¬ 
nified  (~5x)  near-field  picture.  Another  portion  of  the  light  la  removed  by  beamaplltter  B,  and  hlta 
mirror  M 2,  which  can  be  placed  in  or  out  of  poaltlon  depending  on  the  uae  of  the  alignment  laaer. 

From  M_  the  light  either  goea  to  the  Fabry-Perot  Interferometer  or  can  be  partially  reflected  from 
wedged  beamaplltter  W^,  where  It  reaulta  In  a  pair  of  far-fleld  patterna.  A  0.6  neutral  denalty  filter 
la  placed  neer  the  focal  plane  of  the  camera  ao  that  the  far-fleld  pattern  and  the  Febry-Perot  pattern 
can  be  aeen  at  two  different  expoaurea.  The  two  Gian  prlama  are  uaed  ea  a  variable  attenuator  after 
the  amplifier.  Beamaplltter  B,  temple*  the  light  to  photodiode  No.  2,  which  monitor*  the  power  Inci¬ 
dent  upon  the  focualng  lena,  which  waa  dealgnad  for  minimum  apherlcal  aberration  (Special  Optica). 
Photodiode  No.  3  monitor*  the  light  after  the  aampla.  The  algnala  from  the  two  detector*  are  Integrated 
and  diaplayed  on  a  duel-beam  oacllloacopa. 

The  water  cooled  amplifier  ruby  la  6  lnche*  long  x  0.5  lnche*  diameter,  with  one  end  wedged  rel¬ 
ative  to  the  other  by  about  0.5°.  The  Input  end  of  the  amplifier  rod  la  antlref lection  coated  to  min¬ 
imize  the  chance*  of  oaclllatlon  within  the  amplifier  ltaelf.  The  ruby  rod  la  closely  coupled  to  a 
helical  flaahlamp,  which  la  pumped  with  a  power  aupply  capable  of  delivering  8  kj  In  a  3  maec  pulae. 

The  power  aupply  employa  a  pulae  ahaplng  network  of  20  aectlona,  each  aectlon  pumping  for  150  |j*ec. 

The  maximum  gain  obtained  with  the  amplifier  la  about  10  dB. 

In  all  experiment*  daacrlbed  In  thla  paper,  the  amplifier  flaahlamp  pumping  waa  held  conatant,  and 
the  amount  of  light  Incident  on  the  aample  waa  varied  by  rotating  the  flrat  of  the  two  Glen  prlama.  We 
found  that  our  amplifier  acta  a*  a  weak  negative  lena  whoae  focal  length  danands  on  optical  pumping  [3]. 
Thu*  we  fixed  the  amplifier  pumping  to  minimize  the  ahot-to-ahot  variation*  In  the  beam  characteristic*. 
Typical  output  from  the  oaclllator  plua  amplifier  la  about  150  mj  In  pulses  which  can  range  from  about 
15  to  30  nsec.  The  far-fleld  beam  profile  waa  measured  to  be  gauaalan  down  to  8%  of  the  peak  using  a 
modified  multlple-lena  camera  technique  (4]. 

The  atreak  camera  experimental  aetup  la  ahown  In  Fig.  4.  We  uaed  an  STL  Image  converter  camera 
operating  In  the  streaking  mode  to  photograph  the  aelf-lumlnou*  damage  track*  during  their  formation. 
The  camera  waa  triggered  either  by  the  sync  pula*  from  the  Pockela  cell,  when  It  waa  uaed  aa  the  Q- 
awltch,  or  from  the  laaer  light  ltaelf,  when  the  cryptocyanln*  Q-awltch  was  employed.  For  all  the  ex¬ 
periment*  described,  the  light  waa  focused  Inside  the  sample  using  a  lena  (f  ■  19  cm)  which  was  de¬ 
signed  for  minimum  spherical  aberration.  In  the  latter  experiment*  a  portion  of  the  Incident  light  waa 
split  off  and  allowed  to  enter  the  camera  directly  to  give  a  marker  atreak  which  give*  the  relation 
between  the  time  of  formation  of  a  particular  point  on  the  damage  track  and  the  peak  of  the  Incident 
pula*.  In  the  experiment*  described  here,  a  Corning  4-94  filter  was  placed  In  front  of  the  camera  lens. 
This  served  to  block  the  laaer  light  which  waa  acattered  from  tha  damage  sitae  while  paaslng  the  self- 
lumlnoua  light  In  the  blue-green  part  of  the  spectrum.  When  the  acattered  light  at  6943  A  la  allowed 
to  enter  the  camera,  essentially  the  same  behavior  1*  observed  when  It  1*  blocked.  However,  scattered 
laaer  light  from  prevloualy  formed  damage  tracks  In  the  aample  constitute*  an  Inconvenient  background 
which  Interferes  with  the  observation  of  the  track  of  Interest.  After  each  atreak  photograph  waa 
taken,  the  poaltlon  of  the  head  and  tall  of  the  damage  track  was  determined  by  examining  the  cryatal 
using  a  measuring  mlcroacop*.  These  measureaMnts  were  found  to  be  consistent  with  the  location  and 
extent  of  the  damage  Inferred  from  the  streak  photo*. 
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3.  Results  of  Streak  Measurements 


Typical  streak  phocographs  of  the  damage  In  sapphire  are  shown  In  Figs.  5  and  6.  In  Fig.  5  we  see 
the  time  development  of  the  damage  track  when  an  unmodulated  pulse  is  Incident  on  the  sample.  In  Fig.  6 
we  see  the  streak  photo  for  a  modulated  Input  pulse.  (The  750  htU  modulation  occurs  when  the  Pockets 
cell  Q-swltch  la  used.)  The  essential  features  of  these  results  to  be  noted  are  the  following. 

1.  The  track  first  appears  at  or  very  near  the  location  of  the  beam  waist.* 

2.  The  track  grows  In  the  upstream  direction  moving  toward  the  sample  entrance. 

3.  The  track  reaches  its  maximum  length  when  the  Incident  pulse  reaches  Its  peak. 

The  qualitative  behavior  of  the  time  evolution  of  these  damage  tracks  can  be  described  In  terms  of 
moving  self-focus  as  follows.  At  a  certain  critical  power  Pc ,  a  self-focus  first  appears  at  the  beam 
waist  where  the  damage  first  occurs.  As  the  power  increases  the  distance  Zf  required  to  form  a  self- 
focus  decreases,  and  the  damage  track  grows  in  the  upstream  direction.  The  minimum  corresponds  to 

the  maximum  power,  thus  the  head  of  the  track  will  occur  at  the  peak  of  the  Incident  pulse.  The  rate  of 
growth  of  the  track  depends  on  the  temporal  shape  of  the  later  pulse.  The  self-focused  spot  sweeps 
through  the  downstream  part  of  the  track  more  quickly  than  the  upstream  part,  hence  the  extent  of  damage 
Is  greater  at  the  head  than  at  the  tail. 

When  the  Incident  pulse  la  modulated,  we  expect  to  see  local  regions  of  heavy  damage  associated 
with  the  passage  of  local  maxima  on  the  leading  edge  of  the  pulse.  Each  consecutive  peak  has  slightly 
more  power  than  the  preceding  one,  and  therefore  causes  a  self-focus  which  dwells  at  slightly  smaller  z. 
A  glance  at  the  damage  track  formed  from  a  modulated  pulse  (Fig.  6)  shows  the  spacing  between  the 
heavily  damaged  regions  to  decrease  toward  the  head  of  the  track,  as  one  might  expect  from  the  above 
qualitative  explanation.  Thus,  what  might  appear  to  be  evidence  for  siultlple  or  repeated  focusing  could 
be  simply  an  artifact  of  temporal  spiking  on  the  Input  pulse,  which  causes  a  single  focus  to  pause 
occasionally  as  It  sweeps  upstream. 

Streak  photographs  of  damage  tracks  In  ruby  show  the  same  qualitative  behavior  as  that  seen  In 
sapphire,  but  the  comparison  with  the  self-focusing  theory  below  was  not  carried  out.  The  general 
features  of  damage  In  ruby  are  not  the  same  as  those  In  sapphire,  and  the  location  and  lengths  of  dam¬ 
age  tracks  are  not  as  reproducible.  One  may  attribute  this  to  the  absorption  in  ruby  at  6943  X,  but 
the  differences  have  not  been  Investigated  In  detail  (sea  ref.  [3]). 

The  reduction  of  data  taken  from  streak  and  oscilloscope  photographs  was  accomplished  by  tracing 
the  photographs  onto  graph  papar  at  a  convenient  magnification  using  an  opaque  projector  and  picking  off 
values  of  Instantaneous  power  and  damage  location  at  specific  times.  The  two  time  scales  were  synchro¬ 
nized  by  locating  the  center  (most  Intense  region)  of  the  marker  streak  and  allowing  this  to  coincide 
with  the  peak  of  the  incident  User  pulse.  The  center  of  the  marker  streak  was  determined  "by  eye"  and 
was  found  to  be  reproduclbly  locatable  to  within  about  1  nsec  (1  mm  on  the  streak  photographs).  Thus, 
we  obtain  values  of  distance  of  a  particular  point  on  the  damage  track  from  the  entrance  surface  z^  and 
the  corresponding  Instantaneous  laser  power.  An  example  of  the  traces  Is  shown  In  Fig.  7. 


4.  Comparison  of  Experiment  with  Theory 

We  will  now  attempt  to  correlate  the  above  experimental  results  with  the  self-focusing  theory  of 
Marburger  and  coworkera  (5,6).  Many  of  the  details  of  this  theory  ere  presented  In  a  following  paper 
[7],  and  we  will  only  dwell  on  the  results  as  applicable  to  our  experimental  results. 

The  abofe  referenced  theory  assume?  that  the  Induced  refractive  Index  j.n  the  medium  6n  responds 
Instantly  to  local  changes  In  the  optical  Intensity.  The  following  equation  Is  obtained  from  numerical 
solutions  of  the  nonlinear  wave  equation  (for  P  >  Pc).. 

^(P/Pc)l/2  -  0.858^2  -  0.0202  +  0.136  [ka2/zf(-)]2  (1) 

and 

*f_1(-)  -  Sj'V)  +  R'1  •  (2) 

Here,  for  an  Incident  gauaslan  equlphase  beam,  a  Is  the  e"'’  radius  of  the  Intensity  profile  at  the 
sample  entrance,  k  Is  the  wave  vector  In  the  medium,  °R  Is  the  distance  from  the  crystal  entrance  to  the 
low  Intensity  bee*  waist,  z,(R)  Is  tha  distance  from  the  sample  entrance  to  the  self-focus  when  the 
Incident  beam  has  phase  curvature  R  (R  <  0  for  a  converging  beam) ;  P  and  ?  are  the  Incident  power  and 

*In  many  of  the  photographs  the  Intensity  of  the  light  from  the  tell  of  the  track  Is  too  weak  to  record. 
However,  subsequent  examination  of  the  tracks  show  that  the  location  of  tall  Is  very  reproducible  from 
shot  to  shot. 
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th«  critical  powar  for  aelf-focuaing ,  respectively .  The  theory  predicta  that  two  self-foci  era  formed; 
they  coincide  at  the  low  lntenalty  beam  focua  at  P  P  and  at  higher  power  they  aaparate,  one  moving 
downstream  at  the  apeed  of  light  and  the  other  wing  upstream  at  a  rate  which  depende  on  the  temporal 
ahape  of  the  Input  pulse.  We  will  deal  here  with  a  comparison  of  the  theory  and  experiment  for  the 
backward  moving  self-focus. 

Equation  1  Is  a  hyperbola  whose  asymptote  la 

<P/P.)l/2  -  0.858  ♦  0.369  ka‘  ^  j  .  (3)  ’ 

1/2  -V 

Thus,  for  powers  greater  than  about  2P  at  plot  of  P  versus  if  should  give  a  atralght  line. 
Figure  8  ahows  that  this  Is  nearly  so.C  Here  we  present  data  taken  from  a  number  of  etreek  camera  meas¬ 
urements  for  different  peak  powers  and  at.  two  different  pulse  widths  t. 

The  evident  curvature  of  the  plots  In  Fig.  8  and  their  dependence  on  peak  power  end  pulae  duration 
can  be  explained  by  Including  the  slow  response  of  the  Index  change  6n  which  waa  Ignored  In  the  above 
analyale. 

Of  the  varloua  possible  mechanisms  for  8n,  thosa  arising  from  nonlinearity  of  electronic  reaponae 
[8]  and  llbratlon  (91  are  essentially  Instantaneous.  The  largest  nonlnstantaneoue  mechanism  for  in  Is 
electrostrlctlon,  which  leads  to  an  effective  nonlinear  Index  of  refraction  02,  which  decreasea  aa  the 
dimensionless  quantity  x  ■  a/ur  lncraasaa  (10].  Hera,  a  la  the  spatial  beam  radlua,  t  the  ef¬ 
fective  pulse  width,  and  u  Is  the  longitudinal  sound  velocity  In  the  medium.  Thua,  there  la  a  char* 
acterlatlc  time  a/u  for  the  electrostrlctlve  nonlinearity  to  develop,  and  hence  the  critical  power 
required  to  obtain  a  self-focua  will  depend  on  beam  alsa  and  pulaa  duration. 

The  index  change  leading  to  a  aalf-focus  at  a  particular  distance  t.  can  only  ba  induced  by  that 
portion  of  the  pulae  which  has  passed  t.  before  the  self-focus  Is  formed.  Thus,  the  effective  t 
should  ba  considered  a  decreasing  function  of  i.  when  considering  points  along  an  evolving  damage 
track  formed  on  the  leading  edge  of  the  pulse.  For  a  converging  beam,  the  beam  width  a  Increases 
with  decreasing  t.  Thus,  It  becomes  increasingly  difficult  for  electrostrlctlve  self-focusing  to  occur 
as  decreases  for  a  given  Input  pulae;  a  self-focus  formed  upstream  from  the  low  Intensity  beam 

waist  cannot  take  advantage  of  the  small  beam  radius  at  the  natural  focus. 

We  can  crudely  account  for  the  effect  of  electrostrlctlon  by  replacing  Pc'1  In  (1)  by 
FC1  +  [Pc2f  (x)]"1,  where  f  Is  an  Increasing  function  of  x,  Pcl  la  the  critical  power  due  to 
Instantaneous  mechanisms,  and  Pc2f(x)  la  the  critical  power  for  electrlctlon. 

For  our  particular  conditions  of  pulse  lengths  and  focusing,  the  fractional  change  In  a  from  the 
tall  to  head  of  the  track  Is  substantially  larger  than  the  fractional  change  In  t.  Therefore,  we 
expect  f(x)  to  be  a  decreasing  function  of  a.  casulng  positive  curvature  In  a  plot  of  P1'2  versus 
according  to  (1).  This  Is  consistent  wltn  the  data  shown  In  Fig.  8. 

Since  the  first  self-focus  occurs  at  the  low  Intensity  bean  focus,  a  Is  the  eass  for  all  curves 
at  s.  ■  | R I .  Therefore,  we  would  expect  that  x  would  decrease  for  pulses  of  longer  duration.  This 
means*that  the  eelf-focus  forms  mor*  easily  (i.e.,  at  lower  powers)  for  a  longer  pulse,  everything  else 
being  equal.  This  explains  the  vertical  separation  between  the  two  groups  of  curves  In  Fig.  8.  Sim¬ 
ilarly,  one  would  expect  that  pulsea  of  different  peak  power,  but  the  Sana  duration,  would  form  an  Ini¬ 
tial  self-focus  at  different  Incident  powers.  In  this  case  the  critical  power  Is  leas  for  lowar  peak 
power  than  for  higher  peak  power  pulses  of  the  same  duration.  Take,  for  example,  an  Incident  pulse  for 
which  the  self-focus  occurs  just  at  the  peak,  as  coa^ared  with  a  pulse  of  the  same  width  but  much 
higher  peak  power.  One  would  expect  the  self-focus  to  first  occur  at  a  higher  power  for  the  latter 
pulae,  becauae  not  enough  time  would  have  elapsed  for  a  self-focus  to  have  occurred  at  the  same  power 
as  for  the  first  pulae.  This  accounts  for  the  separation  at  the  lover  power  end  between  curves  of  the 
same  pulse  width  In  Fig.  8. 

The  slight  negative  curvature  seen  at  the  high  power  ends  of  the  curves  In  Fig.  8  Is  very  sensitive 
to  the  "synchronisation”  of  the  pulse  profile  with  the  streak  photographs.  Because  of  the  slow  response 
of  the  electrostrlctlve  mechanism,  the  maximum  upstream  excursion  of  the  damage  track  would  be  expected 
to  la>  the  pulse  peak  by  an  amount  somewhat  less  than  the  response  time.  We  did  not  attempt  to  take 
this  Into  account  In  plotting  the  data.  This  could  explain  the  negative  curvature  seen  at  the  high 
power  ends  of  several  of  the  curves  In  Fig.  8.  A  relative  shift  between  the  oscilloscope  traces  and 
the  streak  photographs  of  ths  order  of  nanoseconds  will  resove  this  curvature.  The  direction  of  the 
adjustment  Is  consistent  with  this  explanation. 

1/2  -1 

.  We  .see  from  (3)  that  the  straight  line  plot  of  P  versus  s.  shoulu  have  a  slope  of  0.369 
ka  F  1/2 .  Hence  we  can  obtain  values  for  P  from  the  slopes  of  the  curves  In  Fig.  8.  Values  for  a 
at°thS  crystal  entrance  were  obtained  by  measurement  of  the  beam  profile  at  a  number  of  points  beyond 
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the  focusing  lens.  The  range  of  critical  powers  obtained  froe  the  slopes  In  Fig.  8  Is  180  kw  to 
5.4  HU,  which  Is  consistent  at  one  extreme  with  the  predictions  of  electrostrlctlve  sslf*focuslng  for 
sapphire  (10)  end  at  the  other  with  values  of  n.  Inferred  froe  other  manifestations  of  electronic 
distortion  (ill.  1 

The  author  acknowledges  the  skilled  technical  assistance  of  C.  k.  nickel. 
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Fig.  2  Examples  of  damage  tracks  In  sapphire 

caused  by  a  temporally  modulated  pulse, 
(a)  through  (d)  progressing  from  the 
tall  to  the  head. 


9 


Fig.  3  gchematlc  representation  of  laser  an.: 

a^>llfler  and  aaaoclatad  monitoring 
apparatus. 
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Fig.  4  Strtak  canara  experimental  aatup. 
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Fig.  S  Typical  examples  of  (a)  danaga 

fllaaant,  (b)  atraak  photograph,  (c) 
oacllloacopa  traca  for  a  temporally 
smooth  Incident  pulaa. 


as  am 


Fig.  6  Typical  example  of  (a)  damage  filament, 
(b)  streak  photograph,  (c)  oacllloacopa 
trace  for  a  modulated  750  MRs) 
Incident  pulse. 
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Hg.  7  Juxtaposition  of  oacllloacopa  trace  and  streak  photograph  ahowlng  nathod  of  obtaining 
Inatontonoous  povora  and  corraapondlng  danaga  location. 


foront  Incident  pool  powers  and  pulsa  width*. 


2.  Dynamics  of  Track  Formation  in  Ruby 

The  results  obtained  with  ruby  are  qualitatively  similar  in  some 
respects  to  those  of  sapphire.  However,  there  are  a  few  differ¬ 
ences  which  emphasize  a  definite  contrast  between  the  damage  pheno¬ 
mena  in  both  materials.  A  quantitative  comparison  of  the  dynamical 
data  was  not  made  in  this  program.  Because  the  length  of  the  damage 
track  in  ruby  is  about  one-third  as  long  as  that  in  sapphire  for  the  same 
incident  power  and  focusing,  the  resolution  in  the  streak  photographs  is 
less  for  the  same  magnification.  In  addition,  the  full  length  of  the  track 
in  ruby  generally  is  not  visible  in  the  streak  photographs,  since  the 
downstream  end  does  not  glow  brightly  enough  to  be  detected.  Conse¬ 
quently,  the  visible  portion  of  the  tracks  is  very  short  and  at  the  magni¬ 
fication  used  it  is  very  difficult  to  make  reliable  measurements  from 
the  streak  photographs. 

In  addition  to  the  difference  in  track  length,  another  feature  that 
distinguishes  ruby  from  sapphire  is  that  the  track  in  ruby  is  not  a  con¬ 
tinuous  filament  even  when  the  incident  pulse  is  temporally  smooth. 

An  example  of  this  is  given  in  Fig.  1-1.  This  phenomenon  which  illus¬ 
trates  another  of  the  many  distinctions  between  ruby  and  sapphire  damage 
and  which  might  be  linked  to  the  absorption  in  ruby  remains  unexplained. 

3.  Further  Discussion  of  Theory  Compared  with  Experiment 

The  conclusion  that  the  damage  under  the  conditions  of  our 
experiments  is  the  result  of  a  moving  self-focus  rather  than  a  trapped 
filament  or  a  repetitive  focus  is  an  essential  feature  of  this  work.  A 
close  quantitative  agreement  with  the  self-focusing  theory  of  Marburger 
and  coworkers  would  predict  that  all  the  curves  in  Fig.  8  of  the  pre¬ 
ceding  manuscript  would  lie  on  top  of  one  another  and  be  essentially 
coincident  with  the  curve  of  the  manuscript  eq.  (1),  substituting  the 
proper  values  of  aQ,  R,  and  P^.  Obviously,  this  is  not  the  case,  but 
the  reasons  are  qualitatively  understood,  and  a  plausible  explanation 
for  the  differences  between  the  curves  lies  in  the  consideration  that  the 
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Fig.  1-1.  Example  of  typical  damage  track 
In  ruby  caused  by  temporally 
smooth  pulse  (f  ■  19  cm  lens), 
(a)  through  (e)  progressing  from 
tall  to  head.  Note  breaks  In 
damage  track  In  contrast  to  un¬ 
broken  track  In  sapphire. 
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nonlinearity  does  not  respond  instantaneously  to  the  optical  fields.  The 
curvature  of  the  plots,  the  shift  between  the  two  sets  of  curves  for  dif¬ 
ferent  pulse  widths,  and  the  differences  in  the  low  power  ends  of  the 
curves  with  the  same  pulse  width  can  be  qualitatively  understood  by 
invoking  a  noninstantaneous  response.  In  attempting  to  account  for  this 
empirically,  P  *  in  manuscript  eq.  (1)  may  be  replaced  by 
P  *  +  [p^ffc)]  ”  ,  where  Pcj  is  the  critical  power  for  an  instantaneous 
response,  and  P^ffx)  is  the  critical  power  for  electrostriction.  As 
mentioned  in  the  manuscript,  f(x)  is  essentially  a  function  which  makes 
it  harder  to  get  self-focusing  if  the  pulse  is  short  or  if  the  beam  size  is 
large.  This  says  then,  that  for  focused  beams  self-focusing  becomes 
more  difficult  as  the  track  progresses  back  toward  the  entrance,  at 
least  for  pulses  of  roughly  the  width  and  temporal  shape  with  which  we 
are  concerned. 

Howeve  r,  it  may  be  possible  with  our  present  data  to  determine 
f(x),  thus  obtaining  a  more  satisfactory  quantitative  agreement  than 
exists  at  this  time.  This  possibility  was  not  explored  during  this 
program. 

4.  Summary  of  Streak  Camera  Experiments 

In  summarizing  essential  results  of  the  dynamics  of  formation 
of  the  damage  filaments,  there  are  several  important  features  on  which 
the  theory  and  experiment  agree. 

As  the  incident  power  reaches  a  critical  value,  a  self-focus 
occurs  at  the  low  power  lens  focus.  This  is  the  first  point  where  a 
self-focus  is  expected  to  occur.  As  the  power  is  increased  the  dis¬ 
tance  necessary  to  obtain  a  self-focus  is  less,  and  the  self-focus  moves 
toward  the  entrance  of  the  medium.  Therefore,  the  damage  track  has 
one  end  at,  or  very  near  the  lens  focus  and  the  other  end  between  the 
focus  and  the  crystal  entrance.  The  rate  at  which  the  self-focus  moves, 
ana  consequently,  the  rate  at  which  the  damage  filament  grows  depends 
on  the  rate  of  change  of  intensity  with  time.  The  faster  the  intensity 
grows,  the  faster  the  self-focused  spot  moves  upstream.  When  the 
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incident  pulse  reaches  its  peak,  the  damage  track  reaches  its  maximum 
length.  As  the  intensity  decreases,  the  self-focused  spot  will  tend  to 
move  back  downstream.  However,  since  the  damage  has  already  been 
created,  it  interferes  with  the  propagation  back  downstream.  As  a 
result,  the  head  of  the  track  shows  more  heavy  damage.  If  the  ampli¬ 
tude  of  the  incident  light  pulse  varies  with  time,  the  self-focused  spot 
will  move  in  an  irregular  fashion,  going  rapidly  upstream  when  the 
intensity  is  rapidly  increasing  and  pausing  at  local  peaks.  The  damage 
resulting  from  this  kind  of  pulse  will  show  local  regions  of  more  heavy 
damage,  formed  during  local  maxima,  separated  by  narrow  filaments 
that  form  when  the  self-focus  is  moving  rapidly.  If  the  depth  of  modu¬ 
lation  is  large,  the  connecting  filamentary  characteristic  may  not  be 
observed,  but  only  the  regions  of  high  local  damage  along  a  line  sepa¬ 
rated  by  regions  of  no  damage.  Thus,  what  appear  to  be  evidence  for 
repetitive  or  multiple  focusing  can  be  shown  to  be  a  result  of  temporal 
structure  on  the  rising  portion  of  the  input  pulse. 

5.  Damage  Threshold  as  a  Function  of  Beam  Radius 

In  Semiannual  Report  4  we  presented  the  results  of  some  mea¬ 
surements  of  bulk  damage  threshold  in  ruby  and  sapphire  for  different 
focal  length  lenses.  The  data  are  plotted  in  Fig.  1-2  as  a  log-log  plot  of 
threshold  power  density  versus  beam  radius,  both  evaluated  at  the  beam 
waist  inside  the  sample.  A  slope  of  minus  two  in  these  plots  would  indi¬ 
cate  a  constant  power  threshold,  independent  of  lens  focal  length  that  is 
consistent  with  a  self-focusing  mechanism  in  a  material  with  an  instan¬ 
taneous  nonlinear  response  such  as  that  arising  from  electronic  distor¬ 
tion.  Instantaneous  response  implies  that  the  steady  state  conditions 
for  self-focusing  apply.  In  the  case  of  electrostrictive  self-focusing  as 
a  dominant  mechanism,  this  would  correspond  to  long  pulses  or  small 
beams.  For  pulses  of  constant  length,  the  steady  state  condition  for 
electrostriction  is  approached  as  the  beam  size  gets  smaller.  The  case 
of  large  beams  or  short  pulses  corresponds  to  the  transient  situation 
for  electrostrictive  self-focusing. 
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MINIMUM  BEAM  RADIUS  IN  SAMPLE,  /4m 


Fig.  1-2.  Threshold  power  density  versus  beam  radius 
for  ruby  and  sapphire  samples. 
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In  the  ektreme  of  the  transient  regime,  a  constant  power  density 
threshold  is  expected,  giving  a  slope  of  zero  for1  large  beam  sizes  in 

i  “ 

the  plot  in  Fig.  1-2.  If  self-focusing  by  electrostriction  were  the  sdle 
cause  of  bulk  damage,  we  would  expect  to  see  a  slope  of  minus  two  in 
the  small  beam  radius  region,  tailing  off  to  zero,  slope  at  large  radii. 

A  combination  of  contributions  frorh  both  slow  and  fast  self-focusing 
mechanisms  would,  give  a  slope  less  than  rpinus  two  whiqh  we  observe. 

Of  course,  the  above  discussion  really  applies  to  situations 

where  the  light  incident  on  the  medium  is  unfocused,  or  slightly  focused,1 

1 

because  the  self-focusing  phenomenon  depends  on  a  nonlinearity  which 
takes  place  contiguously  along  the  beam.  That  is,  the  phenomena  which 

i 

take  place  at  a  certain  plane  affect  the  self-focusing  at  points  farther 
downstream  in  the  material.  i  , 

What  will  be  expected  then  for  sharply  focused  beams  compared 
with  gently  focused  beams?  If  the  'beamiis  sharply  focused  inside  a 

1  i 

material  that  is  long  compared  with  the  length  of  the  focal  region,  the 

*  '  ■  1  i 

beam  size  is  rapidly  decreasing  with  distance  on  the  approach  to  the  low 
power  focus.  This  means  that  the  region  over  which  a  phenomenon  such 
as  electrostriction  is  effective  in  creating  selfWocusing  is  less  than  if 
the  beam  radius  were  constant  and  the  same  size  as  at  the  waist.  The 

i  * 

material  responds  to  create  a  kind  of  continuous  ac'oustic  lens  all  along 
the  beam.  But  for  a  certain  range  of  radii  the  ease  with  which  the  lens 
forms  depends  on'the  beam  radius.  Therefore,  for  a  converging  beam, 
electrostriction  might  be  expected  to  be  less  effective  as  the  external 
focusing  becomes  sharper.  In  a  sense,  then,  one  again  enters  a  region 
corresponding  to  transient  behavior.  Here  damage  threshold  would  be 

‘  i 

expected  to  occur  at  a  constant  power  density  independent  pf  l^ns  focal 

i 

length.  For  very  sharply  converging  beams,,  damage  is  expected  to 
occur  before  electr<^strictive  self-focusing  has  a  chance  to  develop. 

The  tailing  off  of  the  curves  at  small  radii  (i.e.  ,  s*harp  focusing)  in 

l  i 

Fig.' 1-2  could  be  explained  in  terms  of  the  behavior  described  above. 

The  power  densities  for  bulk  damage  for  very  short  focal  length 

I 

lasers  where  the  incident  power  is  too  low  for  self-focusing  to  take 

place  should  give  the  intrinsic  breakdown  intensity  for  the  material  of 

!  ! 

'  1  I 
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interest.  A  knowledge  of  the  intrinsic  damage  threshold  intensity 

coupled  with  self-focusing  threshold  powers  can  yield  information  as  to 

the  limiting  size  of  the  self-focused  spot.  The  plots  in  Fig.  1-2  give 

lower  limits  for  the  intrinsic  threshold  intensities  for  sapphire  and  ruby 

of  about  50  and  150  GW/cm  ,  respectively.  Since  the  curves  have  not 

really  leveled  off  at  these  values,  the  intrinsic  thresholds  must  be 

higher.  If  one  takes  typical  self-focusing  powers  of  0.  5  to  1  MW  and 

imposes  the  constraint  that  at  threshold  the  intensity  must  be  at  least 
2 

200  GW/cm  ,  values  of  from  7  to  14  pm  for  a  limiting  self-focus  spot 
radius  is  obtained,  which  is  reasonable  when  one  considers  that  the 
transverse  dimensions  of  the  damage  tracks  for  sapphire  are  about 
6  pm  in  diameter  near  the  tip  of  the  tail. 

6.  Bulk  Damage  and  Self-Focusing  for  Noncircular  Beams 

In  Semiannual  Reports  4  and  5,  some  results  of  experiments  in 
which  an  attempt  was  made  to  generate  internal  damage  in  sapphire 
using  cylindrical  lenses  for  focusing  inside  the  sample  were  discussed. 
Experiments  failed  to  obtain  damage  at  power  densities  ten  times 
higher  than  that  previously  obtained  with  circular  spots  of  the  same 
area.  This  led  us  to  suspect  that  self-focusir.g  might  be  sensitive  to 
beam  shape.  The  following  simple  ideas  were  first  employed  in  an 
attempt  to  understand  self-focusing  with  noncircular  beams.  This  non¬ 
rigorous  treatment  predicts  that  the  threshold  for  self-focusing  of 
elliptical  beams  should  be  higher  than  that  for  circular  beams,  and  the 
quantitative  result  is  not  far  from  that  predicted  using  a  more  correct 
approach  as  discussed  later. 

If  one  considers  an  elliptically  shaped  beam  of  constant  intensity 
entering  a  medium  having  an  index  nonlinearity  An  proportional  to 
power  density  and  define  the  critical  power  Pc  to  be  that  for  which  the 
angle  for  total  internal  reflection  at  the  boundary  equals  the  far-field 
diffraction  angle  characteristic  of  the  smaller  beam  dimension,  we 
find  that  the  critical  power  increases  linearly  with  the  beam  aspect 
ratio  a/b,  where  a  and  b  are  the  long  and  short  dimensions  of  the 


elliptical  beam.  The  critical  angle  for  total  internal  reflection  is 
given  by 
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If  we  expand  the  cosine  for  small  angles  and  the  right-hand  side  of 
eq.  (I- 1 )  for  small  An/n,  we  obtain 
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where  K  is  a  product  of  constants  whose  value  depends  on  the  choice  of 
units. 

The  diffraction  angle  characteristic  of  the  smaller  dimension 

b  is 


edi« 


When  the  two  angles  are  set  equal, 


P 

c 


(1-4) 


(1-5) 


results,  where  P*  is  the  critical  power  for  a  circular  beam  (a  =  b). 
The  diffraction  criterion  for  the  small  be  am  dimension  was  chosen 
because  it  assured  that  the  turning  in  of  the  rays  would  be  accomplished 
for  the  whole  beam  under  those  conditions.  It  can  be  seen  that  at  a 
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lower  power  P'  =  P*(b/a),  there  will  be  a  sufficient  nonlinearity  for 
c  c 

overcoming  diffraction  in  the  direction  of  the  long  dimension  but  the 
remainder  of  the  beam  will  continue  to  diffract. 

After  our  initial  experimental  results  were  obtained,  A.  Yariv 
of  California  Institute  of  Technology  and  J.  H.  Marburger  of  the  Uni¬ 
versity  of  Southern  California  became  interested  in  offering  more 
rigorous  theoretical  treatments  of  the  phenomenon.  As  a  result  of  this 
interest,  a  joint  paper  was  written,  presented  at  the  Seventh  Interna¬ 
tional  Quantum  Electronics  Conference  in  Montreal,  Canada,  and 
accepted  for  publication  in  Applied  Physics  Letters.  A  reprint  of  the 
manuscript  follows. 
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The  power  threshold  for  optically  Induced  hulk  damage  in  sapphire  is  a  sensitive  function  of 
the  elliptieitv  of  the  incident  lieam  shape.  Experimental  re.  nils  arc  consistent  with  a  simple 
self-focusing  theory. 


We  wish  to  report  observations  of  the  strong  influence 
of  departures  from  circularity  of  an  optical  beam  on  the 
power  threshold  for  optically  induced  bulk  damage  in 
sapphire.  This  threshold,  which  is  known  to  arise  from 
the  formation  of  a  catastrophic  self-focus  in  the  medi¬ 
um,1  is  found  to  increase  appreciably  when  the  beam 
cross  section  is  distorted  from  a  circle  to  an  ellipse, 
and  when  vertical  and  horizontal  confocal  parameters 
differ  (leading  to  ellipsoidal  phase  fronts).  These  ob¬ 
servations  are  consistent  with  the  theoretical  predic¬ 
tions  of  Vorob’yev2  and  Shvartsburg,3  which  are  extend¬ 
ed  slightly  here  to  apply  to  our  experiment. 

Our  experiments  were  performed  using  the  output  of  a 
single-mode  Q-switched  ruby  laser  and  amplifier  fo¬ 
cused  inside  sapphire  samples  with  different  lenses. 

The  laser  and  associated  monitoring  apparatus  are  de¬ 
scribed  in  detail  elsewhere.1  The  far-field  beam  pro¬ 
file  was  measured  to  be  Gaussian  down  to  8'%  of  the 
peak,  using  a  modified  multiple-lens  camera  technique. 
Typical  pulse  lengths  are  20  nsec.  The  sapphire  sam¬ 
ples  are  typically  3 -in.  -long  by  J-in.  -square  bars. 

Damage  threshold  powers  for  circular  beams  of  differ¬ 
ent  sizes  are  shown  in  Table  I  and  compared  with  re¬ 
sults  for  elliptical  beams.  The  first  three  entries  in 
Table  I  show  the  increased  threshold  power  with  beam 
size  indicating  the  influence  of  electrostriction.1  The 
dimensions  listed  in  Table  I  are  those  at  the  beam 
waist,  near  which  the  self-focus  first  forms  at  thresh- 
sld  for  circular  beams.  The  remaining  entries  in  Table 
!  show  the  effect  of  noncircular  beam  shapes.  The  fourth 
;ntry  gives  the  threshold  power  for  an  elliptical  beam 
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whose  short  dimension  is  the  same  as  the  first  entry, 
whose  long  dimension  is  the  same  as  the  third  entry, 
and  whose  area  is  the  same  as  the  second  entry.  We  see 
that  the  threshold  power  is  appreciably  higher  for  the 
elliptical  beam  compared  with  any  of  the  circular 
beams.  The  last  two  entries  in  Table  I  show  data  for 
more  elongated  elliptical  beams  where  we  were  unable 
to  reach  threshold  at  the  maximum  power  available  from 
our  system. 

Each  power  threshold  in  Table  I  was  determined  from  a 
series  of  eight  laser  shots  of  differing  peak  powers. 
When  damage  occurred,  the  damaged  region  was  found 
to  possess  a  circular  rather  than  an  elliptical  cross 
section.  In  fact,  the  appearance  of  the  damage  tracks 
formed  above  threshold  was  indistinguishable  from 
those  caused  by  circular  beams.  The  tracks  themselves 
were  confined  to  the  region  between  the  two  line  foci  of 
our  astigmatic  optical  system,  but  there  is  no  evidence 
that  they  would  not  extend  beyond  the  upstream  focus  at 
higher  powers. 

To  correlate  these  observations  with  theory,  we  have 
found  it  necessary  to  resort  to  the  “paraxial-ray— con¬ 
stant-shape”  analysis  of  self-focusing  outlined  in  Ref.  5. 
More  accurate  numerical  solutions  of  the  nonlinear 
wave  equation,  such  as  those  employed  in  Ref.  1,  are 
extremely  difficult  to  obtain  in  this  case,  because  the 
noncircular  beam  shape  requires  an  additional  degree 
of  freedom  in  the  computer  code.  Nevertheless,  the  ap¬ 
proximate  analysis  has  had  some  qualitative  success,' 
and  leads  to  simple  expressions . 
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TABLE  I.  Comparison  ol  damage  thr cuboid  powers  for  circular  ami  elliptical  beams  of  different  sizes. 


1  ype  of  beam 

Radius*  (pm) 

Area  (pms) 

Damage  threshold  power 
(MW) 

Astigmatic  focal  length* 
in  medium  (eni) 

Circular  beams 

14 

6.  3  x  io5 

0.51  *  (1.04 

•  •  • 

(dimensions  at  low- 

:i7 

44  x  10s 

0.  74  1  0.  07 

•  •  • 

intensity  waist) 

ton 

30H  *los 

1.23*  0.  10 

•  •  • 

Elliptical  Ism  ms  (dimen- 

14  xioo 

44  xio*’ 

O.Oi  0.50 

J .  D!) 

sions  at  upstream  low- 

14  *1200 

535  x10s 

>  17e 

«C 

intensity  waist) 

7*1200 

270  x10s 

>17' 

or 

*Thc  beam  radius  is  defined  as  the  1  /c  radius  for  the  intensity.  Iondon,  19(55),  p.  1 7 1 1 . 

•This  Is  the  distance  between  the  two  line  foci  {see  M.  Born  'Unable  to  reach  threshold  at  this  power, 

and  E.  vVolf,  Principles  of  Optics,  :ird  Ed.  (Pergamon, 


Assuming  an  optical  field  of  the  form 


=  0  occurs  only  if  P'PX  exceeds 


£  =  j£0exp(tV)exp[f(*e-w/)|4c.c.,  ncc  =nc+  -  jf'j  ■  (7) 

E.  -  £L(z)exp[-  i(.v»  fa>  +  v2 '&)], 

For  P=  l\ ,  the  self-focus  occurs  at  z  -  zn,  where 


we  find  by  the  method  of  Ref.  5  the  following  equations 
for  the  horizontal  and  vertical  beam  parameters  a  and 

/>: 


'1  -  i  _  JL 

riz*  ~  a*  <Ph  ’ 

,,2{W>  -1  JL 

( U. *  ~  IP  ~  alP  • 


(1) 

(2) 


Here  ii~P/P, |,  where  P  is  the  total  power  and  Px  =n*c/ 

4 ijk2  in  Gaussian  units.  The  first  terms  on  the  right- 
hand  sides  of  Eqs.  (1)  and  (2)  represent  the  effect  of 
diffraction,  while  the  second  terms  arise  from  the  non¬ 
linear  contribution  s.f2EJ  to  the  dietectric  constant  n  is 
the  linear  refractive  index.  Vorob’yev,2  who  derived 
Eqs.  (1)  and  (2)  in  a  different  way,  has  shown  that  they 
imply 

^(n»+ft2)  =  0,  (3) 

which  can  be  integrated  to  give 

4  b2)  =  |*v0 4  b%)  +  (2/oA)(ne  -  q)|z2 

4  2^(«X  +  W*+*Vo +'»?)•  (4) 

Here  n„  and  b0  are  the  axes  of  the  ellipse  formed  by  the 
<’■'  points  of  the  transverse  intensity  profile  at  z  =0, 
and  i)e  =  nn/2bn  4  l>„/2nn.  The  initial  rates  of  change  n„ 
and  b„  are  simply  related  to  the  principal  radii  of  curva¬ 
ture  ft,  and  /i„  of  the  phase  fronts  at  z  =0: 

R«  "o  ‘ "n »  R*  “  ^O'  ^o’ 

Using  Eq  (4),  the  reader  may  easily  find  that  the  on- 
axis  intensity,  which  is  inversely  proportional  to  2 ab 
+/J2,  becomes  infinite  at  the  point 

*/.=*«A(>»/H.  -If"*.  (5) 

for  plane  incident  phase  fronts  and  rj  >  7^ .  It  is  therefore 
clear  that  the  critical  power  for  self-focusing  is 

(ft.  =  «,=-»),  (6) 

which  equals  Px  for  circular  beams.  For  finite  values  of 
ft,  and  ft,,  a  real  solution  for  z  of  the  condition  a2  4  ft2 


-  +  *8)  =«S  'R.  +  AS  (8) 

For  higher  powers  the  self-focus  will  be  found  at 

*' =  i+(«„  'Ll1'2  '  (9) 

For  the  fourth  entry  in  Table  I,  the  beam  parameters 
are  //„  =  0. 179  mm,  />„=0.232  mm,  ft,  =  4. 05  cm,  and 
ft,  =  7.04  cm.  The  measured  power  threshold  for  the 
onset  of  bulk  damage  was  from  five  to  twelve  times 
higher  than  that  for  circular  beams  of  comparable  di¬ 
mensions,  and  the  damage  first  appeared  at  z=  5. 8  cm 
from  the  entrance  face  Assuming  that  the  damage 
should  first  appear  when  a  catastrophic  self-focus  forms 
at  P=1cePx,  we  expect  from  Eqs.  (7)  and  (8)  that  the 
threshold  power  should  be  increased  by  ijCf  —  12.7  and 
that  the  damage  at  threshold  should  appear  at  z0  =  5.51 
cm  The  latter  figure  agrees  well  with  the  measured 
value,  and  the  change  in  critical  power  is  certainly  of 
the  right  order  Unfortunately,  our  knowledge  of  the 
influence  of  the  transient  electrostrictive  contribution 
to  tj  is  too  poor  to  allow  us  to  know  which  of  the  first 
three  entries  in  Table  I  we  should  choose  as  Px.  The 
first  entry  (Px  =0  51  MW)  gives  excellent  agreement 
with  theory,  but  the  electrostrictive  response  for  an 
elliptical  beam  can  hardly  be  as  large  as  that  for  a 
circular  beam  whose  diameter  equals  the  small  radius 
of  the  ellipse.  Nevertheless,  the  agreement  with  theory 
is  sufficiently  good  to  allow  us  to  conclude  that  the 
damage  morphology  and  threshold  are  well  understood 
in  terms  of  the  self-focusing  theory 

The  theory  outlined  above  gives  simple  formulas  for  the 
powers  and  distances  at  which  the  quantity  a 2  +  IP  van¬ 
ishes,  but  the  intensity  can  become  infinite  only  if  the 
area  nab  of  the  beam  vanishes.  A  closer  examination  of 
the  solutions  of  Eqs.  (1)  and  (2)  shows  that  a  and  b  both 
vanish  at  the  same  point.  Prior  to  this  point,  the  ratio 
a  b  approaches  unity.  This  is  consistent  with  the  ob¬ 
servation  that  the  cross  sections  of  the  damaged  region 
have  a  circular  shape  along  the  entire  length  of  the 
damage  track.  We  remind  the  reader  that  at  some  point 
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between  the  astigmatic  tine  foci  the  low-intensity-beam 
cross  section  must  be  circular,  but  elsewhere  it  is 
elliptical 

These  results  imply  that  difficulties  which  arise  from 
bulk  damage  due  to  self-focusing  in  high-power  optical 
systems  can  be  avoided  or  ameliorated  by  the  proper 
choice  of  beam  parameters . 

The  authors  are  grateful  to  V.  Evtuhov  for  helpful  dis¬ 
cussions  and  for  providing  English  translations  of  Refs . 
2  and  3. 
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7. 


Further  Comments  on  Elliptical  Beams 


The  results  of  the  self-focusing  theory  for  gaussian  elliptic 
beams  compared  with  circular  beams  are  summarized  in  Table  1-1. 
Table  1-2  is  a  comparison  of  experimental  results  with  theory  for  both 
critical  power  and  self-focusing  distance.  As  mentioned  in  the  preced¬ 
ing  manuscript,  the  comparison  between  theory  and  experiment  is  of 
the  right  order  for  the  critical  power  and  is  very  good  for  the  self- 
focusing  distance.  In  view  of  the  fact  that  the  theory  is  approximate 
and  that  the  effects  of  noninstantaneous  response  are  ignored  com¬ 
pletely,  we  are  confident  that  the  treatment  is  satisfactory. 

Certainly  a  great  deal  of  work  must  be  carried  out  in  this  area 
before  a  complete  appreciation  of  the  pheonomenon  can  be  realized. 
However,  we  feel  that  these  results  point  to  a  very  promising  solution 
of  bulk  damage  problems  which  arise  from  self-focusing.  Their 
applicability  for  providing  solutions  to  this  problem  in  a  number  of 
high  power  laser  and  amplifier  systems  is  evident  and  should  be 
pursued. 

8.  Use  of  Elliptical  Beams  for  Determining  Intrinsic  Damage 

Threshold  Power  Densities 

As  mentioned  in  Section  I-B-5,  the  measurement  of  damage 
threshold  for  short  focal  length  lenses  can  yield  information  on  intrinsic 
breakdown  power  densities  if  such  damage  can  be  initiated  in  the 
absence  of  self-focusing.  However,  for  this  type  of  experiment,  when 
threshold  is  reached  the  damage  always  occurs  at  the  low  power  lens 
focus  whether  or  not  self-focusing  occurs.  Hence  it  is  difficult  to 
know  whether  in  such  experiments  self-focusing  can  be  avoided 
completely. 

With  focused  elliptical  beams,  however,  one  can  more  easily 
distinguish  between  the  two  possibilities.  The  location  of  the  first 
self-focus  for  externally  focused  elliptical  beams  occurs  at  a  point 
between  the  two  line  foci  as  given  by  eq.  (8)  in  the  preceding 
manuscript  or  in  Table  1-1.  This  point,  however,  does  notin  general 
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TABLE  1-1 


Theoretical  Comparisons  of  Pc  and  Z£ 
for  Circular  and  Elliptical  Beams 


Type  of  B«tm 

Critical  Power  for 

Self- Focusing,  P 

Self- Focusing 

Di  stance, 

a.  at  P  =  P 

f  c 

Collimated 

Circular 

Pc  =  Pj-  nSc/4E2k2 

to 

Elliptical 

P  -  IP, 

'  o  o' 

Ext.ra.lly 

Focused 

Circular 

p  =  p. 
c  1 

VR 

Elliptical 

Pc  •  VP1 

<taobo»2/l  1  \2 

Tcc  \  +  ~ TT"  (tT/M 

-«»o2  +  boZ» 

*1  ‘  1 - 1 — 

•J'VVS 

•  ,  b  «  bum  r.dli  at  entrance  to  medium 

o  o 

Rfc,  «  radii  of  curvature  of  phase  fronts  at  entrance  to  medium. 

T628 


TABLE  1-2 

Comparison  of  Theory  with  Experiment 


Ratio  of  Critical  Power  of 
Elliptical  Beam  to  Circular 

Beam,  Pc/P] 

Self-Focusing  Distance 
at  Threshold, 

Calculated 

12.  7 

5.  51  cm 

Observed 

§  OdD 

12.0  8.1  5.0 

5.  8  cm 

Experimental  beam  parameter*  at  entrance  to  medium: 

a  *  0.179  mm  R  =  -4.05  cm 

o  a 

b  =  0. 232  mm  R,  =  -7. 04  cm 

o  b 

T629 
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correspond  to  the  point  of  highest  intensity  for  linear  propagation.  The 
experimental  conditions  can  be  chosen  so  that  the  location  of  the  first 
self-focus  and  the  location  of  the  maximum  intensity  point  for  a  linearly 
propagating  beam  can  be  made  to  differ  arbitrarily.  Hence  one  can  tell 
whether  the  damage  near  threshold  arises  from  self-focusing  or  by 
intrinsic  breakdown  depending  on  where  it  is  located  in  the  medium.  If 
damage  occurs  at  the  high  intensity  point  predicted  on  the  basis  of 
linear  propagation,  then  the  intrinsic  breakdown  intensity  can  be 
determined  from  a  knowledge  of  the  power  and  spot  size. 

9.  Optical  Pumping  Experiments 

Early  in  the  program  we  presented  the  results  of  a  number  of 
experiments  where  we  optically  pumped  the  sample  while  performing 
damage  threshold  measurements.  (Details  of  this  experimental  setup 
and  procedure  were  given  in  Semiannual  Report  2. )  We  observed  a 
large  amount  of  scatter  in  the  data  because  of  interference  from 
plasma  formation  at  the  entrance  surface  and  were  not  able  to  state 
clearly  whether  or  not  the  threshold  is  affected  by  optical  pumping^-  At 
a  later  time,  we  repeated  the  pumping  experiments  with  better  preci¬ 
sion  and  found  a  slight  decrease  in  threshold  with  optical  pumping. 

This  is  shown  in  Fig.  1-3.  The  upper  graph  shows  the  later  more 
precise  data  for  Cz  ruby  LI 07.  The  lower  three  graphs  are  earlier 
data  presented  for  comparison.  We  see  about  20%  decrease  in  relative 
threshold  for  the  most  extreme  pumping  case  compared  with  the 
unpumped  case. 

We  also  discussed  in  Semiannual  Report  2  an  interesting  effect 
concerning  the  location  of  damage  in  the  s  imple  as  a  function  of  optical 
pumping  of  ruby.  We  noted  that  the  beginning  of  the  damage  track 
(i.  e. ,  the  upstream  end)  is  shifted  downstream  with  optical  pumping. 

(We  know  now  from  the  streak  photographs  that  the  upstream  end  of  the 
damage  track  is  really  the  last  part  formed.  )  From  this  shift  in  the 
damage  track  in  ruby,  we  inferred  that  the  optical  focus  was  shifting, 
but  we  now  know  from  our  beam  profile  measurements  that  the  focus 
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Fig.  1-3. 

Relative  Bulk  Damage  Threshold  as  a 
Function  of  Optical  Pumping  for  Dif¬ 
ferent  Samples.  The  Thresholds  are 
Normalized  to  Unity  for  the  Unpumped 
Sample.  Dashed  Line  —  No  Damage; 
Solid  Line  -  Damage.  The  Three  Lower 
Graphs  are  Reproduced  from  the  Last 
Report;  the  Upper  Graph  shows  Results 
of  more  Reproducible  measurements 
taken  during  this  Reporting  Period. 
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is  located  downstream  from  where  we  thought  it  was.  In  fact,  from 
our  streak  camera  measurements,  we  know  that  the  damage  track  lies 
between  the  natural  focus  and  the  entrance  surface.  In  sapphire,  the 
track  begins  at  or  very  close  to  the  natural  focus  and  progresses 
upstream;  in  ruby  the  points  of  beginning  and  end  of  the  track  occur 
substantially  upstream  from  the  natural  (i.  e. ,  low  power)  focus. 

Optical  pumping  in  ruby  shifts  the  track  toward  downstream  locations 
as  well  as  making  it  longer.  That  is,  it  approaches  the  behavior  for 
sapphire. 

Precise  data  on  this  point  are  not  available  at  this  time. 
However,  we  wish  to  point  out  that  the  shift  in  track  location  is  not 
necessarily  related  to  a  defocusing  due  to  pumping  solely,  but  may  be 
connected  with  the  fact  that  pumped  ruby  behaves  more  like  sapphire 
simply  by  virtue  of  its  being  pumped. 

C.  SURFACE  DAMAGE  STUDIES 

1.  Relation  Between  Surface  Damage  and  Plasma  Formation 

During  the  latter  one-third  of  this  program  more  emphasis  has 
been  placed  on  the  problem  of  surface  damage  rather  than  bulk  damage. 
Our  efforts  in  this  area  have  been  to  measure  the  thresholds  for  ruby 
and  sapphire  under  different  conditions  ard  to  attempt  to  understand 
the  connection,  if  any,  between  the  damage  and  the  surface  plasmas 
which  occur.  A  paper  has  beer,  written  describing  this  work;  it  has 
been  accepted  for  publication  in  the  IEEE  Journal  of  Quantum 
Electronics.  An  expanded  version  of  this  work  was  presented  at  the 
Fourth  ASTM-NBS  Symposium  on  Damage  in  Laser  Materials, 

Boulde  Colorado,  June  1972  and  will  be  published  in  the  Proceedings. 
The  manuscript  is  presented  in  the  following  pages. 
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The  Relation  between  Surface  Danage  and 
Surface  Plasma  Formation* 


Conce t to  R .  Gi ul i ano 

Hughes  Research  Laboratories 
Malibu,  California  90265 


A  number  of  experiments  are  described  in  which  entrance  and 
exit  surface  1 aser- ( nduced  damage  on  sapphire  samples  are  studied. 
Damage  thresholds  are  measured  as  a  function  of  beam  size  and 
divergence  at  the  surfaces.  Temporal  evolution  of  surface  plasmas 
is  studied  using  a  streaking  camera,  and  spatial  di f ferent iu t i on 
of  surface  plasmas  is  accomplished  by  allowtnn  the  liqht  to  strike 
the  surfaces  away  from  normal  incidence.  It  is  found  that  the 
entrance  surface  has  two  plasma  components,  an  air  plasma  which  is 
directed  along  the  light  beam  and  an  "explosion”  plasma  which  is 
directed  normal  to  the  surface.  The  exit  surface  has  only  the 
explosion  plasma.  It  is  concluded  that  surface  plasmas  are  a 
result  of,  rather  than  a  cause  of  surface  damage. 


Key  Words:  Entrance  and  exist  surface  damage,  sapphire, 
spatial  and  temporal  plasma  resolution,  streak  camera 
experiments,  surface  plasmas. 


1  .  Introduction 

It  has  been  known  for  some  time  that  laser-induced  surface  damage  has  distinctly 
different  physical  characteristics,  depending  on  whether  the  entrance  or  exit  surface 
is  damaged.  [1-4]'  To  a  certain  extent  the  conditions  under  which  the  damage  is  gen¬ 
erated  (e.g.,  focusing,  surface  finish)  are  important  in  determining  some  of  the  details 
of  the  appearance  of  the  damage,  but  as  a  general,  rule  entrance  surface  damage  is  less 
extensive  and  more  difficult  to  detect  visually  than  exit  surface  damage,  especially 
when  generated  under  conditions  of  gentle  focusing.  Under  such  conditions  for  sapphire 
the  distinguishing  features  are  that  exit  surface  damage  can  be  most  easily  described 
as  relatively  deep  fractured  areas,  while  entrance  surface  damage  is  relatively  shallow 
crazed  areas.  Scanning  electron  micrographs  of  both  types  of  damage  are  presented  in 
Figures  1  and  2. 

The  occurrence  of  1 aser- 1 nduced  surface  damaqe  is  known  to  be  accompanied  by  a 
bright  flash  of  light  which  comes  from  a  spark  or  plasma  at  the  surface.  [1-5]  (See 
Figure  3)  In  fact,  the  occurrence  of  a  visible  spark  has  been  held  to  be  a  criterion 
for  surface  damage.  [4]  Spectra  taken  of  exit  surface  plasmas  show  that  the  light 

contains  spectral  lines  characteristic  of  the  material  being  Irradiated.  [1]  The  con¬ 
nection  between  the  surface  plasma  and  the  damage  was  not  established  in  the  early  work. 
Whether  or  not  the  plasma  is  a  result  of  or  a  cause  of  the  surface  damage  or  whether  it 
is  an  independently  occurring  phenomenon  that  simply  accompanies  the  damage  phenomenon 
was  not  determined. 


Work  supported  by  the  Advanced  Research  Projects  Agency  under  ARPA  Order  No. 
Air  Force  Cambridge  Research  Laboratories. 


14  34  with 


1 


FI gures 


in  brackets 


Indicate  the  literature  references  at  the  end  of  this  paper. 


29 


Fig.  1.  Scanning  Electron  Nlcrograohs  of  Exit  Damage  Sites  for 
Sapphire.  Pulse  Energy  Increases  from  a  through  d; 
e  and  f  are  Magnified  Views  of  Site  Shown  In  d. 


) 
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Fig.  2. 

Scanning  Electron  Micrographs 
of  Entrance  Damage  for  Sapphire. 
Pulse  Energy  Increases  fr6m  a-c. 


I 

Example  of  Entrance  and  Exit  Surface 
Plasmas  Formed  When  Light  From  a  Q- 
Switched  Ruby  Laser  Passes  Through  a 
jSample.  .Sample  Length  -  3  In.  Light 
Is  Traveling  from  Left  to  Right) 
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Fersman  and  khazov  [3]  have  suggested  that  the  surface  damage  is  caused  by  a 
pi  asr.a-surface  Interaction  and  have  explained  the  difference  between  entrance  and 
exit  surface  damage  using  a  model  in  which  both  the  plasmas  grow  in  an  upstream  direc¬ 
tion  as  they  develop  in  time,  an  extension  of  the  observations  for  air  breakdown  plasmas. 
[6]  The  exit  damage  would  then  be  expected  to  be  more  extensive  than  entrance  damage 
because  the  exit  plasma  grows  into  the  surface,  while  the  entrance  plasma  grows  away 
from  the  surface.  Fersman  and  Kliazov  observe  a  compress  1  ona  1  wave  which  propagates 
inwardly  from  the  exit  surface.  Compress ional  waves  arising  from  damage  sites  in  solids 
have  been  observed  by  others  using  both  schlieren  and  holographic  technigues.  [7.P] 
While  the  above  picture  can  predict  a  difference  between  the  entrance  and  exit  damage, 
it  was  not  experimentally  confirmed. 


Another  possible  explanation  for  the  difference  between  exit  and  entrance  damage 
is  that  exit  damage  might  be  somehow  inf'uenced  by  sel f- focus i ng  of  the  beam  in  the 
material,  while  the  beam  at  the  entrance  would  not  be  so  affected. 

We  have  carried  out  a  number  of  experiments  in  an  attempt  to  determine  more 
explicity  the  connection  between  the  surface  damage  and  the  plasmas  and  to  clarify  the 
nature  of  the  damage  mechanism  or  mechanisms. 

All  the  experiments  described  in  this  paper  were  carried  out  using  a  single  longi¬ 
tudinal  and  transverse  (TEM00)  mode  Q-switched  ruby  laser.  The  laser  and  associated 
monitoring  apparatus  are  described  in  detail  elsewhere.  [Z]  Typical  laser  pulse  width 
and  peak  power  are  20  nsec  ( FWHM)  and  1  MU,  respectively. 

2.  Damage  Threshold  as  a  Function  of  Beam  Size  and  Diveroence 

One  of  the  questions  which  arises  in  the  consideration  of  exit  surface  damage  is 
concerned  with  the  possibility  that  the  damaqe  at  the  exit  surface  might  be  influenced 
by  what  the  beam  does  in  the  bulk  of  the  material.  For  example,  if  some  se  1  f- focus i ng 
or  other  nonlinear  phenomenon  takes  place,  it  might  affect  the  beam  so  that  the  inten¬ 
sity  at  the  exit  face  will  be  different  from  that  expected  on  the  basis  of  linear  prop¬ 
agation  theory.  Also,  some  phenomenon  which  might  cause  the  generation  of  acoustic 
waves  directed  toward  the  exit  surface  would  be  affected  by  how  the  beam  propagates 
inside  the  sample  before  reaching  the  exit  surface.  A  series  of  experiments  on  sapphire 
was  carried  out  in  which  the  sample  was  moved  relative  to  the  focusinq  lens,  so  that 
the  spot  size  at  the  surface  would  vary  over  a  wide  range  and  also  so  that  the  beam 
could  be  either  converging  or  diverging  at  the  surface,  depending  on  whether  the  beam 
waist  lay  downstream  or  upstream  from  the  surface. 

The  result  of  these  experiments  is  shown  in  Figure  A.  In  these  experiments  a  19  cm 
lens  for  focusing  the  light  *rom  the  mode- con  trol  1  ed  ruby  laser  was  used.  The  sapphire 
samples  were  moved  in  5  mm  increments  relative  to  the  lens,  and  surface  damage  threshold 
was  measured  for  each  location.  The  data  on  the  left  half  of  the  figure  correspond  to 
the  case  of  rays  converging  at  the  exit  surface  and  the  data  on  the  right  half  to  diverg¬ 
ing  rays.  The  beam  area  at  the  surface  is  also  plotted  as  a  function  of  distance  in  the 

figure. 

Figure  4  shows  that  the  results  for  exit  surface  damage  threshold  power  In  sapphire 
sample  L134  fall  very  close  to  the  beam  area  curve.  The  fact  that  they  fall  so  closely 
is  a  coincidence  of  the  scale  and  the  units  chosen,  but  the  fact  that  they  have  the  same 
functional  dependence  on  distance  is  significant  in  that  it  strongly  indicates  that  the 
damage  threshold  is  a  power  density  phenomenon.  Also,  the  fact  that  the  threshold  is 
independent  of  whether  the  beam  is  converging  or  diverging  at  the  exit  surface  indicates 
that  the  damage  seems  to  be  strictly  a  surface  or  near-surface'  phenomenon.  The  exit 
surface  damage  threshold  for  this  sample  is  about  1  GW/cm?.  The  data. for  entrance  sur¬ 
face  damage  for  sapphire  sample  C129  are  shown  in  Figure  5.  A  similar  dependence  on  lens 
to-sample  distance  Is  seen  in  the  plot  but  in  this  case  the  scatter  in  the  points  is  some 

what  greater.  Here  the  beam  area  curve  is  plotted  on  a  different  scale  than  in  Fig.  4 

to  show  that  the  dependence  of  threshold  power  with  distance  is  essentially  the  samp  as 
the  dependence  of  beam  area  with  distance,  again  indicating  that  damage  threshold  is  a 
power  density  phenomenon.  The  entrance  damage  threshold  for  this  case  is  about  1.5 
GW/cm? . 

3.  Surface  Damage  without  Plasma  Formation 

He  have  carried  out  a  few  experiments  in  which  surface  damage  was  generated  without 
observing  an  accompanying  plasma.  These  experiments  were  conducted  on  sapphire  samples 
and  the  observations  apply  only  to  the  exit  surface  (entrance  damage  was  not  pursued  in 
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Fig.  4.  Exit  Surface  Damage  Threshold 
Power  for  Sapphire  Sample  as 
a  Function  of  Surface  Distance 
fro*  Beam  Heist  for  19  cm  Lens. 
Beeai  Area  Is  also  Plotted. 


Fig.  S.  Bean  Area  and  Entrance  Surface 
Damage  Threshold  as  a  Function 
of  Surface  Distance  froai  Bean 
Halst  for  19  cm  Lens. 


this  series  of  experiments).  In  these  experiments  we  observed  the  pxlt  surface  during 
laser  Irradiation  both  visually  and  photographically  under  different  conditions  of 
optical  filtering  and  exposure.  It  was  found  that  when  the  sample  was  exposed  to 
focused  radiation  very  close  to  the  threshold,  damage  could  be  generated  In  the  form 
of  a  series  of  micropits  without  the  simultaneous  accompaniment  of  a  detectable  plasma 
or  spark . 

A  series  of  experiments  was  initiated  In  which  the  sample  was  Irradiated  well 
above  damage  threshold,  and  the  plasma  observed  as  the  power  was  gradually  reduced, 
each  time  exposing  a  fresh  area  of  the  surface.  During  this  series,  the  accompanying 
plasma  was  seen  to  decrease  In  size  and  brightness  as  the  laser  power  was  reduced, 
until  a  power  point  was  reached  at  which  the  spark  was  no  longer  visible.  On  a  number 
of  occasions  after  laser  Irradiation  at  a  level  Just  below  that  for  generating  a 
visible  plasma,  it  was  possible  to  detect  damage  in  the  form  of  a  number  of  inicropits 
when  the  sample  was  carefully  examined  through  a  microscope  with  properly  chosen  illum¬ 
ination.  This  ml  cropl  ttlng,  which  (s  seen  only  very  close  to  threshold,  is  less  exten¬ 
sive  than  the  single  fracture  crater  that  occurs  mo'e  often  at  slightly  higher  powers 
when  the  plasma  Is  clearly  visible.  It  should  be  noted  that  at  the  times  when  no  plasma 
Is  visible,  It  is  possible  to  see  small  amounts  of  laser  light  scattered  from  the  tiny 
surface  pits  which  are  being  formed. 

The  results  of  these  observations  suggest  that  It  is  possible  to  obtain  exit  sur¬ 
face  damage  without  plasma  formation.  Of  course,  the  nonobservability  of  a  plasma  does 
not  necessarily  Imply  that  It  Is  not  present,  and  these  results  alone  do  not  allow  us 
to  form  any  definite  conclusions  regarding  the  connection  between  the  damage  and  the 
plasma . 


4.  Time  Evolution  of  Surface  Plasmas 

One  reason  for  the  interest  In  studying  the  temporal  evolution  of  the  surface 
plasmas  has  been  to  determine  the  direction  of  growth  of  the  plasmas  relative  to  the 
crystal  surface,  which  would  lead  to  a  possible  explanation  of  the  differences  between 
exit  and  entrance  damage  and  to  test  the  mechanism  proposed  by  Fersman  and  Khazov. 

He  have  employed  an  STL  Image  converter  camera  operating  In  the  streaking  mode  for 
observing  the  growth  of  the  entrance  and  exit  surface  plasmas  (see  Floure  6).  For  these 
experiments,  the  output  from  our  single-mode  ruby  laser  was  focused  with  a  19  cm  lens 
on  the  desired  surface  of  the  sapphire  sample  being  examined.  Typical  laser  output  Is 
about  20  mj  In  a  20  nsec  pulse,  which  gives  a  power  density  at  the  focus  of  approxi¬ 
mately  10  GW/cm*,  about  10  times  above  damage  threshold.  For  the  streak  camera  experi¬ 
ments,  a  portion  of  the  Incident  laser  light  was  allowed  to  enter  the  camera  directly 
to  serve  as  a  marker  streak  for  synch roni zatlon  purposes. 

Streak  photographs  showing  the  temporal  development  of  the  surface  plasmas  are 
shown  In  Figures  7  and  8.  Because  the  photographic  setup  is  such  that  the  camera  views 
the  surface  at  an  angle  less  than  90°  to  the  normal  (see  Figure  6),  both  the  plasma  and 
Its  reflection  In  the  surface  are  often  seen.  The  fact  that  there  are  Indeed  reflec¬ 
tions  was  confirmed  by  placing  an  appropriately  oriented  Polaroid  filter  in  front  of  the 
camera.  In  which  case  the  reflection  from  the  dielectric  surface  disappears  as  expected. 

Essential  features  of  these  streak  photographs  to  be  noted  are 

•  The  plasmas  grow  away  from  the  surfaces  In  both  cases,  the  entrance 
plasma  growing  1 n  the  upstream  direction  and  the  exit  plasma  In  the 
downstream  direction. 

t  The  entrance  plasma  has  two  temporal  components,  whereas  the  exit 
plasma  has  one  temporal  component. 

0  The  extent  of  growth  and  luminosity  of  the  faster  growing  component 
of  the  entrance  plasma  depends  on  the  Incident  laser  pulse  In  that  It 
ceases  to  expand  and  to  emit  light  when  the  laser  terminates. 

0  The  slower  growing,  longer  lived  component  of  the  entrance  plasma 
expands  at  about  the  same  rate  as  the  exit  plasma.  These  plasmas 
continue  to  expand  outward  and  emit  light  long  after  the  laser  has 
terminated. 

0  Typical  velocities  of  expansion  of  the  luminous  plasma  fronts  are 
•v2  x  10®  cm/sec  for  the  short-lived  entrance  component  and  t-1  x  10® 
cm/sec  for  the  long-lived  entrance  component  and  for  the  exit  plasma. 
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Fig.  6.  Schematic  Diagram  Showing  Setup  used  for  Photo¬ 
graphing  Surface  Platmas  both  for  Still  and 
Streak  Photography. 
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Fig.  7. 

Streak  Photograph  of  Entrance  Plasma 
on  Sapphire  Surface  Showing  Two  Tem¬ 
poral  Components. 

Fig.  8. 

Streak  Photograph  of  Exit  Plasma  on 
Sapphire  Surface. 

Several  ideas  are  suggested  by  these  observations.  First,  the  fact  that  both 
entrance  and  exit  plasmas  grow  away  from  the  surfaces  suggests  that  the  mechanism 
proposed  by  Fersman  and  Khazov  is  Incorrect.  From  the  time  dependences  of  the  plasmas 
relative  to  the  incident  laser  pulse  it  can  be  Inferred  that  the  short-lived  entrance 
component  is  a  plasma  supported  in  the  ambient  atmosphere  and  sustained  directly  by 
the  laser  beam  Itself,  while  the  longer-lived  components  at  both  surfaces  result  from 
a  kind  of  micro-explosion  caused  by  the  deposition  and  subsequent  release  of  energy  at 
or  near  the  surface. 


5.  Spatial  Differentiation  of  Surface  Plasmas 

To  test  these  hypotheses,  further,  we  devised  additional  experiments  that  were 
designed  to  allow  a  spatial  differentiation  of  the  plasmas.  The  reasoning  behind  these 
experiments  is  as  follows.  A  plasma  sustained  directly  by  the  light  beam  would  be 
expected  to  exist  only  where  the  light  is  present  and  to  be  directed  along  the  light 
beam  Independent  of  the  angle  at  which  the  light  enters  the  material.  A  plasma  which 
results  from  a  surface  explosion  caused  by  a  deposition  of  energy  at  the  surface  would 
be  expected  to  extend  in  a  direction  no  rmal  to  the  surface  no  matter  at  what  angle  the 
light  should  enter  the  material. 

Figure  9  shows  time-integrated  photographs  of  the  entrance  plasma  formed  when  the 
light  is  allowed  to  strike  the  surface  at  angles  other  than  normal  incidence.  The 
photographs  show  the  separation  of  the  plasma  into  different  spatial  components.  One 
component  extends  along  the  direction  of  the  light  beam;  this  apparently  is  the  short¬ 
lived  air  plasma.  Another  component  that  extends  normal  to  the  surface  can  be  seen; 
this  is  the  component  resulting  from  the  surface  explosion.  A  third  more  diffusely 
defined  plasma  component  can  be  seen,  extending  in  the  direction  of  the  light  reflected 
from  the  surface  and  is  presumably  additional  air  plasma  sustained  by  the  reflected 
light.  Figure  10  shows  a  t  i  me- i  n tegra ted  photograph  of  the  exit  plasma  formed  when  the 
light  leaves  at  an  angle  of  45°  to  the  surface.  This  photograph  shows  a  sinqle  plasma 
component  directed  normal  to  the  surface,  as  expected  if  a  surface  explosion  were  taking 
place. 

As  a  final  confirmation  of  assignments  for  the  different  plasma  components,  we 
looked  at  the  entrance  plasma  in  a  vacuum.  Figure  11  shows  photographs  of  the  entrance 
plasma  for  a  sample  tilted  26°  from  the  incident  light  both  at  1  atm  pressure  and  at 
about  1  Torr  pressure.  Because  of  the  shape  and  size  of  the  container,  it  was  impossible 
to  perform  the  experiment  at  angles  larger  than  26°.  As  a  result,  the  two  spatial  com¬ 
ponents  are  not  distinguishable  in  Figure  11(a)  and  (b).  However,  note  that  when  the 
air  is  removed  from  the  system  (Figure  11(c)  and  ( d) ) ,  the  plasma  component  which  extends 
along  the  light  beam  is  no  longer  seen;  the  plasma  visible  in  vacuum  Is  directed  normal 
to  the  surface,  i.e.,  only  the  explosion  plasma  remains. 

6.  Conclusions 

As  a  result  of  these  experiments ,  the  following  conclusions  are  reached.  First, 
the  surface  plasmas  are  a  result  rathe'  than  a  cause  of  the  surface  damage.  Second, 
two  kinds  of  plasmas  occur  at  the  entrance  surface,  one  which  is  sustained  oy  the  light 
beam  in  the  ambient  atmosphere  and  another  which  arises  from  a  micro-explosion  at  or 
near  the  surface.  The  exit  surface  plasma  is  of  the  explosion  type.  While  it  may  be 
true  that  some  of  the  features  seen  on  the  surface  of  a  material  that  has  been  subjected 
to  damaging  laser  radiation  might  be  caused  by  the  interaction  of  the  plasma  with  the 
surfice,  the  primary  cause  of  the  damaqe  is  not  the  plasma.  Rather,  the  surface  plasmas 
are  a  result  of  the  process  or  processes,  whatever  they  may  be,  that  qive  rise  to  the 
surface  damage. 

It  is  still  not  clear  why  the  exit  damage  is  physically  so  different  in  appearance 
than  entrance  damage.  If  the  mechanism  of  energy  deposition  were  simply  one  of  surface 
absorption,  It  would  be  expected  that  the  damage  would  be  very  similar  on  both  surfaces. 
For  some  reason,  there  seems  to  be  a  distinct  directionality  to  the  surface  damage 
mechanism.  Since  the  exit  surface  damage  threshold  appears  to  be  power  density  phenom¬ 
enon  and  since  it  is  not  apparently  dependent  on  whether  the  beam  is  converging  or 
diverging  before  it  reaches  the  surface,  one  would  be  inclined  to  conclude  that  the 
damage  mechanism  does  not  depend  strongly,  if  at  all,  on  phenomena  which  take  place  in 
the  bulk  of  the  material.  It  appears  to  be  strictly  a  surface  or  near-surface  phenom¬ 
enon.  The  possibility  that  a  directed  Intense  acoustic  wave  [1]  is  Involved  in  the 
damage  mechanism  is  still  open,  although  Stimulated  Brillouin  Scattering  (SBS)  is  not 
likely  since  there  seems  to  be  no  correlation  between  SBS  thresholds  and  damage 
thresholds  for  a  number  of  different  solids.  A  mechanisms  of  the  type  proposed  by 
Hellwarth  for  bulk  damage  [9]  might  be  a  valid  surface  damage  mechanism.  The  end  result 
of  this  mecharism  is  the  generation  of  intense  acoustic  phonons,  but  whether  one  would 
expect  the  acoustic  wave  to  be  directed  in  any  particular  way  is  not  known  at  this  time. 
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Fig.  9. 

Time  Integrated  Photographs  of  Entrance 
Plasma  at  Sapphire  Surface  Tilted  With 
Respect  to  Incident  Beam  Direction. 
Angle  Between  Light  Direction  and  Sur¬ 
face  Normal  (a)  45°,  (b)  60°. 
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Fig.  10. 

Time  Integrated  Photograph  of  Exit 
Plasma  at  Sapphire  Surface  Tilted 
45°  With  Respect  to  Light  Propaga¬ 
tion  Direction. 
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Fig.  11. 

Time  Integrated  Photographs  of  Entrance  Plasma  at 
Sapphire  Surface  Tilted  26°  From  Beam  Direction 
for  Different  Optical  Attenuations  at  Camera. 

(a)  and  (b)  Taken  at  1  atm  In  Air,  (c)  and  (d) 
Taken  at  1  Torr  Pressure.  Optical  Attenuation 
a)  N.D.  1.0,  (b)  N.D.  1.7,  (c)  N.D.  0.7, 

(d)  N.D.  1.7. 
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2.  Surface  Damage  Threshold  for  Ruby  as  a  Function  of  Beam 

Size  and  Divergence 

A  series  of  experiments  similar  to  those  performed  for 
sapphire  and  described  in  the  preceding  manuscript  were  carried  out 
on  a  ruby  sample.  Both  entrance  and  exit  surface  damage  thresholds 
were  measured  as  a  function  of  lens-to-sample  distance.  While  the 
entrance  damage  threshold  power  shows  a  similar  dependence  on 
distance  as  that  for  the  sapphire  sample  (see  Fig.  1-4),  the  exit 
threshold  shows  a  marked  difference  (Fig.  1-5).  This  phenomenon  is 
not  understood,  but  we  strongly  suspect  that  it  is  related  to  the 
intensity  dependent  beam  distortion  observed  in  ruby.  The  latter 
phenomenon  is  discussed  in  Semiannual  Report  3  and  summarized  in 
Section  I-D-2.  Since  the  beam  distortion  was  observed  to  occur  at 
powers  well  below  the  surface  damage  threshold  power,  it  is  reason¬ 
able  to  suspect  that  the  intensity  at  the  exit  surface  of  a  ruby  sample 
would  not  be  easy  to  predict.  The  fact  that  the  beam  profile  in  ruby 
changes  with  incident  power  and  position  indicates  that  the  power 
density  at  a  given  location  (e.  g. ,  the  exit  surface)  cannot  be  easily 
determined. 

3.  Morphology  of  Entrance  and  Exit  Damage 

Figures  1  and  2  of  the  preceding  manuscript  show  scanning 
electron  micrographs  of  entrance  and  exit  damage  sites  which  illus¬ 
trate  the  striking  difference  in  the  morphology  of  the  two  types  of 
damage.  At  the  time  the  paper  was  written  it  was  not  clear  why  the 
two  kinds  of  damage  have  such  a  different  appearance.  The  crazing, 
which  is  the  distinguishing  characteristic  of  entrance  damage,  is  not 
seen  at  the  exit  except  at  very  high  powers,  and  the  crater  type 
damage  that  occurs  at  exit  surfaces  is  not  seen  at  the  entrance  surface. 
Additional  pictures  of  the  different  features  of  surface  damage  are 
shown  in  a  paper  which  we  published  early  in  this  program/  Recently 
we  have  examined  the  morphology  of  entrance  surface  damage  that 
was  formed  when  the  sample  was  placed  in  a  vacuum  chamber  at 
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Fig.  1-4.  Beam  Area  and  Entrance  Surface 
Damage  Threshold  as  a  Function 
of  Surface  Distance  from  Beam 
Waist  for  19  cm  Lens . 
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~1  Torr  pressure.  Recall  that  the  entrance  surface  plasma  was  seen  to 
have  two  components,  an  air  plasma  and  an  explosion  plasma.  The  air 
plasma  is  absent  when  the  sample  is  damaged  in  a  vacuum.  Figure  1-6 
shows  optical  and  scanning  electron  micrographs  of  entrance  surface 
damage  generated  in  a  vacuum.  Note  that  no  crazing  exists  in  these 
pictures.  The  distinguishing  feature  is  now  a  slight  rippling  of  the 
surface.  Figure  1-7  shows  an  optical  micrograph  of  vacuum  entrance 
surface  damage  for  a  sample  tilted  45°  with  respect  to  the  incident 
light  beam.  Here  we  see  what  appears  to  be  only  a  slight  darkening  of 
the  surface  with  little  or  no  change  in  surface  topography.  Attempts 
to  obtain  scanning  electron  micrographs  of  this  damage  site  have  so 
far  failed,  perhaps  because  there  is  too  small  a  change  in  surface 
features  to  be  detected. 

From  these  observations  we  conclude  that  the  crazing  which 
characterizes  entrance  surface  damage  arises  from  an  interaction  of 
the  air  plasma  with  the  surface. 
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Fig.  1-6.  Optical  and  Scanning  Electron 
Micrographs  of  Vacuum  Entrance 
Surface  Damage  for  Beam  at 
Normal  Incidence. 
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Fig.  1-7.  Optical  Micrograph  of  Vacuum 
•  Entrance  Surface  Damage  for- Beam 

Incident  at  459. 
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4. 


Attempts  to  Increase  Surface  Damage  Threshold  by  Ion  Beam 
Polishing 


During  the  latter  stages  of  this  program  we  have  obtained  some 
very  encouraging  results  which  indicate  that  surface  damage  thresholds 
can  be  substantially  improved.  We  carried  out  a  few  experiments  in 
which  sapphire  samples  were  polished  using  an  energetic  argon  ion 
beam  and  found  that  the  surface  damage  threshold  for  both  entrance  and 
exit  were  increased  substantially  (up  to  ten  times  depending  on  the 
particular  points  compared).  This  work  led  to  a  publication  in  Applied 
Physics  Letters  (1  July  1972)  and  a  presentation  at  the  Fourth  ASTM: 
NBS  Conference  on  Damage  in  Laser  Materials,  Boulder,  Colorado, 
June  1972.  A  reprint  of  this  manuscript,  which  gives  the  essential 
features  of  this  work,  follows. 


Later-induced  damage  in  transparent  dielectrics:  ion  beam  polishing 
as  a  means  of  increasing  surface  damage  thresholds* 

Concetto  R.  Giuliano 

Hugh*  Rttmrch  Laboratoria*,  Malibu.  California  90265 
(Received  14  April  1972) 

PolUhlng  of  sspphlre  surface*  with  energetic  Ar*  ion  beams  ie  shown  to  result  In  a  substan¬ 
tial  Increase  In  laser  damage  threshold  over  that  for  conventionally  polished  surfaces.  Data 
for  both  entrance  and  exit  damage  are  presented.  The  results  are  interpreted  In  terms  of  an 
Increase  In  surface  strength  with  ion  beam  polishing. 


It  is  a  generally  accepted  fact  that  the  threshold  for  ated.  Any  foreign  material  on  a  surface  subjected  to 

laser- Induced  surface  damage  of  transparent  dielectrics  high  optical  fluxes  can  cause  catastrophic  damage  to 

is  a  function  of  the  condition  of  the  surface  being  irradl-  occur  at  levels  well  below  those  which  can  be  sustained 
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G.  1.  Number  of  damage  aitea  aa  a  function  rt  threahold 
wer  denaity  at  entrance  and  exit  aurfacea  of  sapphire  cry  a 
a  for  abraalvely  pollahed  and  lon-beam-pollahed  aamplea. 


clean  aurfacea.  The  quality  of  aurface  fabrication  has 
so  been  recognized  as  being  important  in  influencing 
a  damage  resistance  of  materials.  For  example,  it  is 
own  that  a  frosted  surface  damages  more  easily  than 
surface  with  a  good  optical  polish. 1  Various  chemical 
rface  treatments  have  been  tried  on  different  mate- 
ila,  *•'  some  of  which  have  shown  improvement,  but 
en  the  effect  is  only  temporary. 

t  have  found  that  polishing  of  sapphire  surfaces  with 
i  beams  results  in  a  marked  increase  of  surface 
mage  threahold  compared  with  that  for  abrasively 
ished  samples. 

perlments  were  carried  out  using  a  single  transverse* 
I  longitudinal- mode  ruby  laser  whose  output  is  to- 
led  on  either  the  entrance  or  exit  surface  of  a  sap- 
re  sample.  The  samples  are  sapphire  bars  3  in.  long 
h  0. 25-in.  -square  cross  section.  The  laser  and  asso- 
ted  monitoring  apparatus  are  described  in  detail 
e where. 4  Typical  laser  pulse  width  and  peak  power 
1 20  nsec  (FWHM)  and  1  MW,  respectively.  The  area 
he  focused  spot  has  been  measured  using  a  modified 
tlpe- lens—  camera  technique'  to  be  10*4  cm*  at  the  1/ 
□tensity  points.  The  beam  profile  was  measured  to 
jausslan  down  to  8%  of  the  peak.  The  power  incident 
he  sample  is  varied  by  rotating  the  first  of  a  pair  of 
•spaced  polarizers,  the  second  being  held  at  a  fixed 
mtatlon 

rpical  series  of  measurements  is  carried  out  as 
?ws.  The  laser  is  fired  at  a  particular  spot  on  the 
red  surface  at  a  power  below  damage  threshold.  The 
ipl  is  then  examined  through  a  traveling  microscope 
tout  disturbing  the  sample  between  shots.  If  no  dam¬ 


age  Is  observed  the  incident  power  is  increased  (by 
- 10  —  20%  and  again  the  sample  is  examined.  The  pro¬ 
cedure  continues  until  damage  is  observed;  then  the 
sample  is  moved  and  the  process  repeated  on  an  undam¬ 
aged  spot.  Typically  10—20  damage  thresholds  are  mea¬ 
sured  in  this  way  for  a  given  surface. 

Abrasively  polished  samples  were  used  as  obtained 
from  the  manufacturer  (Union  Carbide).  They  were 
fabricated  with  the  “laser  finish”  specifications  typical 
for  ruby  laser  rods.  Ion-polished  samples  were  ex¬ 
posed  after  mechanical  processing  to  an  Ar'  beam  (7  kV, 
300  pA/cm*)  for  periods  from  2  to  4  h.  The  beam 
strikes  the  surface  at  an  angle  of  20c  from  the  surface 
plane.  An  estimated  200  A/min  Is  removed  from  the 
sapphire  surface  under  these  conditions,  giving  a  total 
of  from  about  2. 5  to  5  m  of  material  removed  In  a  given 
exposure. 

The  results  ot  a  series  of  measurements  on  abrasively 
polished  and  ion-beam-polished  samples  are  shown  in 


Imm— H  H— 
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FIG.  2.  Examples  of  x-ray  diffraction  topographs  taken  for  (a) 
abrasively  polished  sample,  <b)  Ar*  ion  polished  for  2  h  (7  kV, 
300  S A/ cm*),  (c)  same  conditions  as  (b),  but  for  4  h. 
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FI*.  I.  In  Pig*.  1(a)  and  1(b)  th*  result*  lor  entrance 
surface  damage  are  shown.  Por  the  abrasively  polished 
sample  most  of  the  damage  threshold*  occur  between  I 
and  2. 5  GW/cm*  with  a  few  higher  values  occurring.  For 
the  ion-polished  case  [Pig  1(b)]  we  see  a  few  damage 
thresholds  ranging  from  2  to  6  GW/cm1  with  a  large 
fraction  (>90%  above  10  GW/cm*.  [The  dashed  portion 
of  Fig.  1(b)  indicates  that  we  were  unable  to  reach  dam¬ 
age  threshold  with  the  maximum  output  from  our  laser 
under  the  focusing  conditions  of  our  experiment.  ]  Figure 
1(c)  shows  the  exit  damage  thresholds  for  the  abrasively 
polished  sample.  Here  the  bulk  of  the  points  are  clus¬ 
tered  between  1  and  2  GW/cm*.  Foi  the  ion-polished 
case  [Pig.  1(d)]  the  values  are  again  spread  over  a  wide 
range  (up  to  0  GW/cm*),  but  for  this  case  no  damage 
thresholds  lower  than  2  GW/cm*  are  observed .  The  data 
in  Pig.  1(d)  were  obtained  from  samples  which  were  ion 
polished  for  2  and  4  h.  No  noticeable  change  in  the  dis¬ 
tribution  was  observed  between  the  2-  and  4-h  cases, 
and  the  data  are  presented  together  in  the  figure. 

Figure  2  shows  examples  of  x-ray  reflection  topographs 
for  sapphire  samples  which  were  abrasively  polished 
[Fig.  2(a)]  and  ion  beam  polished  [Figs.  2(b)  and  2(c)]. 
The  depth  of  the  region  being  probed  in  these  topographs 
depends  on  the  angle  of  incidence  and  energy  of  the  x 
rays  as  well  as  other  experimental  conditions  and  sam¬ 
ple  properties.  An  estimated  region  10—20  jx  thick  is 
being  probed  in  the  above  topographs.  We  see  in  Fig. 

2(a)  a  number  of  deviations  from  crystallinity  which  ap¬ 
pear  to  follow  scratch  lines  from  the  mechanical  polish¬ 
ing  process.  In  Fig.  2(b)  the  topograph  for  a  different 
sample  shows  that  some  of  the  mechanical  damage  ap¬ 
pears  to  have  been  “erased"  after  2  h  of  ion  beam 
polishing.  Figure  2(c)  shows  a  topograph  taken  after  4 
h  of  exposure  to  the  ion  beam,  and,  although  evidence  of 
damage  still  is  seen,  there  is  a  marked  difference  be¬ 
tween  this  and  the  previous  topographs.  We  present 
these  topographs  to  point  out  that  a  distinct  change  in 
surface  topography  is  brought  about  by  ion  beam  polish¬ 
ing.  A  detailed  correlation  between  the  changes  in  the 
topographs  and  the  damage  thresholds  is  not  possible  at 
this  time. 

The  main  conclusion  reached  by  these  preliminary  re¬ 
sults  is  that  a  distinct  increase  (-2—6  times)  in  surface 
damage  threshold  can  be  realized  by  Ion  beam  polishing. 
Although  there  is  a  wide  scatter  of  results  in  measured 
thresholds  for  ion-polished  samples,  we  see  that  even 
the  low  values  for  the  ion-polished  samples  are  general¬ 
ly  higher  than  those  for  the  abrasively  polished  samples. 
We  interpret  the  scatter  in  the  threshold  data  for  the 
ion-polished  samples  to  indicate  a  variation  of  surface 
strength  from  place  to  (dace  over  the  sample.  R  is  very 
likely  that  the  influence  of  mechanical  damage  due  to 
abrasive  polishing,  although  somewhat  lessened,  has 
not  been  completely  removed  under  the  conditions  of 


our  experiments.  This  is  suggested  by  th*  x-ray  topo¬ 
graphs,  although  as  yet  no  definite  correlation  between 
topograph  details  and  threshold  variations  has  been  es¬ 
tablished.  W*  choose  to  interpret  the  data  as  reflecting 
a  variation  of  threshold  fronf  spot  to  spot  on  the  surface 
rather  than  in  a  probabilistic  manner  as  described  by 
Bass  and  Barrett.'  It  is  apparent  from  the  topographs 
that  the  surface  crystallinity  varies  considerably  from 
place  to  place  over  th*  sample  and  that  the  surface  is  not 
uniform  all  over.  A  possible  Interpretation  of  the  differ¬ 
ence  between  the  abrasive-polish  and  ion-polish  results 
is  that  In  the  former  case  the  surface  is  “uniformly 
weak*’  and  the  damage  thresholds  are  found  to  have  a 
relatively  narrow  range  of  values,  whereas  in  the  latter 
case  a  distinct  toughening  of  the  surface  has  taken  place, 
but  the  Influence  of  mechanical  damage  from  the  abra¬ 
sive  polishing  still  exists. 

A  final  point  worthy  of  mention  is  that  the  damage,  when 
it  occurs  at  a  high  threshold,  is  generally  much  more 
extensive  than  when  it  occurs  at  a  low  threshold.  For 
example,  exit  surface  damage  in  abrasively  polished 
samples  is  generally  characterized  by  a  small  fractured 
crater  or  pit  about  10  v  in  diameter  when  it  occurs  near 
threshold.  On  the  other  hand,  the  high-threshold  damage 
seen  in  the  ion-polished  samples  is  more  severe;  the 
damage  craters  are  usually  several  times  larger  and 
more  extensively  fractured.  This  fact  suggests  that  the 
effect  of  ion  polishing  is  to  toughen  or  strengthen  the 
surface  in  some  way  rather  than  to  change  the  basic 
damage  mechanism,  whatever  it  may  be.  It  seems  that 
the  surface  is  able  to  sustain  a  larger  deposition  of  en¬ 
ergy  before  catastrophic  breakdown  occurs,  but,  when 
the  surface  does  let  go,  it  goes  with  a  bigger  bang. 

The  author  wishes  to  acknowledge  critical  suggestions  by 
V.  Evtuhov,  as  well  as  the  aid  of  H.  L.  Garvin,  who  per¬ 
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provided  the  x-ray  topographs. 
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5. 


CommentB  on  Fluctuations  of  Surface  Damage  Threshold  Values 


We  commented  in  the  preceding  manuscript  that  we  chose  to 
interpret  the  large  scatter  in  our  threshold  results  as  reflecting  a  real 
variation  in  damage  resistance  from  spot-to-spot  on  the  surface  rather 
than  reflecting  an  intrinsic  damage  probability.  To  support  this  point 
further,  we  will  describe  the  experimental  technique  employed  in 
obtaining  the  threshold  values. 

The  site  of  interest  on  the  sample  surface  is  first  irradiated 
with  the  laser  at  a  level  well  below  that  believed  necessary  to  cause 
damage,  and  then  the  sample  is  examined  carefully  for  the  presence  of 
damage.  If  no  damage  is  observed  the  power  is  increased  (by  10  to 
20%)  and  another  shot  is  fired  at  the  same  spot.  The  sample  is  then 
reexamined  and  the  process  continues  until  damage  is  finally  observed 
(or  until,  as  in  some  cases  of  entrance  damage,  the  maximum  power 
available  from  our  laser  is  reached).  In  a  number  of  cases  when  the 
spot  of  interest  was  not  damaged  at  a  relatively  high  power  (say  twice 
the  value  for  some  of  the  lower  thresholds  observed  for  other  sites), 
we  held  the  power  essentially  constant  and  fired  the  laser  for  as  many 
as  fifteen  times  before  continuing  to  increase  the  power  further. 

Hence,  a  number  of  sites  were  subjected  to  up  to  20  nondamaging  shots 
whose  peak  powers  lay  in  the  range  from  70  to  90%  of  the  power  at 
which  damage  was  finally  seen  to  occur.  This  type  of  behavior  strongly 
supports  the  conclusion  that  a  marked  variation  in  damage  resistance 
exists  from  spot-to-spot  on  the  surface  and  that  the  surface  cannot  be 
considered  to  be  uniform.  The  smaller  degree  of  scatter  for  the 
abrasively  polished  samples  suggests  that  the  surface  is  more  uniform 
(in  this  case,  uniformly  bad)  than  in  the  ion  polished  samples,  in  which 
the  effects  of  mechanical  damage  due  to  the  influence  of  the  abrasive 
seem  to  have  been  removed  in  some  spots  but  not  in  others. 
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D.  BEAM  PROFILE  AND  SPOT  SIZE  MEASUREMENTS 
I .  Summary  of  Early  Work 

In  the  early  stages  of  this  program  the  damage  threshold  power 
densities  quoted  were  based  on  beam  sizes  calculated  according  to  the 
following  assumptions: 

•  The  spatial  beam  profile  at  the  laser  output  mirror 
is  gaussian  with  a  1  mm  radius  (1/e  radius  for  the 
electric  field). 

•  The  beam  diffracts  according  to  the  well-known 
expressions  for  gaussian  beam  propagation. 

•  The  beam  suffers  no  distortion,  focusing,  or 
defocusing  on  passing  through  the  ruby  amplifier; 
hence,  the  diffraction- limited  beam  waist  occurs 
at  the  location  calculated  according  to  the  propaga¬ 
tion  equations. 

During  the  period  discussed  in  Semiannual  Report  3,  we  looked 
more  carefully  into  the  validity  of  these  assumptions;  these  results 
will  be  summarized  in  this  report. 

We  undertook  a  series  of  measurements  to  determine  the 
spatial  beam  profile  at  different  locations  and  under  a  number  of 
different  conditions.  Early  experiments  employed  the  measurement 
of  a  pair  of  spots  with  a  known  relative  attenuation  for  both  near  and 
far-field  spot  sizes.  As  mentioned  in  the  previous  report,  we  obtained 
a  near-field  spot  size  of  1  mm  (M5%  accuracy)  by  photographing  a 
pair  of  spots  on  a  ground  glass  screen  and  also  by  measuring  beam 
patterns  on  exposed  Polaroid  film  with  the  relative  energies  of  succes¬ 
sive  shots  known.  Far-field  divergence  measurements  of  the  unampli¬ 
fied  oscillator  using  the  two  spot  technique  gave  a  divergence  of 
0.  3  ±0.05  mrad  (half-angle).  These  results  are  based  on  the  assumption 
that  the  beam  profile  is  gaussian. 

For  the  more  detailed  studies,  we  attempted  to  determine  how 
close  to  gaussian  the  beam  profile  actually  is.  To  accomplish  this  we 
have  used  a  multiple  lens  camera  and  the  technique  described  by 
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Winer  .  The  setup  is  illustrated  in  Fig.  1-8.  We  image  a  particular 
plane  (plane  A)  onto  a  MgO  block  with  the  imaging  lens  choosing  the 
object  and  image  distances  to  get  convenient  magnification,  and  photo¬ 
graph  the  spot  on  the  MgO  block  with  the  multiple  lens  camera.  (This 
is  accomplished  by  placing  a  pinhole  of  known  size  at  plane  A,  imaging 
it  on  the  MgO  block  with  appropriate  magnification,  and  photographing 
it.  )  The  camera  has  nine  lenses  (f  =  ~11  cm,  12  mm  diameter),  each 
of  which  is  backed  by  a  different  calibrated  neutral  density  (N.  D.  )  filter. 
Thus  we  obtain  nine  spots  on  the  focal  plane  of  the  camera  with  known 
relative  exposures. 

The  principle  of  the  multiple  lens  camera  technique  for  meas¬ 
uring  intensity  profiles  is  that  the  film  response  need  not  be  known. 
Densitometer  scans  of  the  different  spots  are  made,  and  the  widths  of 
the  scans  are  measured  at  a  constant  density  determined  by  the  peak  of 
the  curve  vdth  the  least  transmission. 

Beam  profiles  and  spot  sizes  were  photographed  using  the  mul¬ 
tiple  exposure  camera  at  different  places  and  under  different  conditions 
of  amplifier  pumping.  For  example,  we  measured  the  beam  size  at 
the  laser  output  resona.it  reflector,  at  the  entrance  and  exit  of  the 
amplifier,  at  the  focusing  lens  position,  and  at  a  number  of  locations 
beyond  the  focusing  lens,  including  the  focal  plane. 

The  details  and  results  of  these  measurements  are  presented 
in  Semiannual  Report  3.  The  following  summarizes  the  highlights. 

•  Although  the  beam  profile  is  very  close  to  a  gaussian 
distribution,  the  far-field  divergence  is  substantially 
larger  than  that  calculated  for  a  gaussian  beam 

(0.  35  mrad  half-angle  compared  with  a  calculated 
value  of  0.  20  mrad). 

•  The  measured  beam  waist  spot  size  is  also  sub¬ 
stantially  larger  than  that  calculated  using  the 
gaussian  mode  propagation  equations.  The  dis¬ 
crepancy  between  measured  and  calculated  values 
depends  on  the  focal  length  of  the  lens.  Because 

of  these  two  observations,  it  is  clear  that  calculated 
spot  sizes  are  suspect  even  for  beams  which  are 
very  close  to  gaussian. 
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Schematic  Representation  of  Setup  Used  In  Beam 
Profile  Measurements. 


•  The  ruby  amplifier  acts  as  a  negative  lens  whose 
focal  length  depends  on  optical  pumping.  As  a 
result,  the  position  of  the  beam  focus  varies  as 
the  amplifier  pumping  varies.  After  having  dis¬ 
covered  this  fact,  a  particular  set  of  pumping 
conditions  were  chosen,  and  all  subsequent  exper¬ 
iments  were  carried  out  with  those  conditions 
fixed. 

2.  Spot  Size  Measurements  Inside  Samples 

In  this  section  we  discuss  the  results  of  experiments  where  we 
measured  the  spot  size  inside  both  ruby  and  sapphire  samples  by  the 
same  techniques  as  that  described  in  Section  I-D-l.  The  obvious 
extension  of  the  previous  spot  size  measurements  in  air  would  be  to 
continuously  monitor  from  shot-to-shot  the  beam  width  inside  a  sample 
while  performing  the  damage  experiments.  In  this  way  ve  would  have 
a  direct  measure  of  energy  density  rather  than  having  to  infer  it  from 
previous  measurements.  All  the  experiments  described  here  were 
carried  out  using  the  19  cm  lens  focused  inside  the  particular  sample. 
The  pumping  conditions  in  the  amplifier  were  fixed  at  150  pF  and 
7.  5  kV,  and  the  light  intensity  incident  on  the  sample  was  varied  by  a 
pair  of  Glan-Kappa  prisms,  the  first  of  which  was  rotated  to  the 
desired  angle  and  the  second  of  which  was  fixed  so  that  the  polarization 
was  always  the  same  inside  the  sample  (E  perpendicular  to  the  C-axis). 

The  lenses  and  samples  were  placed  so  that  the  waist  of  the 
focused  light  beam  would  occur  about  2  cm  inside  the  exit  surface  of 
the  sample.  The  plane  being  imaged  in  the  photographs  lies  0.  5  cm 
upstream  from  the  beam  waist. 

A  number  of  photographs  taken  with  the  multiple  lens  camera 
are  shown  in  Figs.  1-9  and  I- 10.  The  interesting  result  of  this  probing 
is  that  the  beam  profile  begins  to  change  radically  in  ruby  at  powers 
well  below  the  damage  threshold;  whereas  for  sapphire,  the  beam 
profile  is  essentially  smooth  all  the  way  to  and  beyond  the  threshold 
for  bulk  damage. 

We  wish  to  emphasize  here  that  all  the  photographs  were  taken 
in  the  same  plane  in  the  sample,  and  the  only  parameter  that  was 
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Fig.  1-9. 

Multiple  Lens  Camera  Photo¬ 
graphs  of  Beam  Profile  Inside 
Ruby  Sample  for  Different 
Incident  Energies  and  Arbi¬ 
trary  Relative  Exposures. 
(Constant  amplifier  pumping  - 
150  yF,  7.5  kV.) 
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Fig.  1-10. 

Multiple  Lens  Camera  Photographs  of  Beam 
Profile  Inside  Sapphire  Sample  for  Dif¬ 
ferent  Incident  Energies  and  Arbitrary 
Rel atl ve  Exposures  • 


varied  wai  the  amount  of  energy  incident  on  the  )9  cm  focusing  lens. 

We  see  from  Fig.  I-9(a)  that  at  low  incident  energy  (0.  33  mJ)  the  beam 
profile  is  smooth,  but  that  by  the  time  we  reach  3  mJ  (Fig.  I- 9(b) )  we 
begin  to  see  what  appears  to  be  a  bright  central  spot  with  a  less  bright 
peripheral  halo.  (In  these  experiments,  the  laser  pulse  width  was 
about  30  nsec,  FWHM.  )  At  higher  incident  energies  we  see  a  variety 
of  patterns  in  the  beam  profile.  We  note  also  that  the  threshold  for 
internal  damage  in  this  particular  sample  measured  under  the  same 
conditions  is  30  mJ  incident  on  the  focusing  lens,  and  that  the  threshold 
for  exit  surface  damage  is  about  10  mJ.  Thus  we  see  these  effects  of 
beam  distortion  well  below  the  threshold  for  any  catastrophic  phenomena. 

The  whole  question  of  measurement  of  energy  density  inside  ruby 
samples  by  measuring  spot  sizes  assumes  a  completely  different  aspect 
in  the  light  of  these  drastic  changes  in  the  beam  profile.  Obviously, 
it  would  be  extremely  difficult  to  determine  energy  densities  in  these 
peculiar  rpatial  distributions. 

When  we  examine  sapphire  under  the  same  conditions  as  those 
described  for  ruby,  we  see  no  such  behavior,  f  igure  I- 10  shows 
photographs  taken  under  two  extreme  conditions.  Figure  I-  10(a)  was 
taken  at  low  incident  energy  (M).  3  mJ)  and  Fig.  I- 10(b)  at  an  energy 
higher  than  the  bulk  damage  threshold,  which  was  20  mJ  for  this  sam¬ 
ple.  For  a  large  number  of  shots  over  this  range  of  energies  we  saw 
no  irregularities  in  the  beam  profile.  After  reaching  threshold  for 
internal  damage,  we  saw  a  general  decrease  in  light  reaching  the 
imaging  optics  presumably  because  of  scattering  from  the  damage  sites 
and  increased  beam  spreading  due  to  self-focusing. 

Again  we  see  an  interesting  difference  in  the  behavior  of  rub) 
as  compared  with  sapphire  when  subjected  to  intensi/e  illumination 
at  6943  A.  It  is  reasonable  to  suspect  that  the  beam  distortion  seen 
in  ruby  arises  from  its  absorption  at  6943  A. 


We  do  not  have  sufficient  data  to  locate  the  threshold  for  this  beam 
distortion  effect.  We  only  know  it  lies  between  0.  3  and  3  mJ  for  the 
conditions  described  above. 


Because  the  obvious  difference  between  sapphire  and  ruby  at 
6943  A  is  the  optical  absorption,  we  decided  to  look  at  the  beam  profile 
inside  pumped  ruby  samples  to  see  whether  the  same  distortion  is  seen 
under  these  circumstances.  We  chose  an  incident  power  level  where 
the  distortion  is  pronounced  in  an  unpumped  sample  and  optically 
pumped  the  sample  at  different  levels.  The  results  were  inconclusive, 
because  while  the  beam  profile  definitely  changed  with  optical  pumping, 
there  was  no  marked  improvement  or  lessening  of  the  distortion  that  could 
be  directly  attributed  to  the  pumping.  At  the  time  these  experiments 
u<»re  carried  out,  we  were  having  difficulty  with  oscillator  and  ampli¬ 
fier  alignment  fluctuations,  which  also  gave  rise  to  distinct  beam  pro¬ 
file  irregularities.  It  was  not  possible  at  the  time  to  separate  effects 
arising  from  a  possible  misalignment  variation  from  shot-to-shot. 

This  lack  of  reproducibility  of  results  from  shot-to-shot  and  from 
day-to-day  made  it  impossible  to  reach  any  definitive  conclusion  con¬ 
cerning  the  effects  of  optical  pumping  on  the  beam  distortion  in  ruby. 

At  this  time  the  question  still  remains  unanswered. 

3.  Ring  Structure  in  the  Focused  Single  Mode  Laser  Beam  Beyond 
Beam  Waist  ~ 


Our  plans  to  examine  in  more  detail  the  beam  distortion  effect 
in  optically  pumped  ruby  samples  were  preceded  by  preliminary 
experiments  in  which  the  ruby  sample  was  absent  from  the  beam. 
During  these  experiments  we  obtained  an  unexpected  result  which 
diverted  us  from  our  original  objective  and  which  we  pursued  further 
to  understand  its  origin.  We  will  describe  the  results  of  these  obser¬ 
vations  in  this  section. 

We  wish  to  point  oui  that  these  phenomena  were  observed  with 
only  the  light  from  the  laser.  No  other  elements  are  placed  in  the 
beam  between  the  oscillator  output  mirror  and  the  lenses  except  for 
beam  splitters  and  polarizer  attenuators;  the  removal  of  these  does 
not  change  what  was  observed.  Also,  the  presence  or  absence  of  the 
amplifier  does  not  affect  what  is  observed. 
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The  essential  result  of  the  observations  is  shown  in  Fig.  I- 11, 
in  which  we  see  the  following  features.  The  beam  profile  is  smooth 
and  approximately  gaussian  in  the  back  focal  plane  of  the  lens  and  also 
at  the  bean,  waist. 

At  a  point  somewhat  beyond  the  beam  minimum,  the  beam  begins 
to  show  a  marked  deviation  from  a  smooth  profile.  This  first  shows  up 
as  a  bright  spot  in  the  center  of  the  beam.  A  short  distance  farther 
down  we  encounter  a  minimum  on-axis  with  a  bright  ring  around  it.  As 
we  progress  farther  downstream  we  see  more  rings  developing  in  the 
beam  profile. 

At  first  we  suspected  that  these  peculiar  beam  profiles  might  be 
an  artifact  of  some  sort,  and  attempted  to  examine  more  closely  the 
conditions  under  which  they  are  observed  (i.  e. ,  we  wanted  to  determine 
if  the  phenomenon  was  real). 

We  found  that  the  same  kind  of  behavior  was  observed  for  differ¬ 
ent  focusing  lenses  and  different  imaging  lenses,  showing  that  our 
observation  was  not  a  peculiarity  of  a  particular  lens.  The  pattern  of 
rings,  etc.,  is  reminiscent  c*  what  might  be  expected  if  spherical 
aberration  were  present.  To  check  whether  the  phenomenon  was  a 
property  of  the  lenses  or  the  ruby  laser  itself,  we  passed  the  light 
from  the  He-Ne  alignment  laser  through  the  same  setup.  A  number  of 
regions  along  the  beam  were  imaged  onto  the  MgO  block  and  photo¬ 
graphed.  We  saw  no  evidence  of  rings  or  other  irregularities  over  the 
entire  range  of  positions  on  both  sides  of  the  beam  focus.  This  led  us 
to  conclude  that  what  we  observed  was  a  property  of  the  ruby  laser 
itself  and  not  of  the  lenses. 

Further  evidence  supporting  the  contention  that  the  structure  i& 
a  property,  at  least  in  part,  of  the  ruby  laser  itself  and  not  an  artifact 
of  the  system  of  observation  is  as  follows. 

•  he  structure  is  observed  with  different  focusing  and 

iging  lenses. 

•  The  structure  is  not  simply  a  propagation  phenomenon. 

To  verify  the  latter  point,  we  imaged  the  same  plane  with  lenses 
of  different  focal  length.  Therefore  the  distance  between  focusing  and 
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imaging  lenses  is  different  in  each  case.  The  pattern  observed  is  the 
same  in  the  two  experiments,  indicating  that  the  imaging  lenses  arc 
indeed  imaging  information  from  the  plane  of  interest.  Another  similar 
experiment  was  performed  in  which  the  distance  between  the  two  lenses 
was  kept  constant  and  the  MgO  block  was  moved  to  different  positions. 
Here  we  were  imaging  different  object  planes  at  different  magnifications. 
Patterns  characteristic  of  the  planes  of  interest  when  viewed  at  the 
same  magnification  were  observed.  The  pattern  changed  when  the  MgO 
block  was  moved  from  one  position  to  the  other,  and  the  corresponding 
change  was  the  same  as  if  the  new  plane  was  imaged  by  moving  the 
imaging  lens  rather  than  moving  the  MgO  block  alone  and  keeping  the 
lens-to-lens  distance  fixed. 

With  this  result  and  others  discussed  in  Semiannual  Report  4, 
we  were  convinced  that  the  phenomenon  is  real  and  a  property  of  the 
ruby  laser  or  perhaps  a  combined  property  of  the  laser  and  the  lenses. 
We  took  a  few  beam  profile  pictures  of  the  laser  output  with  no  lenses 
in  the  beam  at  a  number  of  distances  from  the  laser  from  2  to  55  m; 
no  evidence  of  ring  structure  in  the  beam  profile  was  seen. 

The  behavior  described  above  is  reminiscent  of  the  description 

4 

given  by  Innes  and  Bloom  who  discuss  beam  profiles  for  focused 
gaussian  beams  for  different  amounts  of  truncation  and  spherical 
aberration.  However,  the  effects  they  predict  become  significant  only 
for  large  spherical  aberration  (X./2)  and  relatively  severe  truncation 
of  the  gaussian  beam.  Of  course,  we  have  no  assurance  that  the  phase 
fronts  of  our  beam  are  uniform  at  the  output  mirror  of  the  oscillator. 
Nonuniform  phase  across  the  beam  could  certainly  give  rise  to  irreg¬ 
ularities  in  the  beam  profile  such  as  that  observed.  It  is  interesting 
to  note  that  this  behavior  could  well  be  a  general  characteristic  of  all 
single  mode,  solid  state  lasers.  Further  pursuit  of  this  phenomenon 
would  certainly  be  valuable  but  is  probably  outside  the  realm  of  the 
present  program. 

We  have  recently  become  aware  of  a  similar  observation  by 
workers^  at  the  University  of  Southern  California  who  observed  fluc¬ 
tuations  in  the -on-axis  intensity  of  a  focused  ruby  laser  beam  at  points 
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downstream  from  the  beam  waist.  They  also  pointed  out  that  computer 
solutions  for  focused  truncated  gaussian  beams  showed  that  such  fluc¬ 
tuations  are  to  be  expected  and  result  from  the  truncation.  Since  the 
phenomenon  occurs  much  closer  to  the  beam  waist  for  short  focal 
length  lenses  than  for  long  focal  length  lenses,  it  is  important  that  the 
precise  location  of  the  beam  waist  be  known  for  experiments  where  the 
beam  is  focused  on  a  surface,  because  small  fluctuations  in  sample-to- 
lens  distance  could  give  rise  to  a  large  variation  in  beam  properties 
and  damage  thresholds. 

E.  EXPERIMENTAL  APPARATUS  USED  IN  DAMAGE  THRESHOLD 

STUDIES 

The  early  stages  of  the  experimental  phase  of  this  program  were 
devoted  to  improving  and  refining  the  performance  of  the  laser  source 
used  in  the  program.  The  setup  showing  the  laser,  amplifier  and 
other  experimental  apparatus  is  shown  in  Fig.  1-12.  The  single  longi¬ 
tudinal  and  transverse  mode  oscillator  employs  a  4  in.  long  by  1/4  in. 
diameter  ruby  pumped  by  two  linear  lamps  in  a  double  elliptical  pump 
cavity.  The  ruby  crystal  is  water  cooled  by  a  closed  cycle  refrigeration 
system  maintained  at  0°C.  The  high  reflectivity  mirror  is  coated  with 
a  99+%  reflectivity  high  field  damage  coating  from  Perkin  Elmer 
Corporation.  The  original  Q- switch,  a  rotating  Brewster  angle  prism, 
was  replaced  by  a  solution  of  cryptocyanine  in  methanol  in  a  1  mm  path 
length  cell  whose  transmission  is '>-30%  at  6943  A.  The  2  mm  diameter 
aperture  allows  oscillation  in  the  TEMqq  mode. 

The  temperature  controlled  (34°C)  resonant  reflector  that  was 
designed  to  optimize  longitudinal  mode  control  consists  of  two  quartz 
etalons  and  a  quartz  spacer  whose  combined  effect  is  to  enhance  cavity 
modes  separated  by  2  cm  *  and  to  discriminate  against  intermediate 
modes. 

Portions  of  the  laser  beam  are  split  off  in  various  ways  (see 
Fig.  1-12),  so  that  the  power  output,  near  and  far-field  patterns,  and 
Fabry-Perot  patterns  can  be  monitored  for  each  shot.  This  is  accom¬ 
plished  in  the  following  way.  Light  reflecting  from  wedged  beamsplitter 
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Fig.  1-12.  Schematic  Representation  of  Experimental  Apparatus. 
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Wj  gives  two  diverging  beams;  one  of  them  hits  the  magnesium  oxide 
diffuser,  where  the  scattered  light  is  monitored  by  a  biplanar  photodiode 
used  as  our  power  monitor.  The  second  beam  from  Wj  hits  ground 
glass  screen  G,  where  it  is  photographed  through  lens  L  and  the  1  m 
focal  length  camera  focused  at  infinity.  This  gives  a  magnified  (— 5x) 
near-field  picture.  Another  portion  of  the  light  is  removed  by  beam¬ 
splitter  Bj  and  hits  mirror  M^,  which  can  be  placed  in  or  out  of  posi¬ 
tion  depending  on  the  use  of  the  alignment  laser.  From  the  light 
either  goes  to  the  Fabry-Perot  interferometer  or  can  be  partially 
reflected  from  wedged  beamsplitter  W^,  where  it  results  in  a  pair  of 
far-field  patterns.  A  0.  6  neutral  density  filter  is  placed  near  the  focal 
plane  of  the  camera  so  that  the  far-field  pattern  and  the  Fabry-Perot 
pattern  can  be  seen  at  two  different  exposures.  The  two  Gian  prisms 
are  used  as  a  variable  attenuator  after  the  amplifier.  Beamsplitter  B^ 
samples  the  light  to  photodiode  No.  2,  which  monitors  the  power  incident 
upon  the  focusing  lens  which  was  designed  for  minimum  spherical 
aberration  (Special  Optics).  The  signal  from  the  detector  is  integrated 
and  displayed  on  an  oscilloscope. 

The  water  cooled  amplifier  ruby  is  6  in.  long  by  0.  5  in.  diameter, 
with  one  end  wedged  relative  to  the  other  by  about  0.  5°.  The  input  end  of 
the  amplifier  rod  is  anti  reflection  coated  to  minimize  the  chances  of 
oscillation  within  the  amplifier  itself.  The  ruby  rod  is  closely  coupled 
to  a  helical  flashlamp,  which  is  pumped  with  a  power  supply  capable  of 
delivering  8  kJ  in  a  3  msec  pulse.  The  power  supply  employs  a  pulse 
shaping  network  of  20  sections,  each  section  pumping  for  150  p.sec.  The 
maximum  gain  obtained  with  the  amplifier  is  about  10  dB. 

Figure  1-13  shows  an  example  of  what  we  see  with  our  monitor¬ 
ing  camera,  and  Fig.  1-14  shows  a  time  trace  of  the  laser  output  taken 
with  photodiode  No.  1  and  displayed  on  the  Tektronix  519  oscilloscope. 
Figure  I- 14(a)  shows  the  smooth  temporal  shape  seen  over  95%  of  the 
time,  while  Fig.  I- 14(b)  is  an  extreme  example  of  multimode  oscillation. 
In  this  case  two  modes  are  oscillating  with  a  frequency  separation  of 
750  MHz.  This  frequency  corresponds  to  the  spacing  between  the  end 
of  the  ruby  and  the  resonant  reflector.  The  overlapping  spectral  ranges 
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HRL  265-7 


Fig.  1-13. 

Typical  Photo  Monitoring  Fabry- 
Perot  Interferogram,  Near-Field 
and  Far  Field  Beam  Patterns.  The 
Difference  In  Optical  Density 
Between  the  Two  Halves  of  the 
Picture  Is  N.D.  0.6.  The  Free 
Spectral  Range  of  the  Interferom¬ 
eter  Is  1.6  cm" ■ . 


HRL  265-5 


Fig.  1-14.  Tektronix  519  Oscilloscope  Traces  of 
User  Output  with  20  nsec/DI vision 
Sweep  Rate,  (a)  Smooth  Pulse  Ob¬ 
served  ^95%  of  the  Time,  (b)  Modu¬ 
lated  Pulse  (^750  MHz) . 
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of  the  Fabry-Perot  interferometer  and  the  519  oscilloscope  insure  that 
the  oscillation  of  anything  other  than  a  single  longitudinal  mode  will 
be  detected. 

The  characteristics  of  the  laser  are  summarized  in  Table  1-3. 

TABLE  1-3 


Characteristics  of  Ruby  Laser 


Energy  Output4 

12  to  15  mJ  (±15%) 

Pulse  Length 

^20  nsec 

Peak  Power 

0.  6  -  0.  8  MW 

Beam  Radius  (1/e  for  E 
field) 

1  mm  (±10%) 

Beam  Divergence 

0.  6  ±  0.  1  mrad  (full-angle) 

Calculated  Beam  Divergence 

0.44  mrad  (full-angle) 

£ 

Our  energy  monitor  was  checked  against  three  different 
"calibrated"  detectors;  the  error  quoted  (±15%)  reflects 
the  disagreement  of  the  standards. 

T627 


A  typical  set  of  experiments  involves  setting  up  the  sample  relative 
to  the  focusing  lens  for  either  surface  or  bulk  damage.  The  location 
of  the  beam  in  the  sample  is  determined  by  the  use  of  a  He-Ne  laser 
which  is  colinear  with  the  ruby  laser.  The  sample  is  mounted  in  such 
a  way  that  it  can  be  translated  either  vertically  or  horizontally  in  a 
known  way.  After  damage  is  observed  in  a  given  location  the  sample 
is  translated  usually  0.  5  mm  to  a  new  site.  During  the  experiments 
the  sample  is  examii  ad  for  damage  between  shots  using  a  traveling 
microscope  mounted  on  the  optical  bench.  In  this  way  the  sample  i 
not  disturbed  in  relation  to  the  laser  beam,  and  a  given  spot  can  be 

# 

illuminated  many  times  if  desired.  The  amplifier  lamp  is  always 
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pumped  with  a  constant  voltage  monitored  with  a  digital  voltmeter  and 
the  energy  incident  on  the  sample  is  varied  by  rotating  the  first  of  a 
pair  of  polarizing  prisms.  For  each  shot  of  the  laser,  we  monitor  the 
temporal  profile,  the  total  energy  incident  on  the  sample,  near-field 
and  far-field  profiles,  and  Fabry-Perot  spectrum.  Hence,  any  devi¬ 
ation  from  the  desired  conditions  is  noted  and  the  data  handled 
accordingly. 
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SECTION  II 


THEORETICAL  STUDIES: 
OPTICAL  DAMAGE 


The  results  of  our  theoretical  studies  of  optical  damage  during 
the  period  of  this  project  are  summarized  in  this  section.  In  line  with 
the  experimental  study  breakdown  of  Section  I,  this  review  will  be 
divided  into  a  discussion  of  bulk  damage  and  surface  damage. 

A.  BULK  DAMAGE 

1.  Breakdown  Mechanisms 


At  the  time  this  project  began,  the  existing  theories  of  optical 
breakdown  in  ideal  crystals  all  hypothesized  the  creation  by  the  high 
optical  electrical  fields  of  enough  hot  electrons  (of  ionizing  energy)  to 
create  an  ionizing  avalanche.  (For  a  review  of  these  theories,  see 
Appendix  A,  reproduced  here  from  the  second  semiannual  report. )  In 
the  first  semiannual  report  of  this  project,  it  was  shown  (using  the 
same  Frolich  model  for  a  conduction  electron  in  a  polar  lattice  from 
which  other  theories  had  started)  that  the  average  kinetic  energy  of 
such  an  electron  would  not  be  significantly  above  kT  (the  lattice  tem¬ 
perature  times  Boltzmann's  constant)  in  the  optical  fields  at  which  break¬ 
down  was  believed  to  occur.  However,  as  was  stated  there:  "A  pos¬ 
sible  complication  that  we  will  study  further  is  that  a  very  few  electrons 
might  be  created  in  the  high  energy  wings  of  the  electron  distribution.  " 
In  our  fourth  semiannual  report,  the  equations  were  developed 
which  determine  the  electron  energy  distribution  and,  after  several 
simplifying  assumptions,  it  proved  possible  to  cast  them  in  a  form 
that  might  be  solved  on  a  computer.  However,  we  have  not  expended 
the  considerable  effort  that  would  be  required  to  do  this;  efforts  to  see 
by  analytical  methods  whether  a  high  energy  tail  would  be  predicted  in 
the  electron  distribution  have  not  proved  fruitful. 
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Recently,  speculation  has  been  rekindled  by  Bass  and  Barrett 
about  the  plausibility  of  the  existence  of  such  a  high  energy  tail  in  con¬ 
nection  with  their  experimental  results  on  surface  damage.^  They  call 
this  tail  "lucky  electrons"  as  they  have  undergone  just  the  right  kinds 
of  collisions  at  the  right  phase  in  the  optical  field  to  be  continuously 
accelerated  by  it.  However,  in  order  to  appreciate  the  kind  of  theoreti¬ 
cal  effort  required  to  make  even  a  semiquantitative  estimate  of  the 
electron  velocity  distribution,  we  refer  the  reader  to  Appendix  B  of  this 
report  which  reviews  our  development  of  truncated,  but  still  compli¬ 
cated,  equations  for  the  electron  distribution  from  the  FrBlich  model 
that  is  universally  used  in  these  problems. 

In  conclusion,  we  are  certain  that  potentially  damaging  amounts 
of  heat  (—10^  J/cm"'*)  will  be  transmitted  from  the  optical  beam  to  the 
lattice  by  cold  electrons  at  optical  intensities  around  10*^  to  10**  W/cm^ 
if  the  light  has  created  at  least  10^  photoelectrons/cm\  This 
was  demonstrated  in  Semiannual  Report  2  and  reproduced  in  Appendix  A 
here  (and  reported  at  the  Second  Annual  Symposium  on  Laser  Damage, 
Boulder,  Colorado,  June  1970).  Whethe*-  other  damaging  processes 
occur  first  is  not  clear.  However,  it  is  now  known  that  sapphire 
resists  bulk  damage  until  self-focusing  produces  high  enough  intensities 
to  produce  any  of  a  number  of  nonlinear  effects. 

Semiannual  Report  4  documented  the  first  conclusive  experi¬ 
mental  evidence  that  self-focusing  had  actually  been  occurring  when¬ 
ever  bulk  damage  was  produced  in  ruby  or  sapphire.  These  results 
are  reviewed  in  Section  I-B  of  this  report.  Whenever  self-focusing 
sets  the  threshold  for  damage,  it  is  less  interesting  to  know  the  actual 
mechanism  of  rupture  than  the  mechanism  and  threshold  for  self- 
focusing.  Consequently,  after  the  first  demonstration  by  this  project 
that  bulk  optical  damage  in  sapphire  and  ruby  always  follows  self- 
focusing  for  pulses  in  the  10  to  20  nsec  region,  we  have  sought  and 
found  correlations  between  the  tendency  of  a  material  to  self-focus  and 
its  other  material  properties.  In  this  report  we  hypothesize  and  give 
arguments  that  the  tendency  of  crystals  to  self-focus  increases  with 
their  ordinary  Raman  scattering  strength,  provided  the  pulses  are 
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short  enough  to  freeze  out  electrostriction.  This  rule  is  in  contra¬ 
distinction  to  an  anharmonic  oscillator  model,  which  scales  the 
nonlinear  refractive  indices  to  the  linear,  and  which  in  some  instances 
predicts  tendencies  for  self-focusing  opposite  to  those  predicted  by  the 
Raman  effect.  We  have  also  considered  the  use  of  optical  beams  of 
elliptical  cross-section  in  order  to  prevent  e.elf-focusing.  The  experi¬ 
mental  results  obtained  in  this  project  which  show  the  reduced  self- 
focusing  tendencies  of  elliptical  beams  and  the  reason  for  this  effect 
are  discussed  in  Section  I-B.  However  the  detailed  behavior  of  ellipti¬ 
cal  beam  propagation  at  high  optical  powers  is  not  understood  complete¬ 
ly,  and  it  is  possible  that  objectionable  beam  distortions  may  occur 
even  if  self- focusing  is  prevented.  Some  hope  for  reducing  optical 
damage  and  distortion  '  .resholds  in  bulk  lies  in  finding  and  employing 
crystals  with  a  smaller  nonlinear  refractive  index;  our  findings  on  this 
are  reviewed  here.  We  have  also  continued  the  efforts  begun  in  previous 
years  to  understand  the  actual  mechanisms  of  crystal  breakdown  and 
destruction.  The  most  important  of  these  mechanisms  are  linear  and 
nonlinear  heating  and  electron  avalanche  ionization.  Because  surface 
damage  seems  to  be  initiated  by  these  mechanisms  rather  than  by  self- 
focusing,  these  latter  studies  are  probably  more  pertinent  to  the  sur¬ 
face  damage  problem  and  therefore  will  be  discussed  in  Section  II- B 
on  surface  breakdown. 

2.  Mechanisms  Governing  Self -Focusing 

From  our  own  self-focusing  measurements  in  ruby  and  sapphire, 
it  is  now  evident  that  sapphire  has  a  nonlinear  refractive  index  n^  (in 
the  absence  of  electrostriction)  that  is  an  order  of  magnitude  smaller 
than  that  of  laser  glasses.  It  also  appears  that  this  index  is  significantly 
smaller  still  in  YAG  laser  crystals.  One  cannot  understand  this  trend 
by  the  philosophy  of  the  well-known  "Miller's  rule,  "  which  predicts 
that  the  material  with  the  higher  linear  refractive  index  will  have  the 
higher  nonlinear  refractive  index  n^  insofar  as  the  latter  arises  from  the 
nonlinear  polarizability  of  the  electronic  orbits.  The  nonlinear  index 
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for  self-focusing  of  linearly  polarized  light  is  defined  in  terms  of  the 
power  dependent  change  6n  in  the  refractive  index  by 

6n  =  n2  E2  ,  (II- 1) 

where  E  is  the  rms  electric  field.  A  rudimentary  model  of  anhar- 
monically  bound  electrons  predicts  that  n^  scales  like 

n2  «  (n2  -  1  )3  .  (II- 2 ) 

The  linear  indices  of  laser  glass,  ruby,  and  YAG  are  about  1.  53,  1. 76, 
and  1.82,  respectively  from  which  eq.  (II-2)  would  predict  n^  (glass): 
n^  (sapphire):  n^  (YAG)  —2.4:  9:  12.  It  does  not  matter  that  (II- 2 '  is 
crude;  any  anharmonic  n.jde1  would  predict  that  n^  increases  with  n, 
the  opposite  of  what  is  observed.  This  suggests  that  the  fast  self- 
focusing  indices  are  not  entirely  due  to  nonlinear  electronic  distortions, 
but  also  to  nuclear  nonlinearities,  i.  e.  ,  the  field-induced  alterations 
in  local  nuclear  arrargements  which  give  a  lowering  of  the  crystal- 
field  energy  via  the  linear  electronic  response  about  the  altered  nuclear 
positions.  This  nuclear  contribution  to  n^  is  much  slower  responding 
in  time  to  field  amplitude  changes  than  are  electronic  nonlinearities, 
but  much  faster  responding  than  elcctrostrictive  changes  that  arise 
from  a  macroscopic  density  change.  The  nuclear  response  time  is  of 
the  order  of  an  inverse  nuclear  vibration  frequency,  typically  10*^ 2  - 
10  sec.  Therefore,  even  picosecond  pulses  would  experience  self- 
focusing  from  this  nuclear  mechanism;  it  is  presently  difficult  to  dis¬ 
tinguish  nuclear  from  electronic  nonlincarities  by  time  dependence. 

If  nuclear  motions  arc  mainly  responsible  for  n^  in  laser 
crystals,  then  one  might  roughly  and  arbitrarily  guess  how  n^  scales 
by  replacing  n2,  the  electronic  dielectric  constant,  by  (q  -  n^,  the 
nuclear  contribution  to  the  dc  dielectric  constant  in  ( II - 2 ) .  The  values 
of  *o  for  laser  glass,  ruby,  and  YAG  arc  about  6,  10,  and  12, 
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respectively;  the  relation  uI-2)  thus  revised  would  predict  n?  (glass)  : 

'(sapphire)  :  n,  (YAG)  —  1  :  5  :  10,  in  even  worse  disagreement  with 

2  2  1  'i  ,  ' 

experiment.  As  an  alternative  hypothesis,  with  much  more  physical 
basis,  we  propose  here1  that  the  size  of  n ^  is  proportional  to  the  ten¬ 
dency  of  the  crystal  to  scatter  (low  intehsity)  light,  at  least  insofar  as 
electronic  nonlinearities  can  be  neglected.  iWe  have  demonstrated  that 
such  a  relation  is  true  Within  the  linear  dipole  approximation  (LDA)  for 

i  .  7  . 

iaotropic  materials  (glasses  and  fluids),  iJe.  , 

•  *1 

i  1  •  '  ' 

'  S  i 

n,  (nuclear)  =  ,  (II-3) 

i  n  *  i  , 

!  ,  >  •  I  , 

where  S  is.  the  total  scattering  cross-section  per  unit  solid  angle  per 

1  ‘i  I 

unit  volume  of  material  to  scatter  plane  polarized  light  into  parallel 

'  O  1  '  ! 

polarized  light  at  90  *  *  •=  w  0/c>  *  is  Boltzmann's  constant,  and  T  the 
absolute  temperature.  If  Brillouin  scattering  is  included  in  S,  then  n^ 
includes  the  contribution  of  electrostriction  as  well  as  the  nuclear 
redistribution;  if  the  Brillouin  scattering  contribution  to  S  is  omitted,  * 

then  (U-3)  gives  the  fast  nuclear  contribution  only.  We  have  recently  , 

shown  that  (II-3)  is  also,true  under  more  general  assumptions  than  the 
LDA. 

I  ‘  *  , 

Unfortunately,  a  comparison  of  the  inelastic  light  scattering 

strengths  S  for  laser  glasses,  sapphire,  and  YAG  are  not  yet  available. 

However,  it  appears  that  most  glasses  have  a  broad  Raman  band  from 

50  to  500  cm  *  shift  that  is  about  as  intense  as  that  of  liquids  like 

8 

benzene  which  exhibit  large  Rayleigh  wing  scattering.  It  is  alsp 
known  that  Raman  scattering  intensities  are  not  correlated  with  the 

•  I 

refractive  indices  or  dielectric  constants  of  materials;  our  hypothesis 

i  1  i  i 

is  not  unreasonable  on  physical  grounds. 

If  (II-3)  is  valid  for  predicting  self-focusing  in  crystals,  then 
crystals  like  the  alkali  halides,  which  have  no'allowed  first-order 

l 

Raman  scattering,  and  hcncc  a  small  scattering  cross-section  S, 
should  have  uncommonly  high  thresholds  for  self-focusing  of  short 
pulses.  Measurements  on  alkali  halides  to  test  this  hypothesis 

I  ' 

.  :  I  : 

7  3 

i  1 


are  currently  underway  at  other  laboratories.  As  part  of  this  project, 
we  have  estimated  n^  for  the  alkali  halide  crystals,  using  eq.  (15a)  of 
Ref.  7  for  S  in  (II- 3 )  above  and  the  known  ionic  polarizabilities,  masses, 
and  restrahl  frequencies.  In  this  model  the  change  6n^  in  refractive 
index  results  essentially  from  the  field-induced  change  in  the  local 
field  factor  at  an  ion  due  to  the  slight  statistical  anisotropy  induced  in 
the  neighborhood  of  the  ion  by  the  lattice -field  interactions.  The 
results  of  these  calculations  for  some  typical  alkali-halide  crystals 
are  given  in  Table  II - 1  along  with  a  recently  measured  estimate  of  n^ 
for  fused  quartz  for  comparison.  Here  the  nuclear  contributions  to  n^ 
are  so  small  that  one  suspects  that  the  electronic  contribution  must 
dominate. 


TABLE  II- 1 


Calculated  Nuclear  Contribution  to  the  Fast  Nonlinear  Index  n^  for 

—  ^ 

Optical  Self-Focusing  of  Certain  Alkali  Halide  Crystals. 


Material 

n^  (esu)  xlO17 

Rbl 

20 

KI 

10 

KCI 

6 

LiF 

0.  4 

NaCI 

0.  3 

Fused  SiO^ 

10,  000b 

No  electrostriction  or  electronic  hypcrpolarizability  effects  arc 

included.  Total  fast  index  for  fused  quartz  is  given  for 

comparison. 

bF rom  Ref.  9. 

T637 
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Whether  the  nuclear  contribution  dominates  n^  or  not  in  any 

given  material  can  be  determined  directly  by  comparing  the  Kerr  effect 

9 

and  the  intensity  dependent  change  in  polarization  in  the  material. 
Attempts  to  do  this  in  glass,  sapphire,  and  YAG  are  currently  in  pro¬ 
gress  elsewhere.  We  feel  that  with  the  knowledge  gained  in  such 
studies,  there  is  a  real  possibility  of  finding  materials  for  which  bulk 
damage  thresholds  for  short  laser  pulses  are  much  higher,  and  non¬ 
linear  optical  distortions  of  the  beam  are  much  smaller,  than  in  exist¬ 
ing  laser  glass  and  sapphire. 

B.  SURFACE  DAMAGE 

Damage  by  an  unfocused  optical  pulse  in  sapphire,  as  in  most 
materials,  appears  first  at  its  surfaces  and  at  about  the  same  flux  local 
optical  field  densities  for  the  entrance  as  for  the  exit  surfaces.  Our 
investigations  described  in  Section  I-C  show  that,  for  a  given  pulse 
length,  the  damage  threshold  of  sapphire  is  reached  when  the  peak  power 
per  unit  area  exceeds  a  certain  value  at  the  surface,  and  the  threshold 
does  not  depend  on  the  characteristics  of  the  beam  i.'side  the  crystal. 
Mainly  from  these  facts,  we  have  narrowed  the  possible  mechanisms 
for  the  surface  damage  to  three  types.  In  order  of  descending  likeli¬ 
hood,  these  are  (1)  heating  by  normal  linear  absorption  in  a  surface 
layer,  (2)  heating  from  nonlinear  absorption  by  photoexcited  electrons 
in  a  surface  layer,  and  (3)  electron  avalanche  ionization  of  the  surface 
material.  We  discuss  what  each  of  these  mechanisms  predicts  and 
various  experimental  ways  of  verifying  or  discounting  them  below. 

1 .  Heating  from  Linear  Absorption 

Suppose  a  thin  uniform  surface  layer  of  thickness  £  and  absorp¬ 
tion  coefficient  a  is  irradiated  by  an  optical  pulse  of  U  J/cm^  energy 
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density.  If  the  pulse  length  t  is  less  than  i  *  (the  velocity  of  sound*, 
then  the  temperature  rise  AT  in  the  layer  is 

AT  —  q  rr—  (II-4) 

where  Cy  is  the  average  specific  heat  (J  cm'^  °K  * )  at  constant  volume 
in  the  given  temperature  range.  For  sapphire  at  room  temperature, 

C  ~  3  and  a  potentially  damaging  surface  temperature  rise  of  10^  °K 

v  .  i 

would  result  at  U  —  30  J  cm  (a  typical  observed  threshold  value  for 
-  8 

t~-10  seel  if  the  layer  absorption  coefficient  satisfied 

a  >  1  00  cm  * .  ( II -  5 ) 

For  a  layer  thickness  l  —  10  ^  cm,  a  fractional  linear  absorption  of 

_4 

only  10  at  a  surface  could  therefore  be  strongly  suspected  to  give 

rise  to  surface  damage  at  typically  observed  thresholds. 

Layers  of  such  fractional  absorptivities  F  —  10  *  to  10  ^  may 

well  exist  on  even  the  most  carefully  prepared  surfaces.  Suppose,  for 

example,  that  a  certain  surface  density  N/cm^  of  absorbing  cento-  '■ 

(e.  g.  ,  dangling  bonds,  oxygen  impurities,  etc.  )  having  optical  absorp- 

- 1  8  2 

tion  cross-sections —1 0  cm  were  present  in  a  surface  layer.  Then 
only  10  ^  to  10  ^  such  impurities  per  cm^  would  give  rise  to  dangerous 
absorption,  a  not  unreasonable  number. 

This  mechanism  can  be  distinguished  from  the  other  two  leading 
candidates  by  the  dependence  of  the  threshold  power  density  on  the 
pulse  time  r  for  a  carefully  controlled  single  mode  beam.  This  study 
is  difficult  and  has  not  yet  been  done.  Therefore,  we  have  considered 
various  methods  for  detecting  and  measuring  the  a  and  t  of  a  surface 
layer  nondcstructivcly.  The  most  promising  methods  seem  to  be  the 
following. 
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a.  Calorimetry 

For  a  thin  sample  of  material  (<1  mm),  more  of  the 
optical  power  would  be  absorbed  in  the  surface  layer  than  in  bulk,  if 
surface  absorption  is  important  in  surface  damage.  The  surface  absorp¬ 
tion  in  this  case  could  be  measured  by  monitoring  the  temperature  rise 

(after  thermal  equilibration)  produced  in  such  a  thin  sample  by  a  laser 

-4 

pulse.  For  a  fractional  surface  absorption  F  —  10  ,  a  temperature 

rise  of  1/300°C  would  be  produced  by  each  J/cm^  incident  over  a 
sapphire  slab  0.  1  mm  thick.  Thermal  decay  times  could  be  made  such 
that,  for  example,  a  1  W  cw  Nd:YAG  laser  might  easily  deposit 
1000  J/cm^  during  a  thermal  decay  time  and  produce  a  measureable 
temperature  rise  that  would  then  give  the  fractional  absorption  F  =  a l 
for  the  surface  layers. 

b.  Measurement  of  a  Dielectric  Reflection  and/or 

Transmission  Anomaly 

It  is  a  straightforward  matter  to  derive  the  fractions  of 
power  reflected  and  transmitted  T^  from  a  plane  wave  incident, 
through  a  lossless  medium  (of  dielectric  constant  <  =  n^),  upon  a  uni¬ 
form  absorbing  layer  of  complex  dielectric  constant  and  thickness 
g,  that  is  followed  by  a  lossless  semi-infinite  crystal  medium  having  a 
real  dielectric  constant  <c>  The  results  show,  for  example,  that  for 
normal  incidence  on  a  thin  layer,  the  power  reflectivity  Rg  will  deviate 
from  the  reflectivity  Ro  for  the  case  of  no  absorbing  layer  by  a  fraction 


_  ( Rg  -  Rq> 


-(£)(■•?)" 


Im  «£  +  0(0,  (II-6) 


where  wi*  the  optical  (angular)  frequency.  The  fraction  Fg  of  the  inci¬ 
dent  beam  absorbed  in  the  layer  is  (for  Dg  <<  1) 


i-.cA 

T-rrr< 


( II- 7 ) 
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and  is  seen  to  be  of  order  1/2  D|  for  sapphire.  It  is  difficult,  though 
possible,  to  see  a  deviation  in  reflectivity  from  nominal  of  the  order 
of  a  percent.  Therefore,  reflectivity  measurements  are  a  feasible  way 
to  measure  layer  absorptions  of  the  order  of  a  percent  or  greater. 

(The  fractional  change  in  transmission  is  generally  much  smaller  and 
more  difficult  to  interpret  because  it  must  involve  two  surfaces.  ) 


Reflection  at  Nominal  Index  Matching 


One  might  vary  the  index  of  the  incident  medium  (by 

placing  the  crystal  in  a  liquid  of  variable  composition)  in  contact  with 

the  surface  of  the  crystal.  As  the  nominal  reflection  passed  through 

zero  ( t  =  t  ),  the  actual  residual  reflection  due  to  the  layer  would  be 
c 

(at  e  =  0) 


(II-8) 


which  is  roughly  the  square  of  the  fraction  absorbed.  Though  small, 
this  reflection  might  be  measurable  from  a  layer  absorbing  as  little  as 
10  ^  to  10  *  (depending  on  complications  from  birefringence)  because 
there  is  no  normal  background  reflection. 


Residual  Reflection  at  Brewster's  Angle 


Rather  than  vary  the  incident  medium  index,  one  could 
observe  and  plot  the  reflection  versus  angle  of  incidence  and  polariza- 
tion  near  Brewster's  angle  to  discern  the  small  residual  reflection, 
which  again  would  be  of  the  order  the  square  of  the  fraction  absorbed 
in  the  beam.  Therefore,  fractional  layer  absorptions  less  than  10"^ 
could  be  observed  by  this  method. 

Reflection  measurements,  unlike  calorimetric  measurements, 
offer  the  possibility  of  determining  the  layer  thickness  as  well  as 
absorption.  Comparing  reflectivities  at  enough  different  angles  and 
polarizations  would  determine  each  of  Re  <f  ,  Im  <f  ,  and  l.  The  fact 
that  an  absorbing  layer  might  not  be  uniform  with  depth,  or  laterally, 
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would  complicate  the  interpretation,  though  reasonable  estimates  of 
an  effective  F  and  £  could  still  be  made  for  a  given  angle  of  incidence. 


e.  Measurement  of  Reflectivity  Change  Due  to  Predamage 
Heating  of  Surface  Layer 

If  a  layer  is  being  damaged  by  linear-absorption  heating, 
then  a  very  short  lived  but  large  temperature  rise  must  take  place  in 
the  surface  layer  at  pulse-energy  levels  just  below  that  required  for 
damage.  This  large  temperature  rise  should  cause  a  large  change  in 
the  complex  optical  index,  in  and  near  the  layer,  that  should  produce  a 
marked  change  in  reflectivity,  calculable  for  example  at  normal  inci¬ 
dence  by  (II-6).  Perhaps  easier  to  observe  would  be  the  heating  by  a 
polarized  pulse,  incident  upon  the  surface  at  Brewster's  angle.  The 
light  pulse  itself  should  show  an  increase  in  reflection  and  scattering 
during  its  latter  phase,  as  Brewster's  angle  will  have  changed  with  the 
heat-altered  dielectric  constant.  This  increase  could  be  correlated 
with  the  layer  absorption,  and  it  would  be  at  least  as  sensitive  an  indicator 
of  a  surface  layer  as  the  previous  detection  methods. 

f.  Interferometric  Measurement  of  Heating:  Induced  Change 
IrT  Sur  face  Refractive  Index 

The  increase  in  phase  of  a  wave  after  it  has  traveled  a 

z 

distance  z  into  a  medium  is  (in  the  WKB  approximation)  I  u>n'(z')dz'/c, 

o 

where  n'  is  the  real  part  of  the  refractive  index.  The  change  6n'  in  n' 
could  well  be  of  order  unity  in  a  layer  that  is  heated  to  just  below 
damage.  An  interferometric  measurement  of  this  phase  shift  (by 
Michelson,  Mach-Zender,  holographic,  or  other  such  interferometric 
technique)  would  probably  allow  detection  of  a  hundredth  of  a  wavelength 
fringe  shift,  which  would  occur  if  the  absorbing  layer  were  roughly  a 
hundredth  of  a  wavelength  thick  or  of  the  order  of  60  A. 

g.  Measurement  of  Layer  Heating  by  Polarization  Change 
of  Reflected  or  Transmitted  Light 

If  a  beam  of  arbitrary  elliptical  polarization  is  incident 
at  an  arbitrary  angle  on  a  surface,  upon  which  is  a  layer  whose  index 
has  been  changed  by  optical  heating,  then  the  reflected  and  transmitted 
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beams  will  experience  a  change  in  their  polarizations.  Polarization 
analyzers,  which  detect  only  the  polarization  orthogonal  to  that  reflected 
from  the  preheated  surface,  would  sense  the  changed  index  of  refraction 
against  a  very  loW  background.  Straightforward  calculation  would 
reveal  the  incident  polarization  and  angle  that  would  give  the  largest 
effect  for  a  given  material.  The  polarization  measurement  would  be 
sensitive  to  heat-induced  changes  in  the  birefringence  of  the  material  as 
well  as  to  the  changes  in  absolute  index.  The  sensitivity  of  this  method 
would  depend  critically  on  the  degree  to  which  the  polarization  analyzer 
can  be  made  to  reject  the  nominal  reflected  polarization. 

h.  Measurement  of  Optical  Heating  of  Layer  Through 

Change  in  Surface  Resistance 

The  electrical  resistivity  of  any  small  volume  element  at 

or  near  the  surface  of  an  insulator  will  decrease  roughly  as 

A  exp  (T  / T)  with  increasing  absolute  temperature  T.  (For  bulk 
°  4  7  5 

sapphire  Tq~  5  x  10  K  and  A  —  1 0  n  -  cm).  If  any  observable 

resistance  change  could  be  induced  by  a  beam  falling  between  two 
closely  spaced  electrodes  on  a  surface,  one  could  correlate  its  initial 
value  with  (the  absorbed  fraction)  ♦  (layer  thickness  affected).  Uncer¬ 
tainties  in  space  charge,  contact  resistance,  and  surface  conductivity 
dependence  on  temperature  would  weaken  the  correlation.  The  time 
dependence  of  the  resistance  after  a  short  optical  heating  pulse  would 
shed  much  light  on  these  uncertainties  and  on  whether  the  nonlinear 
absorption  of  the  type  discussed  next  was  present. 

L.  Heating  from  Nonlinear  Absorption 

The  mechanism  we  have  proposed  in  this  project  as  responsible 
for  bulk  damage  (sec  Appendix  A  and  Ref.  10)  is  also  likely  to  be  opera¬ 
tive  at  surfaces.  This  mechanism  is  the  nonlinear  absorption  of  the 
optical  beam  by  electrons  excited  to  the  conduction  band  by  the  beam 
itself  through  various  single  or  multistep  excitation  processes.  (We 
refer  to  this  mechanism  as  APE  —  absorption  from  photo  electrons.  ) 

This  absorption  is  expected  to  increase  rapidly  with  beam  power 
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density  P  by  a  relation  that  depends  on  the  nature  of  the  photoexcitation 
process.  It  must  increase  at  least  as  rapidly  as  P;  we  have  argued  that 
it  ought  to  increase  according  to  an  avalanche  type  equation,  once  the 
photoelectron  density  is  high  enough.  *  ® 

Shallow  electron  states  may  abound  near  material  surfaces,  and 
electron-lattice  scattering  is  greater  near  a  surface,  both  effects  tend¬ 
ing  to  make  the  absorption  coefficient  a  from  this  process  much  greater 
near  a  surface  than  in  bulk.  We  have  calculated  that  for  a  free  carrier 
density  around  10^  cm-'*,  a  would  be  at  least  10  *  cm"*  in  a  perfect 
sapphire  crystal;  it  is  therefore  conceivably  orders  of  magnitude 
greater  at  an  irradiated  surface. 

To  establish  the  existence  of  this  mechanism  independently  of 
damage,  one  could  measure  the  surface  photoconductivity  of  a  sample 
irradiated  by  a  laser  pulse  so  short  that  no  space  or  electrode  charges 
would  have  time  to  build  up  during  the  conductivity  pulse.  (Ordinary 
giant  pulses  would  perform  this  function  well.  )  Recombination  times 
are  expected  to  be  much  less  than  a  nanosecond;  therefore,  this  con¬ 
ductivity  pulse  would  end  with  the  optical  pulse,  unlike  the  thermal 
change  in  conductivity  mentioned  in  the  preceding  subsection. 

The  dependence  of  damage  threshold  on  pulse  length  would  dis¬ 
tinguish  this  effect  from  linear  heating.  In  the  nanosecond  pulse  region 
at  least,  APE  would  tend  to  have  a  power  density  threshold,  whereas 
heating  effects  would  cause  an  energy  density  threshold.  The  first 
published,  and  rather  tentative,  evidence  on  pulse  length  variation  of 
ruby  surface  damage  favors  the  latter. 

3.  Electron  Avalanche  Ionization 


From  their  studies  of  optically  induced  surface  damage  of 
various  crystals,  Bass  and  Barrett  have  proposed  that  an  ionizing 
electronavalanche  is  the  responsible  mechanism.^  We  had  previously 
argued  in  this  project  that  their  coupling  to  the  optical  phonons  would 
cool  conduction  electrons  too  fast  to  allow  any  of  them  to  achieve  ener¬ 
gies  much  above  the  ambient  lattice  temperature.  We  attempted  to 
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verify  this  by  examining  the  equations  that  govern  the  fraction  of 
electrons  that  would  have  enough  energy  to  ionize  a  lattice  site  in  bulk, 
but  did  not  pursue  the  extensive  computer  program  that  we  found  would 
have  been  required  to  solve  them  (see  Semiannual  Report  4  of  this  con¬ 
tract,  July  1971  and  Appendix  B  of  this  report);  no  other  investigators 
have  accomplished  more.  The  surface  problem  is  further  complicated 
by  lack  of  reliable  physical  parameters. 

Like  the  APE  mechanism,  this  effect  would  tend  to  have  a  con¬ 
stant  power  density  threshold.  Ionizing  electrons,  if  they  existed, 
might  be  detected  as  high  energy  photoelectrons  ejected  from  the  surface 
into  a  surrounding  vacuum. 
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SECTION  III 


SUMMATION  OF  EXPERIMENTAL  AND  THEORETIC  A  L  STUDIES 


The  work  carried  out  during  this  program  has  led  to  a  number 
of  fruitful  results.  Although  the  emphasis  has  been  placed  on  damage 
in  ruby  and  sapphire  at  6943  A,  many  of  the  phenomena  studied  are 
generally  applicable  to  all  transparent  material  for  visible  and  infrared 
radiation.  The  highlights  of  our  findings  arc  summarized  in  the 
following. 

•  We  have  demonstrated  the  extreme  importance  of  using 
a  well-controlled  and  well-characterized  laser  source 
for  performing  studies  of  this  type. 

•  The  initiating  mechanism  for  bulk  damage  in  inclusion - 
free  sapphire  and  ruby  is  self-focusing.  Good  semi- 
quantitative  agreement  has  been  obtained  between  our 
experimental  results  and  the  theory  for  a  moving  sclf- 
focuscd  spot. 

•  We  have  shown  experimentally  and  theoretically  that  the 
self-focusing  threshold  power  is  a  function  of  beam 
shape,  elliptical  beams  having  a  higher  threshold  than 
circular  beams.  A  solution  to  the  problem  of  self- 
focusing  damage  is  suggested  by  proper  choice  of  beam 
shape. 

•  Spatial  and  temporal  studies  have  shown  that  the  surface 
plasmas  are  a  result  of  rather  than  a  cause  of  surface 
damage.  At  power  levels  well  above  threshold,  however, 
the  surface  damage  morphology  can  be  influenced  by 
interaction  with  the  plasma. 

•  The  threshold  for  surface  damage  of  sapphire  can  be  con¬ 
siderably  enhanced  when  the  surfaces  are  polished  with 
inn  beams.  Marked  changes  in  both  surface  smoothness 
and  crystallinity  take  place  on  ion  polishing.  It  is  not 
yet  known  to  what  extent  these  changes  arc  responsible 
for  the  improvement  in  damage  thresholds. 

•  We  have  developed  a  new  model  for  the  optical  break¬ 
down  of  a  material  based  on  estimates  of  the  relative 
likelihood  of  many  nonlinear  photoelcctron  processes  in 
polar  crystals. 
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•  A  theoretical  framework  has  been  established  for  the 
stjdy  of  the  momentum  distribution  of  conduction  elec¬ 
trons  in  a  polar  lattice  being  driven  by  a  strong  optical 
field.  Although  we  have  not  expended  the  considerable 
computational  effort  to  obtain  numerical  estimates  of 
numbers  of  high  energy  electrons,  we  have  been  able  to 
show  that  many  ad  hoc  hypotheses  that  have  been  put 
forward  are  without  theor'  tical  basis. 

•  We  have  developed  models  of  the  mechanisms  responsible 
for  the  self-focusing  that  triggers  bulk  damage.  We  have 
found  evidence  that  electronic  hyperpolarizability  is  not 
the  dominant  mechanism  in  laser  glass,  sapphire,  and 
YAG  but  probably  is  in  the  alkali  halides.  We  have 
pointed  out  ways  of  establishing  the  mechanisms  by 
standard,  nondamaging,  optical  measurements. 

•  We  have  studied  the  three  types  of  mechanisms  most 
likely  to  be  central  in  surface  damage:  heating  by  linear 
surface  absorption,  heating  by  photoelectron  absorption, 
and  avalanche  ionization.  We  have  isolated  eight  methods 
which  we  show  can  achieve  nondamaging  measurements 
of  surface  absorption  and  distinguish  among  the 
mechanisms. 
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'  SECTION!  IV  i 

i 

! 

RECOMMENDATIONS  TOR  FUTURE  WORK  1 

1 

As  is  typical  in  research  programs  of  this  type,  more  questions 
are  raised  during  the  course  of  the  involvement  than  are, answered.  , 

Because  of  the  scope  of  the  work,  certain  areas  were  pursued  more  , 

actively  than  others.  Also,  termination  of  such  programs  occurs  for 

i  i 

a  number  of  practical  reasons  according  to  a  predetermined  time  inter- 

)  J 

val  and  not  necessarily  when  the  pertinent  questions  have  been  satis- 

*  *  I  i 

factorily  answered.  As  a  result,  we  will  discuss  briefly  a  number  of 
areas  which  deserve  further  attention.  Sorrie  of  these  have  been  pur¬ 
sued  fairly  actively  on  the  present  program  and  should  be  pursued 
further  for  additional  clarification.  Others  are  interesting  areas,  )  i 

t 

which  were  uncovered  during  the  course  of  this  work,  but  for  which  > 

,  I 

there  were  neither  sufficient  capacity  noV  immediate  justification  to 

» 

pursue  in  detail.  :  ,  1 

The  following  pertinent  comments,  not, necessarily  in  order  <?f 

i  ! 

importance,  evolved  from  our  investigations  and  observations  during 

the  course  of  this  program  and  are  based  on  the  results  of  these  studies.  i 

« 

•  Work  of  the  type  carried  out  bn  this  program 
should  be  extended  to  other  materials  and 
wavelengths  (as  is  being  done  already  to  some 

i 

extent)  with  initial  emphasis  op  whether  the 
damage  is  intrinsic  or  extrinsic.  '  I 

•  Self-focusing  experiments  of  the  type  explored  ! 

:  1  i 

in  this  program  with  the  streak  camera  studies 

should  be  continued  in  more  depth  fbr  both  i 

sapphire  and  ruby  bnd  for  other  laser  materials 

'  *  J  I 

as  well.  Pertinent  results  for  pulses  of  different 

duration  and  ,for  different  focusing  conditions  will  (  , 

I  i 

allow  appropriate  modification  of  the  present 

t  •  ,  i 

theory  to  account  for  the  relative  contributions  i 
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of  slow  and  fast  response  mechanisms  to  the 
nonlinear  index  of  refraction. 

•  More  experimental  data  should  be  taken  on  self- 
focusing  with  beams  of  noncircular  spatial  dis¬ 
tributions,  especially  with  elliptical  beams.  A 
comparison  of  the  results  with  theory  in  more 
detail  than  was  possible  in  this  work  is  essential 
in  determining  the  applicability  of  this  phenomenon 
to  avoiding  self-focusing  in  high  power  lasers  and 
amplifiers.  Along  this  same  line,  experimental 
and  theoretical  work  should  be  pursued  in  the  area 
of  optical  resonators  with  elliptical  modes  of 
oscillation. 

•  Intrinsic  breakdown  thresholds  in  the  absence  of 
self-focusing  should  be  studied  for  a  variety  of 
optical  materials.  (A  good  method  for  doing  this 
type  of  experiment  is  described  in  Section  I-B-8.  ) 
This  type  of  measurement  will  yield  information 
about  the  limiting  size  of  self-foci  as  well  as 
further  insight  into  the  nature  of  the  ultimate 
breakdown  mechanisms. 

•  A  number  of  the  differences  between  damage 
properties  of  ruby  compared  with  sapphire 
pointed  out  in  this  report  are  suspected  to 
arise  from  ruby's  absorption  at  6943A  and  pos¬ 
sibly  because  the  laser  and  sample  are  at  differ¬ 
ent  temperatures.  Detailed  exploration  into  the 
nature  of  the  above  differences  should  be  carried 
out  in  order  to  clarify  them. 

•  Further  studies  are  needed  in  the  area  of  surface 
preparation  techniques  such  as  ion  beam  polishing 
or  related  methods  to  surface  damage  resistance, 
not  only  in  the  materials  studied  on  this  program 


but  in  lase:  materials  in  general.  The  reasons 
for  threshold  enhancement  should  be  understood 
so  that  proper  direction  can  be  given  for  future 
work 

•  As  implied  in  the  above  comment,  all  the  avail¬ 
able  surface  diagnostic  techniques  should  be 
brought  to  bear,  where  applicable,  to  studying 
the  effects  of  various  surface  treatments. 

•  The  question  to  be  answered  as  far  as  surface 
damage  is  concerned  is:  How  do  surface  finish, 
crystallinity,  and  material  purity  (e.g.,  presence 
of  surface  states)  affect  the  surface  damage  thresh¬ 
old  of  a  given  material  and  what  are  the  mechanisms? 

•  Some  of  the  techniques  outlined  in  Section  II  for 
measuring  surface  light  absorption  should  be  pur¬ 
sued  in  future  work.  This  type  of  data  will  not 
only  yield  information  as  to  the  purity  of  the  mate¬ 
rial  but  also  will  give  some  clues  to  the  pertinent 
damage  mechanisms. 

•  The  equation  derived  in  Appendix  B  on  the  theory 
of  electron  velocity  distribution  in  a  crystal  irra¬ 
diated  by  a  strong  optical  beam  should  be  solved 
numerically.  Although  considerable  effort  and 
computer  time  will  be  involved  in  accomplishing 
this,  the  results  promise  to  be  quite  pertinent 
for  exploring  the  contribution  of  avalanche  break¬ 
down  to  optical  damage. 
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SECTION  V 


PRESENTATIONS  AND  PUBLICATIONS 


HRL  work  on  this  program  has  resulted  in  a  number  ol 
presentations  and  publications  which  are  listed  here. 

A.  PRESENTATIONS 


•  "Damage  Threshold  Studies  in  Ruby  and  Sapphire" 

C.R.  Giuliano  and  L.  D.  Hess,  presented  at  the 
ASTM-NBS  Symposium  on  Damage  in  Laser 
Materials,  Boulder,  Colorado,  June  1970. 

•  "Role  of  Photo-Electrons  in  Optical  Damage" 

R.W.  Hellwarth,  presented  at  the  ASTM-NBS 
Symposium  on  Damage  in  Laser  Materials, 
Boulder,  Colorado,  June  1970. 

•  "Optical  Damage  in  Polar  Crystals" 

R.W.  Hellwarth,  presented  at  the  Gordon  Con¬ 
ference  on  Nonlinear  Optics,  Issaquah,  Wash¬ 
ington,  August  1970. 

•  "Damage  Threshold  Studies  in  Laser  Crystals" 

C.R.  Giuliano,  L.  D.  Hess,  and  E.S.  Bliss, 
presented  at  the  Sixth  International  Quantum 
Electronics  Conference,  Kyoto,  Japan, 
September  1970. 

•  "Optical  Damage" 

R.W.  Hellwarth,  presented  at  the  Physics 
Colloquium  of  the  Clarendon  Laboratory,  Ox¬ 
ford,  November  1970. 

•  "Time  Evolution  of  Damage  Tracks  in  Sapphire  and 
Ruby" 

C.R.  Giuliano,  presented  at  the  ASTM-NBS 
Symposium  on  Damage  in  Laser  Materials, 
Boulder,  Colorado,  May  1971. 
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•  "Self-Focusing  of  Gaussian-EUiptic  Light  Beams" 

C.R.  Giuliano,  J.H.  Marburger,  and  A.  Yariv, 
presented  at  the  Seventh  International  Quantum 
Electronics  Conference,  Montreal,  Canada, 

May  1972. 

•  "The  Relation  between  Surface  Damage  and  Surface 
Plasmas" 

C.  R.  Giuliano,  presented  at  the  ASTM-NBS 
Symposium  on  Damage  in  Laser  Materials, 
Boulder,  Colorado,  June  1972  (to  be  published 
in  the  Proceedings). 

•  "Ion  Beam  Polishing  as  a  Means  of  Increasing  Surface 
Damage  Thresholds  in  Sapphire" 

C.  R.  Giuliano,  presented  at  the  ASTM-NBS 
Symposium  on  Damage  in  Laser  Materials, 
Boulder,  Colorado,  June  1972  (to  be  published 
in  the  Proceedings). 

•  "Fundamental  Absorption  Mechanisms  in  High-Power 
Laser  Window  Mate  rials" 

R.W.  Hellwarth,  presented  at  the  ASTM-NBS 
Symposium  on  Damage  in  Laser  Materials, 
Boulder,  Colorado,  June  1972  (to  be  published 
in  the  Proceedings.) 

•  "Laser-Induced  Damage  in  Optical  Materials" 

C.  R.  Giuliano,  presented  at  the  Gordon  Confer¬ 
ence  on  Nonlinear  Optical  Effects,  Beaver  Dam, 
Wisconsin,  July  1972. 


PUBLICATIONS 


•  "Damage  Threshold  Studies  in  Ruby  and  Sapphire" 

C. R.  Giuliano  and  L.  D.  Hess,  National  Bureau 
of  Standards  Special  Publication  341,  U.S.GPO, 
Washington,  D.C.  (Dec.  1970)  p.  76  (ed.  by 
A.J.  Glass  and  A.H.  Guenther). 

•  "Role  of  Photo-Electrons  in  Optical  Damage" 

R.W.  Hellwarth,  National  Bureau  of  Standards 
Special  Publication  341,  U.S.GPO,  Washington, 

D. C.  (Dec.  1970)  p.  67  (ed.  by  A.J.  Glass  and 
A.  H.  Geunther) . 


•  "Time  Evolution  of  Damage  Tracks  in  Sapphire  and 
Ruby" 

C. R.  Giuliano,  National  Bureau  of  Standards 

Special  Publication  356,  U.S.GPO,  Washington, 

D.  C.  (Nov.  1971)  p.  44  (ed.  by  A.  J.  Glass  and 

A.  H.  Geunther). 

•  "Observations  of  Moving  Self-Foci  in  Sapphire" 

C.R.  Giuliano  and  J.  H.  Marburger,  Phys.  Rev. 

Letters  27,  905  (1971). 

•  "Laser-Induced  Damage  in  Transparent  Dielectrics: 

Ion  Beam  Polishing  as  a  Means  of  Increasing  Surface 
Damage  Thresholds" 

C.  R.  Giuliano,  Appl.  Phys.  Letters  2_1_,  39(1972). 

•  "Enhancement  of  Self-Focusing  Threshold  in  Sapphire 
with  Elliptical  Beams" 

C.R.  Giuliano,  J.H.  Marburger,  and  A.  Yariv, 

Appl.  Phys.  Letters  21 ,  58  (1972). 

•  "Laser-Induced  Damage  in  Transparent  Dielectrics:  The 
Relationship  between  Surface  Damage  and  Surface 
Plasmas" 

C.R.  Giuliano,  to  be  published  in  IEEE-J. 

Quantum  Electronics  8,  749  (1972). 
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APPENDIX  A 


Role  of  Photo- Electrons  in  Optical  Damage 


1 


R.  W.  Hellwartn 


Department  of  Electrical  Engineering 
University  of  Southern  California 
Los  Angeles,  California  90007 
and 

Hughes  Research  Laboratories 
Malibu,  California  90265 


Fr^vious  workers  have  discussc  1  the  possibility  that  electrons  photo- 
excited  into  the  conduction  band  of  a  p  >lar  crystal  by  a  high-intensity,  short 
optical  pulse  may  gain  enough  energy  to  damage  the  lattice  by  initiating  an 
ionizing  avalanche.  V/e  discuss  here  other  processes  by  which  these 
conduction  electrons  may  damage  the  lattice  even  before  they  gain  enough 
energy  to  ionize  their  surroundings,  specializing  our  calculations  to  the  case 
of  sapphire  and  ruby.  First,  we  show  that  the  energy  that  the  conduction 
electrons  absorb  linearly  from  the  optical  beam  is  deposited  almost 
immediately  in  the  lattice  without  significant  heating  of  the  electrons.  ,*t 
electron  densities  (  —  10*"  cm*3)  and  optical  intensities  (  ~  1 0 •  0  W/cm2) 
likely  to  exist  at  sapphire  damage  thresholds,  this  deposited  energy  is  found 
to  be  of  the  order  of  what  one  might  expect  would  be  required  to  form  a 
rupturing  shock  wave.  We  also  show  that  the  photo-excitation  of  both  bound 
and  unbound  impurity  levels  is  enhanced  by  the  presence  of  conduction 
electrons,  so  markedly  so  for  the  former  that  the  promotion  of  electrons 
into  the  conduction  band  may  be  significantly  "bootstrapped,"  thereby 
increasing  the  optical  absorption.  The  presence  of  conduction  electrons  and 
excited  impurities  in  the  expected  numbi  rs  is  likely  to  alter  the  refractive 
index  significantly  and  affect  thereby  the  focusing  (self-  or  external)  of  the 
beam  in  a  complicated  way.  For  the  si*  .plest  model  these  nonlinear  index 
contributions  would  tend  to  produce  repeated  focal  regions  along  the  beam. 
Implications  of  these  results  for  raising  damage  thresholds  arc  discussed. 


Keywords:  Crystals,  electrons,  optical  damage,  photo- absorption, 
photo-conductivity,  photo-electrons,  polar  crystals, 
ruby,  sapphire,  self- focusing. 


1.  Introduction 


The  physical  processes  responsible  for  the  bulk  damage  caused  in  various  transparent  crystals 
by  short  optical  pulses  (causing  negligible  clcctrostriction)  have  not  yet  all  been  identified.  Here  we 
argue  that  several  processes  not  considered  previously  arc  likely  to  be  important  in  the  optical 
damage  of  inclusion-free  polar  crystals,  especially  sapphire  and  ruby.  For  their  initiation,  these 
processes  would  all  seem  to  require  on  the  order  of  1 0 1 conduction  electrons  per  cc  to  be  present  at 
the  point  of  maximum  optical  intensity,  a  number  widely  suspected  to  be  present  in  sapphire  and  ruby 
at  peak  intensities  ~  10*0  w/cm*  just  below  damage  threshold.  We  will  argue  here  '.hat  such  densities 
of  cold  electrons  can  a)  transfer  damaging  amounts  of  heat  from  the  beam  to  the  lattice;  b)  accelerate 
the  rates  of  photo- excitation  of  bound  and  free  impurity  electron  slates,  thus  significantly  increasing 
the  supply  of  conduction  electrons,  possibly  even  in  "avalanche";  and  c)  significantly  change  the  local 
refractive  index,  possibly  in  a  way  which,  in  conjunction  with  the  normal  nonlinear  index,  could  cause 
repealed  beam  focusing  along  its  axis.  In  none  of  these  processes  do  the  electrons  become  hot  enough 
(  —  8  eV)  to  ionize  the  lattice,  in  contradistinction  to  the  process  suggested  by  previous  workers  ;i-4]2. 


*  Work  supported  in  part  by  the  Joint  Services  Electronics  Program  (U.  S.  Army,  Navy,  and  Air  Force), 
under  Grant  No.  AF- AFOSR-69-  1 622A.  and  in  part  by  the  Advanced  Research  Projects  Agency  through 
Air  Force  Cambridge  Research  Laboratories. 

^Figures  in  brackets  indicate  the  literature  references  at  the  end  of  this  paper. 
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Hence,  we  will  refer  to  them  aa  "cold"  conduction  electrons.  How  the  damage  proccasca  discusaed 
here  depend  on  crystal  temperature,  optical  wavelength  and  pulse  duration  ia  a  complex  function  of  the 
energies  and  wavcfunctions  of  the  states  of  the  impurities  supplying  and  trapping  conduction  electrons, 
the  electrons'  effective  mass(es),  the  longitudinal  optical  (LO)  phonon  structure,  and  the  electron- 
phonon  coupling.  Although  only  the  latter  two  are  fairly  well  known  for  ruby  and  sapphire,  plausible 
hypotheses  about  the  former  lead  one  to  single  out  the  processes  we  propose  as  potentially  important 
to  damage.  It  would  appear  that  reducing  donor  concentrations  or  increasing  their  electron  binding 
energies  might  predut  e  more  damage-resistant  crystals.  However,  if  the  conduction  electron  density 
were  proportional  to  optical  intensity  (which  appears  not  to  he  the  ease  it  such  high  intensities),  it 
would  tend  to  reduce  self-focusing  (or  produce  a  n«.t  defocusing)  and  thereby  lower  damage  thresholds 
Other  conclusions  and  suggested  experimental  checks  of  electron  processes  arc  mentioned  in  the  final 
Section.  In  the  following  Section  2  we  summarize  the  parameters  and  physical  regimes  we  encounter 
ir.  ruby  and  sapphire  and  we  outline  our  rough  criteria  for  when  damage  is  expected.  In  Section  3  we 
define  the  approach  we  take  for  evaluating  the  required  parameters  and  discuss  difficulties  with 
previous  calculations  of  the  conditions  under  which  hot  ionizing  electrons  might  be  produced.  The 
calculations  of  the  enhancement  by  conduction  electrons  of  lattice  healing  and  photo-excitation  follow 
in  Sections  4  and  5.  Before  proceeding,  we  first  summarize  our  assessment  of  often-conflicting 
reports  on  photoconductivity  in  ruby  and  sapphire. 


Hochuli  [5]  has  measured  the  low-frequency  photoconductivity  at  various  wavelengths  and  temper¬ 
atures  of  ruby  and  sapphire  at  optical  intensities  I  —  0.  2  to  5  W/cm^  and  using  applied  voltages  at 
frequencies  zero  Hz,  100  Hz,  and  9-  39  GHz.  He  found  no  essential  differences  in  results  for  the  two 
crystals  The  observed  conductivities  (at  5  W/cm^)  were  4  x  1 0"  1  3,  10-11,  and  4.  2  x  10‘#  (Km)'l, 
respectively,  at  the  above  frequencies.  Since  all  these  frequencies  are  well  below  the  electron  col¬ 
lision  frequency  (which  from  the  theory  of  Section  4  we  estimate  to  be  ~  10*4  rpS),  the  wide  variations 
in  results  indicate  the  magnitude  of  experimental  pitfalls  in  any  attempt  to  measure  such  low  bulk 
conductivities.  Anomalies  at  the  two  lower  frequencies  probably  arise  from  surface  effects,  from  the 
inability  of  electric  contacts  to  inject  charges  into  the  crystal,  and  from  space  charge  buildup.  There¬ 
fore,  the  9.  39  GHz  value,  obtained  in  a  microwave  cavity  with  a  small  electron  drift  excursion 
amplitude  (of  the  order  of  an  Angstrom),  is  probably  the  most  reliable.  This  conductivity  was  linear  in 
the  optical  intensity  within  the  range  of  observations.  Its  wavelength  dependence  suggested  that  the 
electrons  were  supplied  from  donor  levels  between  0.  6  and  1.  2  eV  below  the  conduction  band,  in 
agreement  with  what  one  would  conclude  by  studying  the  normal  bulk  conductivity  variation  with 
temperature  at  higher  temperatures.  The  wavelength  dependence  also  suggested  more  donor  levels 
appearing  around  3  eV.  Although  Hochuli  was  able  to  observe  a  Hall  voltage,  it  did  not  vary  when  the 
light  was  turned  on  or  off.  The  nominal  Hall  mobility  value  p  =  0.  052  cm^/V  sec.  derived  therefrom 
has  nevertheless  been  used  without  question  by  other  workers  whenever  a  value  was  required  in 
calculations.  This  value  corresponds  to  a  collision  time  r  of  4  x  10'*^  sec.  ,  or  an  electron  mean  path 
of  the  order  10“*®  cm,  and  must  be  considered  unphysical,  corroborating  the  difficulty  of  making 
electrodes  on  sapphire.  Using  the  standard  Frohlich  theory  of  a  conduction  electron  in  a  polar  crystal 
adapted  to  the  large  coupling  constant  of  sapphire  [6’,  one  estimates  r  (at  room  temperature)  to  be 
between  lO"*’-^  and  see,  depending  on  where  the  effective  band  mass  lies  between  1/10  and  10 

electron  masses.  Assuming  a  free  electron  mass  for  which  t  =  7  x  10-15  sec  (and  p  =  1  2  cm^/V  sec), 
one  would  infer  from  Hochuli's  microwave  observations  that  he  produced  —  2  x  10^  photo-electrons  per 
cc  with  5  W/cm^  of  broadband  Hg  lamp  exciting  intensity.  We  see  from  this  that  10*6  electrons  per  cc 
would  very  likely  be  excited  near  damage  thresholds  of  10'®  W/cm^  provided  that  10*®  donor  levels  per 
cc  were  occupied,  a  provision  which  is  unknown. 


Except  possibly  for  the  experiments  of  Belikova,  et  al.  [7,8],  other  reported  observations  of 
photoconductivity  in  ruby  appear  to  be  plagued  by  surface  and  contact  effects.  Belikova,  et  al.  , 
observed  the  photoconductivity  to  be  highly  nonlinear  with  the  (6943  A)  optical  intensity  I  for  intensities 
near  10*®  W/cm^  [7].  One  can  use  our  value  of  -  above  with  their  data  to  estimate  very  crudely  (since 
their  beam  geometry  was  unspecified)  that  they  were  observing  well  over  10*6  electrons/cc  at 
lOl®  W/cm‘  (provided  that  the  electron  recombination  time  was  short  compared  with  the  optical  pulse 
length).  Belikova,  et  al .  ,  also  observed  optical  emission  bands  near  2,  2.  7,  and  3  eV  (they  could  not 
see  below  2  eV)  which  arc  suggestive  of  some  impurity  level  transitions.  The  foregoing  is  about  all 
one  can  conclude  about  photo-electrons  from  existing  data  for  ruby  and  sapphire. 


2.  Description  of  Electron  Interactions 


In  studying  the  role  of  photo-excited  electrons  in  crystal  damage,  we  will  constantly  apply  our 
formulae  to  the  ease  of  sapphire  and  ruby.  The  phonon  structure  of  these  two  crystals  are  essentially 
identical,  as  arc  the  low-intensity  photoconductivities.  Although  distinct  differences  in  damage 
behavior  have  been  observed,  the  gross  damage  thresholds  arc  statistically  indistinguishable  for 
sapphire  and  ruby.  Therefore,  most  pertinent  aspects  of  these  two  crystals  can  be  studied  together. 
For  brevity,  w».  shall  refer  hereafter  only  to  sappni.-c,  but  shall  intend  our  remarks  to  apply  also  to 
ruby  unless  stated  otherwise.  We  shall  also  assiunc  throughout  that  the  optical  beam  wavelength  is 
6943  A.  that  of  the  room  temperature  ruby  laser. 
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A  free  electron  gas  of  density  lO1^  cm*  3  has  a  plasma  frequency  <*>„  6  x  1012  rps  that  is  much 

lower  than  ».)  the  optical  beam  frequency  of  interest  y  (2.  71  x  10*5  rps),  o)  the  phonon  frequencies  <J|, 

(  ~  10- 1  .  os),  and  c)  typical  electronic  excitation  frequencies  (  ~  1  O'  6  rps). 

We  shall  call  a  conduction  electron  "cold"  when  its  energy  is  much  less  than  interband  energies 
and  its  wavelength  it  much  longer  than  the  scale  of  potential  variations  of  the  field  lattice.  Other 
electrons  we  will  refer  to  as  "hot."  The  time  for  cold  electrons  in  the  1/40  to  3  volt  range  to  equili- 
briatc  among  each  other  by  collisions  at  a  density  10*6  cm- 3  is  of  the  order  of  10*'®  to  10"'  3  sec  [9], 
much  longer  that,  the  inverse  electron-lattice  collision  time  (~  lf>-'4  sec  in  ruby).  Therefore,  at 
least  before  any  large  ionizing  avalanche  occurs,  the  cold  cond-  tion  electrons  may  be  studied  as 
independent  particles,  not  interacting  among  themselves  but  in  .-racting  only  with  the  fie’ds  and  other 
particles  in  the  crystal. 

When  a  cold  electron  moves  in  a  polar  crystal  such  as  ruby  or  sapphire,  it  is  often  called  a 
"polaron.  "  Its  interactions  with  the  crystal  phonons  arc  known  to  be  well-described  at  room  temper¬ 
ature  and  above  by  the  Hamiltonian  derived  by  Frohlich  [10]  and  used  by  all  previous  workers  in 
calculating  electron  motions  leading  to  damage: 


HF  =  ho  +  VF  (,) 

where 

ho  =  *[£  *  AW2  +  zk  “V  ak  ak  (2a) 

and 

VF  =  £k  Ck  ak  C~'2  +  h  e.  (2b) 


Here  j>  is  the  electron's  momentum,  and  x  the  electron's  position  coordinates,  ak  is  the  annihilation 
operation  for  a  LO  phonon  of  frequency  uj?'.  Optical  wavelengths  are  long  cnouglv^so  that  the  vector 
potential  A(t)  may  be  taken  as  a  function  of  time  only  oscillating  at  a  frequency  y.  We  use  dimension¬ 
less  "polaron"  units  in  which  the  effective  band  mass  of  the  electron  m^  is  unity  as  is  Planck's  constant 
A .  All  energies  arc  measured  in  units  of  an  effective  LO  phonon  energy  aWq  (840  cm"'  or  1200s  K  for 
sapphire  and  ruby).  Frequencies  are  measured  in  units  of  u.  (=  1.6  x  10*  ‘•rps)  so  that  y  =  17  for  a 
6943  A  beam  in  sapphire.  Lengths  arc  measured  in  units  of  (A/m^  <«>0)£  which  is  8.  5  x  10-8  cm  for  ruby 
if  we  take  mb  to  be  the  free  electron  mass  m..  Frohlich  showed  that  in  those  units  the  coupling 
coefficients  Cjj  are  well  approximated  by  V"*2^'4n*'2a  k"‘  where  V  is  the  crystal  volume,  a  is  the 
dimensionless  polaron  coupling  constant  (e"'  -  cj  )  (Ryd'/Au0)$.  Here  Ryd'  is  the  Rydberg  for  an 
electron  with  mass  m0,  and  c  and  e0  arc  the  optical  and  static  dielectric  constants  (3.  1  and  10).  For 
mb  "  mOI  a  ~  3  for  ruby.  In  ruby  an  electron  is  hot  if  k  is  of  order  10  or  more  (in  polaron  units).  Hot 
conduction  electrons  have  neither  the  simple  forms  of  kinetic  energy  or  lattice  interaction  energies 
found  in  cq  (1)  but  obey  complex  equations  in  which  exchange  must  be  accounted  for  and  which  have 
never  satisfactorily  been  approximated  in  usable  form  for  a  dynamic  lattice.  There  would  seem  to  be 
no  way  at  present  to  make  a  reliable  estimate  of  how  strong  an  optical  field  would  be  required  to 
produce  enough  interband  electronic  transitions  (i.  e.  ,  lattice  ionization)  to  cause  damage. 


From  the  discussion  of  the  previous  section,  it  is  evident  that  the  electrons  photo-excited  (below 
damage  thresholds)  in  ruby  and  sapphire  crystals  come  from  donor  impurities  to  which  they  are  bound 
initially  with  much  less  energy  than  the  valence-conduction  band  gap  energy  (  ~  8  eV).  The  effects  of 
these  impurities  on  the  supply  and  motions  of  conduction  electrons,  as  well  as  on  the  excitation  of  the 
donors  by  conduction  electrons,  can  be  studied  by  considering  the  Coulomb  interactions  vt  between  a 
conduction  electron  at  x  and  the  i1*1  particle  of  the  alh  impurity  having  charge  cqSand  located  at  ^9  . 


v 


t 


(3) 


where  the  electronic  charge  c  is  4  in  our  polaron  units,  and  is  screened  by  the  dielectric  constant  es 
appropriate  to  the  frequencies  wf  motions  under  consideration. 

To  deal  with  the  clastic  scattering  of  the  electrons,  one  needs  only  the  matrix  clement  of  eq  (3) 
diagonal  in  the  impurity's  ground  state.  This  gives  an  effective  classical  scattering  potential  seen  by 
the  conduction  electron.  To  study  inelastic  scattering,  we  will  assume  that  the  impurities  may  be 
described  by  an  unperturbed  Hamiltonian  h,  =  T.„  hto  with  "ionized"  or  unbound  as  well  as  bound 
electronic  eigenstates.  Calculations  will  therefore  start  from  a  total  Hamiltonian  II  comprised  of 
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Frohlich's  Hamiltonian  eq  (I)  (which  describes  cold  clcctrona,  phonona,  and  their  interactional  plua 
the  electron  impurity  interactiona  of  eq  (3)  and  with  h(  determining  impurity  atatea: 


H  =  ho  +  ht  *  VF  +  Vt' 

We  will  assume  that  damage  to  the  crystal  occurs  when  the  interactions  vj-  and  vt  combined 
cause  the  conduction  electrons  to  absorb  an  energy  U  per  cubic  centimeter  with  a  short  optical  pulse 
that  exceeds  a  threshold  value  UQ.  That  is,  we  shall  strive  to  determine  at  what  optical  beam  intensity 
we  may  expect 


U  a  U0.  (5) 

Present  evidence  indicates  that  UQ  is  of  order  10^  to  10^  J/cm'  for  ruby  and  sapphire  for  20  nsc 
pulses  in  a  beam  of  order  one  mm  diameter.  A  shock  wave  of  this  energy  surrounding  the  beam  would 
develop  stresses  on  the  order  of  the  static  yield  stress  *3]  Obviously  a  large  part  of  understanding 
optical  damage  in  crystals  entails  developing  an  accurate  theory  of  the  amount  of  energy  deposition 
required  to  produce  rupturing  strains  as  a  function  of  the  amount  and  spatial  and  temporal  distribution 
of  this  deposition.  However,  we  shall  content  ourselves  here  with  the  above  crude  estimate  for  UQ. 


3.  Approaches  to  the  Calculation  of  Damage  Thresholds 


Previous  treatments  of  the  role  of  electrons  in  damage  [1,4]  started  from  Frohlich's  Hamiltonian 
cq  (l)  and  sought  to  calculate  the  optical  intensity  at  which  electrons  v-ould  gain  enough  energy  u  to 
excite  an  avalanche  of  interband  transitions.  The  expected  energy  u  of  an  electron  in  the  optical  field 
was  estimated  from  the  equation 


du 

dt 


'fe 


ef 


(6) 


where  Rfc  is  the  average  rate  at  which  the  optical  field  does  work  on  the  electron  and  Re£  is  the 
average  rate  at  which  the  excited  electron  "radiates"  phonons,  i.  e.  ,  transfers  its  energy  to  the 
lattice.  Previous  treatments  have  all  taken  for  Rfc  the  usual  linear  conductivity  loss  per  electron  with 
reasonable  order  of  magnitude  estimates  from  cq  ( 1 )  for  the  optical  conductivit  (or  "electron 
collision  rate").  However,  they  have  made  approximations  for,  or  statements  about,  Rct  that  are  not 
directly  based  on  a  quantum  average  of  the  appropriate  operator.  Wasserman  [1]  took  for  RC£  the  rate 
at  which  an  electron  having  a  Maxwellian  distribution  would  lose  energy  to  a  cooler  lattice  in  the 
absence  of  the  optical  field  (A  =  0).  He  predicted  breakdown  intensities  orders  of  magnitude  above 
those  actually  observed.  Svercv,  ct  al.  *4],  used  for  RC£  the  approximate  rate  at  which  an  electron 
would  lose  energy  to  the  lattice  starting  from  a  momentum  eigenstate  |p)  for  which  £p^  =  u  also  in  the 
absence  of  the  optical  field.  His  predicted  thresholds  also  exceeded  observed  thresholds  by  over  two 
orders  of  magnitude.  We  feel  it  is  quite  probable  that  a  precise  evaluation  of  the  spectral  energy 
distribution  of  an  electron  obeying  eq  (I)  would  show  that  the  electrons  do  in  fact  remain  too  cool  to 
ionise  impurities  or  the  lattice  at  optical  fields  below  damage  thresholds.  In  the  next  section  we  show 
that,  indeed,  collisions  with  phonons  dominate  an  electron's  motion.  Ve  have  checked  that  it  follows 
rigorously  from  eq  (1)  that,  to  lowest  order  in  the  optical  intensity  and  coupling  parameter  a,  an 
electron's  expected  kinetic  energy  is  simply  the  sum  of  its  coherently  oscillating  energy  plus  3kT/2 
where  T  is  the  ambient  lattice  temperature.  The  coherent  energy  is  very  much  smaller  than  kT  for 
optical  intensities  of  interest.  Therefore,  we  now  proceed  to  examine  whether  cold  electrons  may  not 
mediate  the  deposition  of  damaging  amounts  of  energy  in  a  sapphire  lattice. 

As  we  mentioned  in  the  previous  section,  our  approach  will  be  to  calculate  the  beam  intensity  at 
which  the  energy  U  deposited  per  unit  volume  in  the  lattice  from  the  optical  field  via  the  electrons 
exceeds  a  threshold  value  UQ.  If  the  maximum  optical  beam  energy  per  unit  area  is  S,  then 


U  =  x  S  (7) 

v 

where  xu  is  the  absorption  per  unit  length  of  an  optical  beam  of  frequency  u-  Following  FHIP  [6],  we 
write  the  quantum  expectation  value  of  the  amplitude  of  an  electron's  coordinate  at  frequency  v  as 
Re  Ee*v,/(-» ^+Xu)  when  the  electron  is  in  an  optical  electric  field  Re  Ee*tft.  With  this  definition  we 
have 
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(8) 


k#*  =  i<*>p  /  c)  Im[y2  -  xuj*‘ 

where  X  ia  the  wavelength  of  the  light  in  the  cryatal  +  2n  and  w  is  the  electron  plasma  frequency  (in 
units  of  u>o).  In  cgs  units  =  4n3c2/m|J  with  P  being  the  numbBr  density  of  electrons. 

The  evaluation  of  k  falls  naturally  into  two  parts:  first,  the  determination  of  the  rate  at  which  a 
single  conduction  electron  takes  energy  from  the  field  (and  delivers  it  to  the  lattice)  by  a  calculation  of 
X„.  and  second,  the  estimation  of  the  number  density  o  of  these  electrons.  We  consider  these  two 
questions,  respectively , in  the  following  two  sections. 


4.  Optical  Response  of  Single  Electron 


We  examine  first  an  electron's  optical  response  in  the  case  that  it  is  affected  only  by  electron- 
phonon  interactions  as  described  by  Frohlich's  Hamiltonian  eq  (1)  with  coupling  parameters  appropriate 
to  sapphire.  We  then  find  that,  at  least  at  room  temperatures  and  above,  collisions  with  impurities 
arc  relatively  less  important.  When  the  electron  response  is  linear  in  the  applied  field  for  a  given 
lattice  temperature,  we  may  use  the  determination  of  xv  by  FHIP  who  calculated  the  quantum  expected 
value  of  the  electron  position  to  first  order  in  the  optical  field.  Their  expression  is  exact  for  all 
temperatures  at  small  coupling  (a  <  1)  and  gives  an  accurate  solution  even  when  as  for  ruby.  In 

this  case  it  shows  that,  for  y  =  17,  both  the  real  and  imaginary  parts  of  xv  are  much  less  than  o2  so 
that  we  need  only  calculate  Imxv  to  estimate  the  absorption  constant.  (The  Re  xu  can  of  course  be 
determined  from  im  x„  by  Kramers- Kronig  relations.  ) 


Im  x 


sinh(Pv/2)  v_ 


3  sinh(fi/2) 


V  w 


yr 


cob(ju)  cosu  du _ 

r  2  ‘  2  T  '  7737 2  • 

[a  +  a  -  b  cos(vu) J 


(9) 


Here  v  and  w  are  parameters  to  be  chosen  from  a  variational  principle,  a2  =  p2/4  +  Rp  coth(pv/2); 

R  h  (v2  -  ’v2)  /  (w2  v);  and  b  »  Rp  sinh(pv/2). 

In  th-:  eak  coupling  limit  v  =  w  =  3,  R=b  =  0  and  a  =  p/2.  Then  the  integral  can  be  evaluated 
exactly  in  ’.orms  of  modified  Rcssel  functions.  At  temperatures  much  lower  than  the  Ecbyc  temperature 
(P  large)  this  result  reduces  to 


Im  x  -  2a(v  -  1)^/3  :  a  <  1,  p  >>  1,  v  >  1.  (10) 

For  a  =  3  and  y  =  17,  this  is  8.  The  leading  correction  to  eq  (10)  for  finite  temperatures  multiplies  it 
by  (1  4  p*l),  increasing  it  by  30%  at  room  temperature  (p  =  3.  3). 

When  a  =  3,  the  variational  principle  gives  v  =  3.4  and  w  =  2.  55  in  the  zero  temperature  limili.6]. 
Using  these  parameters  with  p  =  3.  3  gives  R  =  0.  229,  a  -  1. 866  and  b  =  5.  53  x  10-3.  We  have  evalu¬ 
ated  eq  (9)  numerically  for  these  parameters  and  have  obtained  Im  xy  =  18.  We  have  also  evaluated  eq 
(9)  at  frequencies  in  the  neighborhood  of  y  -  17  and  found  that,  as  in  eq  (8),  Im  Xy  varies  slowly  with  v. 
Therefore,  electron  transients  associated  with  the  rise  and  fall  of  even  picosecond  optical  pulses  are 
very  small,  and  it  is  a  good  approximation  in  practice  to  assume  that  at  any  instant  the  electron's 
motion  is  the  same  as  for  a  purely  nnusoidal  field  of  amplitude  appropriate  to  the  intensity  at  that 
instant.  The  fast  response  of  the  electron's  motion,  as  evidenced  by  the  nonresonant  character  of  Xy  > 
together  with  the  knowledge  that,  at  intensities  of  interest,  the  electron  is  gaining  negligible  energy 
(staying  "cold"),  implies  that  the  absorbed  energy  is  being  passed  on  essentially  instantaneously  to  the 
phonons  via  the  coupling  cq  (2b).  Therefore,  i'  is  valid  to  obtain  the  rate  of  deposition  of  energy  in  the 
lattice  by  determining  the  rate  at  which  the  electrons  absorb  energy  from  the  optical  beam,  and  we 
return  to  this  problem. 

When  a  room  temperature  correction  is  made  to  v  and  w,  we  guess  that,  as  for  small  a,  the 
result  is  raised  by  around  30%  so  that  the  Frohlich  model  may  be  taken  to  yield 


Im  x  ~  20 


(11) 


for  a  ruby  laser  beam  driving  an  electron  of  band  mass  equal  to  a  free  electron  mass  in  ruby  or 
sapphire  at  room  temperature.  The  main  uncertainty  in  cq  (II)  comes  from  the  uncertainty  in  the 
effective  band  mass  m^.  Im  Xy  varies  roughly  as  m£'2. 
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If  the  temperature  goes  higher  than  the  phonon  energy  but  is  still  small  compared  with  t>,  then  eq 
(9)  yields 


Im  x.j  -  4o  «  1  ^2/(3p) 


(12) 


showing  that  the  heated  lattice  causes  the  electrons  to  absorb  even  more  light.  Note  that  the  "Q"  of  the 
electron's  oscillation,  2  Im \v),  is  of  order  7  for  our  estimate  cq  (11),  showitg  that  indeed  the 
optical  frequency  is  much  larger  than  the  effective  electron  collision  frequency  as  we  assumed. 


One  can  see  at  this  point  that  absorption  by  conduction  electrons  may  cause  damage,  for,  if  the 
electron  density  were  10*6  cm*3  at  an  opi*cal  intensity  of  1  0 1 0  W/cm2  (which  the  evidence  outlined  in 
Section  1  suggests  as  a  conservative  estimate),  and  if  cq  (11)  were  valid,  then  cl  ~  10*  3  and  eq  (8) 
gives  —  0.  012  cm**.  For  a  30  nsc  pulse  this  absorption  would  deposit  40J  peF  cc  in  the  lattice, 
somewhat  less  than  the  value  UQ  expected  to  cause  damage.  However,  Belikova,  et  al.  [7],  gave 
convincing  evidence  that  the  electron  density  is  increasing  very  rapidly  with  optical  intensity  at  these 
levels,  an  occurrence  which  we  expect  for  reasons  given  in  the  following  section.  Before  considering 
factors  which  affect  the  electron  density,  however,  we  consider  contributions  other  than  those  of  the 
linear  response  eq  (9)  to  electron  damping. 

First,  the  question  arises  as  to  whether  there  are  significant  intensity-dependent  corrections  to 
the  electron  response  at  the  optical  frequency  as  represented  by  cq  (9).  (Responses  at  harmonic 
frequency  multiples  do  not  contribute  to  the  average  work  done  by  the  optical  field  on  the  electrons  and 
so  we  need  not  consider  these  here.  )  There  are  two  types  of  nonlinear  corrections  to  eq  (9)-  first, 
those  arising  from  the  nonlinear  reactions  of  the  phonons  on  the  electron,  and,  second,  those  arising 
from  the  intensity-dependent  deviations  of  the  momentum  distribution  of  electrons  from  a  thermal 
distribution.  We  have  estimated  both  effects  by  expanding  the  quantum  expression  for  the  expected 
electron  position  (and  velocity)  to  the  third  order  in  the  electric  field.  We  have  found  that  both  effects 
tend  to  diminish  the  optical  absorption  from  its  linear  response  value,  a  result  expected  on  the  physical 
grounds  that  electrons  of  high  energy  interact  less  with  the  optical  phonons  than  do  low-energy 
electrons.  For  both  corrections  the  natural  dimensionless  expansion  parameter  is  (kQ  r0)^  where  kD 
is  the  wavevcctor  of  an  electron  having  energy  h'J  and  rQ  is  the  classical  amplitude  of  a  free  electron 
oscillating  in  the  applied  optical  field.  In  our  case  this  parameter  is  —  I(W/crTi  )  x  10'*  3.  Since 
optical  intensities  have  not  yet  been  observed  to  approach  even  10**  U'/cm*  before  damage,  we  con¬ 
clude  that  the  nonlinear  corrections  to  eq  (9)  arc  too  small  to  be  of  importance  i.-;  present  considerations. 


Next,  we  see  why,  under  conditions  of  interest  here,  phonon  scattering  of  conduction  electrons 
dominates  scattering  by  ionized  and  un- ionized  lattice  impurities.  For  each  conduction  electron  there 
is  in  the  lattice  a  heavy  positive  ion  whose  charge  is  shielded  roughly  by  the  static  dielectric  constant  of 
the  lattice  (e0~  10).  The  linear  absorption  by  10*6  electrons  per  cc  in  such  a  two-component  plasma 
at  room  temperature  is  —  10*6  per  cm  *11].  This  is  so  small  compared  with  phonon  effects  that  one 
can  readily  appreciate  that  even  if  10*^  per  cm3  of  electrically  neutral  scattering  centers  were  added 
to  the  lattice,  one  would  not  expect  them  to  scatter  electrons  as  effectively  as  do  phonons  at  room 
temperature  and  above.  This  predominance  of  phonon  scattering  at  or  above  room  temperature  has 
been  observed  in  various  ways  even  in  crystals  whose  electron-phonon  coupling  is  much  smaller  than  in 
ruby  and  sapphire  '12]  .  Having  estimated  the  optical  response  of  a  single  conduction  electron,  we  now 
proceed  to  consider  how  their  numbers  are  affected  by  their  interactions  with  lattice  impurities. 


5.  Conduction  Electron  Interactions  with  Impurities 


From  incomplete  existing  absorption  [13]  and  luminescence  [8]  data,  it  is  evident  that  sapphire 
has  a  complex  impurity  level  structure  with  energy  differences  ranging  up  to  at  least  3  eV.  Light  of 
any  single  optical  wavelength,  such  as  6943  £.,  will  not  be  resonant  with  a  significant  fraction  of  the 
bound-bound  transitions  and,  of  course,  will  only  excite  transitions  to  the  conduction  band  from  bound 
levels  within  I.  8  eV  of  the  band  edge.  If,  however,  the  light  became  effective  in  exciting  a  large  frac¬ 
tion  of  the  bound-bound  transitions,  then  cascade  transitions  of  electrons  into  the  conduction  band  from 
levels  lying  more  than  1.  8  cV  below  the  band  edge  might  dominate  direct  photo-excitation.  We  shall 
now  show  that,  at  conduction  electron  densities  1  0 1  6  cm*3,  the  light  does  become  so  effective  in 

exciting  bound-bound  transitions  with  which  it  is  not  resonant  that  cascade  photo-excitation  may  become 
important  in  increasing  the  conduction  electron  density  near  optical  damage  intensities.  These  non¬ 
resonant  transitions  are  possible  because  the  conduction  electrons  readily  absorb  the  difference  in 
energy  j  -  <^a  between  an  absorbed  photon  and  the  energy  of  transition  between  the  upper  and  low<  e 
states  |  b>  and  |  a).  We  now  show  that  the  rate  R,k  for  a  typical  impurity  to  make  such  a  transition, 
multiplied  by  a  reasonable  impurity  density  (  ~  1  (JT9  cm*3),  may  exceed  the  total  direct  photo¬ 
excitation  rate,  and  thereby  initiate  cascade  ionization  of  impurities. 

• 

To  calculate  the  rates  Rab  it  is  convenient  to  rc-cxpress  the  interaction  of  eq  (3)  as  a  Fourier 
transform  (the  impurity  index  a  is  now  omitted  since  we  arc  considering  only  one  impurity): 
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where  \k  is  the  impurity  operator  4wc2  V** 


:>  k-2r 


j  expik  ■  rj  and  V  ie  a  fiducial  interaction  volume. 
The  sum  ie  confined  to  the  bound  impurity  electrons  a^s  only” electronic  transitions  will  be  considered. 
For  the  transitions  of  interest  we  may  take  Cs  to  be  the  optical  dielectric  constant  c.  With  the  form  eq 
(13)  it  is  a  straightforward  matter  to  rc-dcrivc  the  quantum  "Golden  Rule”  for  calculating  trai  sition 
rates  but  with  an  unperturbed  Hamiltonian  hQ  =  £(g  -  A0  cosut)2  instead  of  $p  •  To  lowest  order  in  the 
electron-impurity  interaction,  the  rate  for  an  electro?!  to  scatter  from  momentum  state  |j>)  to  |jg  4£> 
while  the  impurity  goes  from  |  a)  to  |b>  is 


-  £  4  £;  *  ”*  b)  s  2"  l<bl  vk|a>|2  ^  «(k2/2  +  £•  k  +  +  nu)  •  [Jn  (k-  rj] 


(M) 


n=-oo 


where  we  have  retained  the  terms  for  n-photon  emission  and  absorption  for  future  comparison  with  the 
n  =  1  term  that  we  arc  now  considering.  In  the  polaron  units  used  here  the  classical  electron  oscilla¬ 
tion  amplitude  rD  equals  K/u2,  where  the  real  amplitude  E  of  the  optical  electric  field  equals  the 
amplitude  in  esu  divided  by  Ec  =  (ug/e)  (ftm^/u0)3.  For  our  parameters  E  equals  (the  field  in  V/ cm)  * 

(1.  24  x  106).  The  Dcssel  function  coefficients  Jn  in  cq  (14)  are  the  exact  amplitudes  for  absorbing  (or 
emitting)  n  optical  photons  in  the  process.  Their  arguments  squared  arc  of  order  1 0~  2  3  i(w/Cm2)  and  so  the 
approximation  J  (x)  ~  (x/2)n/n!  is  accurate  here.  Consider  now  the  rate  rab  for  the  impurity  inter¬ 
acting  with  one  electron  in  the  optical  field  to  make  a  transition  from  a  state  |  a)  to  a  state  |b),  for 
which  <  v.  In  terms  of  the  rate  cq  (14) 


ab 


f(£)  r(£  -  £  +  k:  a  -  b) 


(15) 


where  f(p)  is  the  initial  conduction  electron  distribution  and  only  the  n  =  1  term  in  T  is  needed.  To 
evaluate  raj,  it  remains  to  estimate  the  dependence  of  the  matrix  elements  of  vj^  on  k.  We  shall  use  the 
"dipole  approximation"  in  which  (b|  expik- rj|  a)  ~  ij<  •  (b|  rJ  a).  For  definiteness  we  will  assume  that 
the  impurity  is  spherically  symmetric  so  tliat  |£j  k‘"-”(b|-r'jfa>|  2  =  k2  Xga  where  Xba  the  dipole  matrix 
element  of  the  impurity  transition.  Finally,  before  evaluating  eq  (15),  we  average  over  all  directions 
of  E  so  that  the  Jj2  terms  becomes  (|cE/2b2)2/3.  The  n  1  contribution  to 'the  sum  over  k  in  cq  (1  5)  is 
then  ~ 


where 


rab  ■  V'1  8.  o4(rJ/12)  c*2  X2fe  L(p)  f(p) 


Km4P 

L(P)  =  f  dk  k2n_,/p. 
Km-P 
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Generally  Kjj,  f  [p2  +  2(u  -  Uba)]^  1*  much  larger  than  the  magnitude  of  p  allowed  in  f(p).  In  this  case, 
it  is  sufficiently  accurate  to  approximate  eq  (16b)  by 


L  ~  2[2(v  -  «^)]* 


(17) 


and  we  sec  that  the  impurity  excitation  rate  is  essentially  independent  of  the  initial  small  conduction 
electron  momentum £.  Since  Ep  f(p)  -  1,  this  fact  makes  the  sum  in  eq  (16a)  trivial.  The  total  rate 
Rab  for  an  impurity  to  absorb  a  photon  is  rab  times  the  number  of  conduction  electrons  in  the  interaction 
volume  V-  It  is  convenient  to  express  Rab  in  terms  of  a  photon  absorption  cross-section  oab  which 
equals  Rab  f  (the  incident  photon  flux).  Recalling  that  2b  r2  is  numerically  10*13  I(W/cm2),  we  have 
from  cq  (15)  for  sapphire 

oab~  7-5  x  10-11  I  xjb  e  2  o*1  (b  -  s^)*  Pe  (18) 

in  polaron  units  where  the  electron  density  pe  is  6.  15  x  10  22  times  the  number  of  electrons  per  cc  p. 
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To  •••  th«  ImpoHanc  •  of  tho  tranaition  cron«i*ctloi)  eq  (II)  relative  to  direct  photo-excitation,  i 
we  eetimate  the  deneity  of  impuritiea  Nj  cm"  )  that  would  be  required  to  attenuate  the  optical  beam  more 
than  doee  the  direct  photo* excitation  proceaa.  Of  courae,  the  latter  attenuation  ia  not  known,  but  it 
muat  be  leea  than  the  total  attenuation  in  aapphire.  and  that  ia  leaa  than  10'*  per  cm  For  argument, 
euppoae  that  direct  photo-ekeitation  cauaea  leaa'than  10*)  per  cm  optical  attenuation  (it  probably  cauaaa 
much  teoa).  For  our  eatimate,  let  ua  aaaume  in  eq  (IS)  that  1/^,  —  8.  and  I  *  1,0* 0  W/cm‘ 

and  #  >  10*0  aa  in  previoua  cxamjJlc..  Thi*  givea  —  1 0“  cm*'ao  that  Nj  —  10*^  impuritiea  per  cc 
would  be  required  to  make  theae  conduction-electron- mediated  impurity  tranaitiona  attenuate  the 
optical  beam  by  10**  per  cm.  Conaidering  all  the  kinda  of  dialocationa  and  foreign  iona  which  might 
take  part  aa  impuritiea  in  thia  proceaa,  we  feel  that  it  ia  not  unreasonable  that  euch  ah  impurity 
denaity  could  exiat.  Of  couVae,  there  are  ao  many  ordcr-of-magnitude  eatimatea  of  unknown  quantltiee 
in  the  foregoing  that  the  reault  can  only  be  conaidered  aa  kuggeative  of  a  kind,of  proceaa  that  muat  be  . 
conaidered  aa  competitive  with  other  poaaibilit^ea  at  thia  time.  However,  the  experimental  obaervation 
of  a  rapidly  increaaing  dc  photoconductivity  near  damage  threaholda  [7]  doea  auggeat  a  "bootatrap" 
proceaa  auch  aa  the  above  by  which  conduction  electron^  can  create  more  of  theii*  kind,  even  before 
they  become  ao  energetic  aa  to  ionire  the  lattice. 

We  have  examined  the  enhancement  of  the  aingle-photon  bound-free  photo- excitation  tranaitiona 
that  ariae  from  cq  (14),  and  have  eatimated  that  excitation  croaa-aection  alao  to  be  of  order  10*22  cm* 
under  the  aamc  conditiona.  Thia  ia  probably  leaa  than  the  croaa-aection  for  direct  photo -excitation. 

If  multi-photon  excitation  ia  important  in  producing  eleclrona,  then  the  higher  n  terma  of  eq  (|4) 
probably  dominate  direct  abaorption  near  damage  lovela.  The  expanaion  parameter  in  eq  (14)  ia  of 
order  10*) 1  l|W/cm2),  whereaa  it  ia  —  I0*17  I  for  direct  multi-photon  abaorption. 

1  I 

I  ,  I 

6.  Focualng  Fffecta 

i 

In  a  region  of  the  cryatal)in  which  conduction  electrona  have  been  excited  there  ia  a  change  in  the 
refractive  ihdex  6n  given  by  ,  1 

'  I  ; 

«n  »-«*/. 2 .  1  (19) 

P 

,  t  i 

At  a  denaity  10**  cm’*,  electrona  (with  froe  electron  band  maaa)  hgve  a  plaema  frequency  of  1012  Ht, 
which  givea  uu/v  w,  0.  0023  and  ~  S.  3  x  10*".  Thia  index  changc'ia  comparable  to  that  required  to 
aelf-focua  a  one  milliradian  gauaaian  beam.  If  the  elbetron  denaity  werp  proportional  to  the  optical  i 
intenaity  I,  and  thia  index  change  were  produced  when  I  «  IQ1"  W/em2,  then  the  index  change  would 
conatitute  a  aecond  order  nonlinear  index  (often  written  nj  E*)  whoae  nonlinear  coefficient  n*  would  ba 
-2  x  I0'1]  eau.  Thia  ia  negative  and  an  ordcf  of  magnitude  larger  than  what  one  would  gueaa  ia  the 
intrinaic  nonlinear  index  of  the  bound  electrona  in  aapphire.  However,  we  auapact  that  the  electron  ( 
denaity  ia  increaaing  much  more  rapidly  than  I  at  high  intenaitiea  [7]  and  ao  the  effact  of  the  index  of 
eq  (19)  ia  not  aimplet  Furthermore,  impurity  electrona  which  Ire  excited  to  bound  atatea  more  than  v 
below  the  band  edge,  auch  aa  by  the  proceaa  diacuaaed  in  the  previoua  aection,  have  an  inCreaaed 
polariaability  and  would  tend  to  cancel  the  effect  eq  (19)  of  conduction  electrona. 

Suppoae  the  net  effect  of  impurity  tranaitiona  were  to  produce'a  negative  and  highly  nonlinear  ; 
refractive  index  change,  aa  would  be  the  caae  if  the  conduction  electron  effect  eq  (19)  dominated.  Then 
there  would  be  atrong"aelf-defocuaing  of  a  baam  in  rogiona  of  vary  high  intenaity,  but  only  the  intrinaic 
aclf-focuaing  index  would  be  operative  in  regiona  of  amaller  Intenaity.  Tl^ia  would  evidently' reault  In- 
aucceaaive  regiona  of  high  intenaity  along  the  beanvaa  each  effect  producea  in  turn  a  focualng  tendency 
that  bringa  the  other  effect  into  dominance.  Damage  tracka  in  aapphire  and  other  materiala  conalat  of 
a  aequence  of  damage  bubbiea  that  are  auggeative  of  auch  a  proceaa.  However,  we  have  not  calculated 
the  apacing  of  the  aucceaaive  high-intenaity  regiona  on  the  baaia  of  tome  reaaonable  model  to  aea 
whether  there  ia  any  poaaible  correapondence  between  theae  ,apeculationa  and  obaervation/ : 

•  ‘  i 

7.  Diacuaelon 

Although  there  are  alternate  explanationa  for  the  behavior  of  aamage  tracka  obaervediin  aapphire 
'14],  it  aeema:likely  that  aelf-focuaing  ia  occurring  and  playing  a  role  in  the  damage.  (It  ia  curioua 
that  the  appearance,  lungthe,  and  poaitiona  of  damage  tracka  differ  between  ruby  and  aapphire  [14].)  If 
aclf-focuaing  ia  occurring,  it  makea  eepccially  difficult  the  determination  of  the  actual  optical  intenaitiea 
at  which  the  material  rupturea,  and  thia  in  turn  lenda  an  extra  uncertainty  to  the  parametere  we  hava 
uaed  in  diacuaaing  the  role  of  conduction  electrona  in  damage.  'Probably  the  moat  promiaing  way  to 
determine  more  precisely  what  the  electrona1  role  may  be  ia  to  determine  the  parametere  of  theae 
conduction  electrona  below  damage  intenaitied.  It  would  be  capeciatly  deairablc  to  have  reliable  dc 
photoconductivity  meaaurcmenta  between  the  intenaitiea  $  W/c m*,  the  highcat  uaed  by  llochuli  [5],  apd 
1010  w/c  m*.  the  lowest  used  by  Dclikovat  ct  il.  [7],  in  their  measurements.  With  t|ic  high  light 
intensities  available  with  lasers*  it  may  now  be  possible  to  excite  enough  electrons  to  obtain  a 
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meaningful  llall  effect.  Perhaps  the  time  decay  of  pholoconductivity  and  luminescence  after  an  exciting 
pulse  could  be  observed  and  yield  information  on  carrier  lifetimes,  trapping  times,  etc.  Unfortunately, 
the  important  and  pervading  parameter  —  the  effective  mass  in  the  conduction  band  (and  its  anisotropy)  — 
appears  to  be  extremely  difficult  to  determine  experimentally. 

Should  the  photo-excited  electrons  be  found  to  play  significant  roles  in  damage,  such  as  those  we 
have  discussed,  then  it  may  be  possible  to  fabricate  more  damage-resistant  sapphire  and  ruby.  For 
example,  if  self-focusing  is  presently  instrumental  in  damage,  then  it  mig'it  be  reduced  by  adding 
donor  impurities  that  produce  a  linear  photoconductivity  at  high  intensities  and  thus  reduce  the  total 
nonlinear  index  via  the  index  change  of  eq  (19).  On  the  other  hand,  those  donors  that  give  a  highly 
'nonlinear  photoconductivity  and  a  sudden  production  of  high  electron  densities  that  result  in  rupture  by 
photo-absorption  may  perhaps  be  eliminated  from  the  crystals  or  compensated  by  other  additives.  In 
'any  case,  it  is  evident  that  the  availability  of  high- intensity  laser  sources  will  greatly  facilitate  the 
Study  of  the  conduction  bands  and  impurity  levels  in  very  weakly  photoconducting  materials  such  as 
slpphire  and  ruby. 
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COMMENTS  ON  PAPER  BV  ROBERT  HELLWARTH 

The  medet  employed  in  thi*  papex  m*  ipecHically  developed  (ox  cxyitalline  toUdt  containing 
4 fpuActin.  ThU  model,  involving  definite  impunity  level*  imposed  on  a  well- defined  band  *txuctuxe 
doe*  not  apply  to  giant*  in  it*  pxetent  (oxm.  * 
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APPENDIX  B 


THEORY  OF  ELECTRON  VELOCITY  DISTRIBUTION  IN 

CRYSTAL  IRRADIATED  BY  A  STRONG  OPTICAL  BEAM 

A.  INTRODUCTION 

Described  here  is  the  theoretical  effort  during  the  project  period 
that  was  devoted  to  an  attempt  to  understand  in  a  fundamental  way  the 
heating  of  electrons  in  polar  crystals  by  intense  optical  fields.  Several 
theories  have  been  proposed1 '  c  in  which  the  electrons  gain  enough 
energy  from  the  laser  radiation  to  overcome  the  energy  loss  due  to  the 
radiation  of  optical  phonons  and  are  "heated"  to  energies  exceeding  the 
valence-conduction  band  energy  gap.  These  hot  electrons  ionize  the 
lattice,  each  hot  electron  producing  two  slow  electrons  and  a  hole,  which 
in  turn  are  heated  by  the  optical  field.  This  results  in  an  avalanche  of 
carrier  multiplication  which  results  in  a  rapid  increase  in  the  absorption 
of  laser  energy.  The  absorption  of  a  large  energy  in  a  small  volume  of 
the  crystal  is  assumed  to  be  the  primary  cause  of  the  observed  damage. 

The  theories  of  Zverev,  et  al.  ?  and  Was serman  were  of  the  type 
described  above.  While  appearing  reasonable  at  first  glance,  they 
suffer  from  the  weakness  of  being  based  on  rather  heuristic  arguments. 
Both  theories  predict  an  avalanche  threshold  about  an  order  of  magnitude 
greater  than  that  observed  in  sapphire.  Recent  experiments  at  Massa¬ 
chusetts  Institute  of  Technology  using  CO,  laser  radiation  in  alkali 

£  3 

halide  crystals  seem  to  agree  with  the  estimates  of  Zverev's  theory. 
Hellwarth*  conducted  a  study  in  which  he  found,  on  the  basis  of  a  more 
fundamental  approach,  that  electrons  could  not  be  heated  sufficiently 
by  an  optical  field  in  ruby  or  sapphire  to  trigger  an  avalanche  process. 

He  suggested  some  alternate  absorption  mechanisms. 

Because  of  the  apparent  controversy  concerning  the  possibility 
of  electron  heating  and  avalanche  in  an  intense  optical  field  in  ruby  or 
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sapphire,  we  have  begun  an  independent  study.  We  have  succeeded  in 

deriving  a  kinetic  equation  describing  the  time  evolution  of  the  electron 

momentum  (i.  e.  ,  lattice  momentum)  distribution  function.  The  Frolich 

model  of  electrons  interacting  with  longitudinal  optic  phonons  in  a  polar 

crystal  is  used  as  in  all  the  previous  studies.  The  electron -phonon 

interaction  is  treated  in  lowest  order  perturbation  theory.  This  not  a 

good  approximation  for  ruby  or  sapphire  for  which  the  polar  coupling  is 

4 

strong  and  real  polaron  effects  must  be  taken  into  account.  However, 
it  is  the  basis  for  the  work  of  Zverev^  and  Wasserman*  and  is  hoped 
to  be  of  qualitative  value  for  determining  the  feasibility  of  electron  heat¬ 
ing.  In  the  kinetic  theory,  which  we  are  developing,  the  interaction  of 
the  electromagnetic  field  of  the  laser  with  the  electron  is  treated  exactly 
to  all  orders. 

Unfortunately,  at  the  time  of  this  report,  the  research  is  not 
complete.  The  detailed  consequences  of  this  kinetic  theory  are  still  not 
understood.  Nevertheless,  we  believe  the  kinetic  equation  to  be  described 
here  offers  a  fundamental  basis  to  settle  the  questions  of  electron  heating 
and  work  will  continue  on  this  problem. 

B.  KINETIC  EQUATION 

We  will  present  a  detailed  derivation  of  the  kinetic  equation  for 
electrons  interacting  with  phonons  in  an  intense  ac  electric  field  in 
another  publication.  This  derivation  is  based  on  the  Green's  function 
formulation  of  nonequilibrium  quantum  statistical  mechanics  which  has 
been  used  by  the  author  in  several  previous  studies.  5  Here,  we  will 
start  from  this  kinetic  equation,  carry  out  some  simplifications  and 
discuss  its  physical  content  and  conditions  for  validity.  Let 
n^(t)  dpx  dpy  dpz  be  the  average  number  of  electrons  with  lattice 
momenta  in  the  range  p^  to  px  +  dp^,  etc.  Then  the  equation  for 
the  time  rate  of  change  for  this  momentum  distribution  functior  is  shown 
to  be 
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-00  o  't-t' 

[Vk(T) '  hP(T'  •  “]) 

f1  +  Vk  (‘  ■  i)|  np  (*  ‘  *)  D^a)> 

•  Vk  (‘  *  i)  [l  +  np  (*  •  ?)J  Dk(Ui)  • 


(4) 


Here,  E  is  a  constant  with  the  dimensions  of  an  electric  field  intensity 
c 

which  measures  the  strength  of  the  electron- phonon  coupling  in  the 
Frttlich  theory 


where  cu^q  is  the  LO  phonon  frequency  (assumed  independent  of  k), 

€ _  and  c  are  the  high  frequency  and  static  dielectric  constants, 

00  O 

respectively,  and  m  is  the  electron  effective  mass. 

Approximation  1.  We  emphasise  again  that  this  equation  is  valid 
only  to  lowest  order  in  the  electr<  i- phonon  coupling.  The  LO  phonon 
spectrum  is  contained  in  the  functions 

D*(cd)  =  [no  6(0)  -  u)^)  ♦  (Nc  ♦  1)  6(o>  ♦  o)^)]  (6) 


Dj[o)) 


MO) 


(7) 
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where  Nq  is  the  thermal  equilibrium  phonon  occupation  number 


N  = 
o 


LO  .  , 
e  +  1 


(8) 


Approximation  2.  The  phonons  are  assumed  to  remain  in  equi¬ 
librium  at  the  lattice  temperature  T  =  (kg  =  Boltzmann's 

constant).  This  approximation  will  fail  once  a  considerable  amount  of 
heating  occurs  but  should  be  sufficient  for  the  calculation  of  avalanche 
thresholds. 

The  function  hp(T)  is  the  kinetic  energy  part  of  the  Hamiltonian 
for  an  electron  in  a  spatially  uniform  time-dependent  electric  field 
represented  by  the  vector  potential  A(t) 

h  (r)  =  -L_  (p  -  |  A(t))2  .  (9) 

^  2m 

The  factors 

t 

exp±  i  J  hp(r ) 

t -t ' 

which  make  up  the  cosine  term  in  eq.  (4)  arise  from  the  exact  solution 
of  the  time- dependent  Schrodinger  equation  in  a  spatially  uniform  ac 
field 

ih  ii  *p(t)  =  hp(t)  Vt}  ‘  (10) 

Thus,  the  kinetic  equation  treats  the  motion  of  free  electrons  in  the  ac 
field  exactly. 

Approximation  3.  We  have  assumed  a  parabolic  energy- 
momentum  relation  for  electrons  in  the  conduction  band.  This 
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approximation  is  easily  relaxed  to  take  into  account  nonparabolic ity 
which  is  undoubtedly  important  for  electrons  excited  high  in  the  band. 

For  the  ac  field  represented  by  the  laser,  the  vector  potential 
A(T)  can  be  taken  to  be 


-c  E 

A(t)  =  — — -  s  in  o>o  t  .  (11) 

o 

Approximation  4.  The  electromagnetic  field  is  taken  to  be 
spatially  uniform.  This  is  not  bad  so  long  as  the  magnitude  of  the 
classical  excursion 


d 


e  E 
~  o 


*  2 

m  CD 

o 


(12) 


of  an  electron  in  the  field  is  small  compared  with  a  wavelength  of  the 
actual  physical  field.  This  requirement  is  always  well  satisfied. 
Using  this  form  for  the  vector  potential  we  obtain 


l. 

r  /  dT  [Vk<t>  -  VT,1  =  (V*  -  ep)  r 

t-t' 


k  •  d 


[■ 


Mil 


.  ,  COS  CO  t  -  COS  CD  (t-t')  i  , 

ft  L  o  o  J 


(13) 


2  * 

where  d  is  given  by  eq.  (12)  and  =  p  /2m  . 

Approximation  5.  We  assume  the  electron  density  remains  low 
enough  that  degeneracy  can  be  neglected,  i.  e.  ,  n^(t)  <  <  1.  Then  the 
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expression  in  curly  brackets  in  eq.  (4)  can  be  written  as 

|D<(0»  [np  (t  -  £)  -  Vk<»-.'/2)  «P*"]|  (14) 

where  (7)  has  been  used. 

Next,  we  use  the  well-known  Bessel  function  identity. 


ia  sin  X  \  T  .  .  -iuX 
i  =  )  J  (a)  e 

Z_,  u 

U 


to  write 


cos 


[(■“  -  V* +  O  ?  +  Hr  ( 


, .  .  v,  ,  v-  .  .  ,  cos  co  t  -  cos  a) 

p+k  p/  h  h  \  o  o 


(.-«•))] 


i(  co-e  ,,+£  )t/ti  r-  ,v  ,, 

- Re  ptk  p  I  ^  (V)  \  (V) 


»v 


pw 


-  lfX  CX>o(t -t  )  ^i(u  J 


(15) 


(16) 


To  consider  the  periodicity  impressed  upon  the  system  by  the 
strong  ac  field,  it  is  useful  to  write  the  distribution  function  as  a  Fourier 
series  with  fundamental  frequency  coq,  whose  coefficients  are  slowly 
varying  functions  of  time  on  the  scale  of  co~*. 

—  -in co  t 

np(t)  =  1  Fn(p>t)  e  ■  (17) 

n 

On  substituting  eqs.  (17),  (161,  (13),  and  (14)  into  eq.  (4),  we  obtain  the 
coupled  set  of  equations  for  the  Fourier  coefficients 
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jji.m 


[Fm(p+k,t)  ep,K0-  Fm(p.«)]  • 


(18) 


Approximation  6  (Weak  Coupling).  This  is  an  extremely  com¬ 
plicated  infinite  hierarchy  of  equations  in  general.  However,  if  we 
assume  that  the  electron-phonon  interaction  as  measured  by  Eq  is 
weak,  then  the  assumption 

Fq  >>  Fn  for  n  4  0  (19) 

is  valid.  This  follows  since  9F  /dt<<  noo  F  for  n  i  0.  Then  if  we 

n  on 

adopt  an  iterative  procedure  starting  by  keeping  only  the  m  =  0  term 
on  the  right  side  of  eq.  (18),  we  see  that  Fn  is  proportional  to  Ec/n03Q 
times  Fq  (n  4  0)  which  we  assume  to  be  small.  We  repeat  that  this 
procedure  is  valid  only  in  the  weak  coupling  limit  for  the  electron- 
phonon  approximation.  In  fact,  the  kinetic  equation  (18)  is  only  valid 
in  this  case.  This  approximation  fails  for  the  parameters  of  ruby  and 
sapphire.  We  use  it  mainly  because  it  is  the  only  known  procedure  for 
the  nonequlibrium  polar  problem. 

Thus,  to  lowest  order  in  the  electron-phonon  coupling,  the 
fundamental  kinetic  equation  for  the  dc  component  of  the  electron 
distribution  function  is  taken  to  be 
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dFo(p,t)  eEc  c  d3 
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2irm 


CD 


j  j  dmD<<“> 


-CD 


I  jS(V)6('i“+»*“o- Vk+s) 
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[Fo(p+k,t)  e?*®  -  Fo(p,t)] 


(20) 


We  have  not  yet  succeeded  in  obtaining  a  complete  solution  of 
this  equation.  Not  the  least  complication  is  the  involved  dependence  on 
the  radiation  intensity  occurring  through  the  Bessel  functions  J  (k  •  d/ft). 
Examination  of  the  requirements  imposed  on  k  by  the  delta  function 
shows  that  values  of  k  as  large  as  /2m*u u)q  can  be  important  in  the 
integrand.  In  intensity  units  then,  (k  .  d/ti)2  <  v  1 0“ 1 3  I  ,  where  I  is  in 
watts/cm2.  For  the  laser  powers  necessary  for  damage  (I  ~  10^®  W/cm2), 
this  parameter  is  small  and  the  Taylor  expansions  of  the  Bessel  functions 
can  be  used. 
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Vk+ep) 
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+  6  (*oi  -  m>o  .  +  ep)] 


[F0(p+k,t) 


,(P.')] 


(21) 
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In  this  limit  the  various  terms  in  the  equation  have  a  simple  interpre¬ 
tation  in  terms  of  the  transition  rates  for  elementary  processes 
calculated  by  Fermi's  Golden  Rule.  The  first  term  in  curly  brackets 
of  order  1  can  be  shown  to  produce  the  difference  of  the  rates  for 
scattering-in  and  scattering-out  of  the  state  p  by  emission  or  absorption 
of  LO  phonons.  The  first  term  of  order  (k«  d)/ft))^  is  the  correction 
to  these  rates  for  the  coherent  possibility  of  absorbing  and  emitting  a 
photon  in  this  process.  The  second  term  of  order  (k*d)/ft)  produces 
the  difference  of  the  rates  of  scattering-in  and  scattering-out  of  the 
state  p  by  either  emitting  or  absorbing  a  photon  and  emitting  or 
absorbing  an  LO  phonon.  The  latter  processes  are  essentially  those 
of  electron- Dho non  Bremstrahlune  emission  and  absorption. 

The  kinetic  equation  (20)  is  not  complete  for  this  problem  be¬ 
cause  it  does  not  contain  terms  which  account  for  the  process  of  impact 
ionization.  We  will  not  try  to  formulate  these  terms  precisely  here, 
but  rather  will  discuss  their  qualitative  features.  The  missing  terms 
consist  of  a  high-energy  "sink"  or  loss  term  which  accounts  for  those 
electrons  whose  energy  exceeds  the  ionization  threshold  Ej  >  ^gap 
and  lose  their  energy  by  impact  ionization  collisions  with  the  lattice. 
There  is  also  a  low-energy  "source"  term  which  takes  into  account 
the  two  low-energy  electrons  which  appear  as  a  result  of  the  impact 
ionization  processes  (i.  e.  ,  the  high-energy  electron  loses  its  energy 
to  produce  another  low-energy  electron  plus  a  hole).  We  may  consider 
that  eq.  (20)  is  valid  for  an  intermediate  band  of  energies  ^max  <  C  <  £j 

where  £  is  a  measure  of  the  maximum  energy  at  which  new 

max  67 

electrons  are  created.  It  appears  reasonable  to  take  £  <fcco  . 

max  o 

C.  APPROXIMATE  CONSIDERATIONS 

In  an  effort  to  gain  some  feeling  for  the  predictions  of  eq.  (20) 
we  have  taken  the  case  k  •  d  <  <  1,  which  is  well  satisfied  by  the  experi¬ 
ments.  As  we  have  seen,  the  Bessel  functions  can  be  expanded  in  this 

case  as  in  eq.  (21).  It  is  convenient  to  transform  this  equation  to  new 

2  .  * 

variables,  £  =  p  /2m  and  p  —  the  cosine  of  the  angle  between  the 
radiation  polarization  E  and  jp. 


113 


dFo(e,n,t)  e  Ec 
5t  p 


J  dfi'  J  d cd  D* 

-  1  -  CD 


(CO) 


where 


1  -  <l/2ft2)  d^pfi-p^')2 


oo 


d2  (pH  -  pn')2 

—  - ^ —  me 

4*T  ±  L 


[F(e+oo,p')  e*3^  -  F(e,n)J 


+  OO  ±  00  ,  p')  e 


|3ftoo 


-  F(e,fi)"| 


p  =  V2m 


(22) 


raj 


=  y/Zrn'  f 


(e+co) 


p  =  '4m'(£+oo±oo  ) 
*±  o 


(23) 
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(24) 


t  =  -L  /_£_  +  £i£\  t  =  -  (-2-  + 

“  2  VPco  P  /  ±  2  \p  p  ) 


(25) 


The  last  term  in  eq.  (22)  is  a  shorthand  notation  for  the  sum  of  both 
the  +  and  -  indexed  quantities. 

Next  we  write  FQ(£,p)  =  and  make  use  of  the  fact 

that  X  will  be  a  more  slowly  varying  function  of  £,p  than  Fq.  We 
divide  the  equation  through  by  Fq.  Differences  such  as 


X(£  +  oo±ooo,p')  -  X(£,n) 
are  Taylor-expanded  in  the  form 

(oo±ooo)  xe(e,p)  +  (p'-n)  x^e.p) 
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and  higher  derivative  terms  dropped.  Here  X_  =  3X/5£  and  X  =  dX/dp. 

c  f* 

The  resulting  equation  can  be  written 
oo 
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where  we  have  assumed.  u>  >  >  co(=  ±  coT  0)  and  where  we  define 
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(27) 


(28) 


In  the  region  £^  >  £  >>  goq  >>  co^q  it  is  clear  that  the  parameters 
t  defined  in  eq.  (25)  all  approach  unity  so  that  A_tuT  ~  -•  A+qw  and 
B+  -•  B  .  In  this  case  the  integration  over  co  (see  eq.  (6))  is  trivial  and 
the  first  factor  in  square  brackets  on  the  right  side  of  eq.  (22)  can  be 
written  as 

M  A  f(  -PW1*  +  P^LO  (  P^V1*  ,)] 

No  ^LO  Xp^  |le  -  lj  +  e  \e  -  1/J 


+  (2Nq+  1)  [Aojlq^p’.X^)  -  A^^fi.p.O)] 


(29) 


sNo  AoLO(e^’Xp)  (P“LO»-  X£<V1,+  2Nc 


[aoolq^P.X^)  -  AajLO(£,p,  0)].  (30) 


115 


In  the  last  member  of  this  equation  we  have  assumed  pu)  LQ«  1  and 
expanded.  This  is  a  good  approximation  for  crystals  above  room 
temperature.  In  this  limit  NQ  *  (po)^)’*.  In  the  second  term  in  square 
brackets  on  the  right  side  of  eq.  (22),  we  note  that  since  pco^  a;  80  for 
the  cases  of  interest,  the  term  with  the  positive  exponential  is  likely  to 
dominate  so  we  neglect  the  other  terms. 

Finally  we  note  that  we  are  expecting  exponential  growth  in  the 
number  of  carriers  during  avalanche.  Thus  we  assume  that 

F0<P.‘>  =  Fo(P»  e*/T  •  <3I) 

where  t  is  a  characteristic  growth  time.  Then  dFo/ot  =  F  /t.  The 
growth  is  a  result  of  the  high-energy  "sink"  for  £  >  £j  and  the  low- 
energy  "source"for  £  <u)q  which  arise  from  the  impact  ionization 
process.  These  terms  do  not  manifest  themselves  explicitly  in  the 
kinetic  equation  in  the  region  u>o  <£  <£j  except  through  the  time 
derivative  term.  Putting  all  these  assumptions  and  approximations 
together  leads  to  the  equation 

(P^LO)p  2 

tTT-  =  a-lo  <«•*•  V  xe<xc-l) 


r-  -  -.2  2  P^ 

+  2  •  A<e^’°)]  4  2"  C  ° 

(32) 

We  have  thus  reduced  cq.  (22),  which  is  an  integro-differcntial- 
difference  equation,  to  a  nonlinear  partial  differential  equation.  Unfor¬ 
tunately,  wc  also  have  not  been  able  to  find  completely  satisfactory 
analytic  approximations  to  the  solution  of  this  equation.  It  is  probably 
solvaolc  by  computer  methods,  but  wc  have  had  neither  time  or  support 
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for  this.  The  boundary  conditions  for  {his  equation  must  be  determined 

t 

by  examining  the  physic#  of  the  high-  and  low-energy  sink  and  source 
terms  representing  the  impact  ionization  process.  This  problem  is 
currently  tinder  ‘study. 

One  very  approximate  conclusion  can  be  obtained  from  eq.  (32). 
First,  we  note  that  for  zero  radiation  intensity  and  steady  state,  i.  e.  , 
d  =  0  and  t  =  oo,  this  equation  has  the  solution  X,  =  1,  X  =  0, 
corresponding  to  a  Maxwellian  distribution  at  the  lattice  temperature. 
Now,  in  the  presence  of  the  radiation  the  electrons  will  be  "heated"  most 

•  1  I  i 

strongly  in  the  p  =  ±  1  directions  (this  follows  from  the  dependence  of 

i  i 

the  field  arising  as  k  •  d).  Thus,  we  might  expect  the  tail  of  the  distri¬ 
bution  to  be  drawn  out,  i.  e.  ,  Xg  <  1  in  the  r  t  1  directions,  bpt  to 
fall  off  faster  in  other  directions.  .This  implies  pX  <  <  -1  for  p  =  +1. 

i  I*.  i 

From  eqs.  (27)  and  (28),  it  can  be  shown  that  for  p  fi:  1  and 

li»Vl>>i.  ■  ,  , 


AojlqIC,  !;Kp)  =  -y  -  log 


where  y  =  0.  5772  (Euler's  constant). 


B  (fc.l.jt  )  =  — I—*  .  '  (34) 

*  “  <px/ 

•  ,  .  I  I 

To  a  good  approximation,  an  asymptotic  relation  for  X.  can  be  obtained 

|  ^ 

from  eq.  (32). 

1  '  (PV2  13<“  fx  '  h(«)  '  'll 

xe  =  i>ra-  *»«  ttJ-*  '?  r  {-*£-  t»Y4i.g  |2px|1|)j| 


where  t  .(C)  =  Jim* l  /  N  E  is  the  well-known  collision  time  for 
e  •  pn  i  o  c  i 

electron-polar -optic -phonon  collision.  This  relation  is  valid  as  long  as 
the  argument  of  the  log  is  >  >  1.  If  |pX  |  >>  1,  we  see  that  Xp  is  a 
decreasing  function  of  the  radiation  intensity.  In  fact,  since 

&  80,  Xc<  <  1  for  any  reasonable  assignment  of  parameters.  Note 
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that  for  the  term  involving  the  e-folding  time  r  to  be  important  require* 

using  m*  =  m  and  E  s  2  x  10^  e.  *.  u. 
c  c 
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t<<t  .  or  t  <  <  10  *ec 
e-pn 


for  sapphire  and  l  = 

If  r  >  >  t  .(C).  we  can  neglect  the  T-dependent  term  which 
e  -  ph 

implies  that  the  steady  solution  of  the  kinetic  equation  (i.  e.  ,  eq.  (4) 
with  dF/dt  =  0)  is  all  we  need  to  consider.  In  this  case,  the  solution 
is  actually  independent  of  the  magnitude  of  Ec,  as  long  as  it  is  small 
enough  so  that  the  weak  coupling  approximation  which  we  have  been  using 
is  valid.  The  determination  of  T  as  a  function  of  IQ  awaits  the  com¬ 
pletion  of  the  solution  satisfying  the  boundary  conditions  imposed  by  the 
impact  ionization  process. 

This  very  rough  estimate  of  the  behavior  of  the  solution  of  the 
kinetic  equation  makes  plausible  the  following  model  of  the  distribution 
function:  The  population  of  electrons  is  highly  nonisotropic,  being  drawn 
out  into  high-energy  tails  in  the  (|*  =  ±  1)  direction  of  the  radiation 
polarization.  These  high-energy  tails  may  be  sufficient  to  trigger  the 
strong  avalanche  of  impact  ionization  events  necessary  for  the  Wasserman- 
Zverev  model.  However,  considerably  more  analysis  of  the  detailed 
solution  of  the  kinetic  equation  is  necessary  to  establish  this  mode. 

Note  that  the  average  kinetic  energy  in  the  distribution  will  be  much 
lower  than  the  ionization  energy  Cj  but  there  may  still  be  a  sufficient 
number  of  hot  electrons  to  trigger  the  avalanche.  This  is  consistent 
with  Hellwarth's  observation  that  the  mean  kinetic  energy  of  electrons 
is  much  lower  than  the  band  gap  energy. 
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